LIBRARY 

NAVAL  POSTGRADUATE  SCHOOL 

MONTEREY,  CALIF.  93940 


TH7657 


POROSITY  DETERMINATION  BY  ELECTRICAL  MEANS 


DOCTOR  OF  PHILOSOPHY  THESIS 


OF 


RONALD  ANTON  ERCHUL 


LIBRARY 

NAVAL  P    IRADUATE  SCHOOL 

MONTEREY,  CALIF.  93940 


THE  USE  OF  ELECTRICAL  RESISTIVITY  TO  DETERMINE 


POROSITY  OF  MARINE  SEDIMENTS 


BY 

RONALD  ANTON  ERCHUL 
II 


A  THESIS  SUBMITTED  IN  PARTIAL  FULFILLMENT  OF  THE 

REQUIREMENTS  FOR  THE  DEGREE  OF 

DOCTOR  OF  PHILOSOPHY 

IN 

OCEAN  ENGINEERING 


UNIVERSITY  OF  RHODE  ISLAND 
1972 


£"5twi- 


ABSTRACT 

A  method  using  electrical  resistivity  measurements  to  deter- 
mine the  in  situ  porosity  of  marine  sediments  was  investigated  in 
the  laboratory,  and  equipment  for  this  purpose  was  designed,  fabri- 
cated and  tested.   Formation  Factor-porosity  relationships  deter- 
mined in  the  laboratory  for  three  clays  (kaolinite,  illite,  and 
montmorillonite) ,  Providence  silt,  four  sands,  and  four  marine  sedi- 
ments showed  that  porosity  was  predicted  within  +  2  percent.   The 
Formation  Factors  ranged  from  1,1  to  5,9  while  porosity  ranged  from 
26  to  93  percent.   The  particle  size  and  distribution  influenced  the 
electrical  resistivity  of  these  sediments  independent  of  porosity 
while  particle  shape  did  not.   The  laboratory  equipment  was  econo- 
mical, safe,  easy  to  operate  and  could  be  used  to  determine  perme- 
ability, tortuosity  and  void  ratio-log  consolidation  pressure  of 
sediments. 

The  Electrical  Resistivity  Measuring  System  designed  for  the 
Deep  Ocean  Sediment  Probe  consists  of  three  interchangeable  electrode 
arrays,  the  electronic  circuit  and  the  FM  telemetry  data  link  with 
6000  feet  of  coaxial  cable.   The  predicted  porosity  values  obtained 
with  the  System  in  the  laboratory  were  within  +  2  percent  of  the  line 
of  best  fit  obtained  with  the  earlier  laboratory  equipment.   Using 
the  Inner  corer  ring  electrode  array,  the  System  was  tested  in  the 
shallow  water  of  Narragansett  Day  and  most  of  the  data  agreed  well 


with  the  value  obtained  later  with  the  laboratory  equipment  on  the 
same  sediment. 

The  inner  corer  ring  electrode  array  read  continuous  conducti- 
vity through  the  water/sediment  column  thus  permitting  a  value  of 
salinity  to  be  obtained.   The  bottom  water  salinity  obtained  in  this 
way  was  in  good  agreement  with  the  interstitial  water  salinity  of  the 
cores.   The  system  detected  a  core  loss,  touch  down,  and  sample  dis- 
turbance.  This  relatively  low  cost  system  allows  accurate,  quick 
in  situ  porosity  determination  and  shows  promise  in  other  areas  such 
as  the  monitoring  of  pollution,  the  prediction  of  other  bottom  sedi- 
ment properties  and  as  a  warning  system  on  ocean  structures. 
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PREFACE 

This  dissertation  is  organized  according  to  the  standard  plan. 
The  text  is  divided  into  two  sections.   The  first  section  includes  a 
report  of  laboratory  experiments  on  an  indirect  method  of  determining 
porosity  of  sediments  and  a  review  of  the  literature  involving  similar 
work.   The  instrumentation  used  in  the  laboratory  is  described  in 
this  section. 

The  second  section  describes  the  design,  fabrication,  testing, 
and  evaluation  of  an  in  situ  system  used  to  determine  porosity  of 
marine  sediments  and  the  equipment  used  by  others  for  similar  work. 
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THE  USE  OF  ELECTRICAL  RESISTIVITY  TO  PREDICT 
THE  POROSITY  OF  SATURATED  SEDIMENTS 


INTRODUCTION 

The  purpose  of  this  work  is  to  investigate  the  accuracy  of 
using  an  electrical  resistivity  method  to  predict  the  porosity  of 
soil-water  mixtures.   Porosity  in  a  two  phase  soil  water  mixture  is 
the  ratio  between  the  volume  of  water  to  the  total  volume.   The  deter- 
mination of  porosity  is  important  to  ocean  engineering  projects  such 
as  underwater  acoustics,  foundation  engineering,  submarine  soil  me- 
chanics and  bottom  survey  work.   At  present,  porosity  is  determined 
in  the  laboratory.   If  electrical  resistivity  is  found  to  be  an  accu- 
rate method  of  predicting  porosity,  a  practical  method  of  in  situ 
testing  would  be  available.   This  method  would  be  quicker,  more  eco- 
nomical, and  allow  more  data  to  be  taken  than  the  present  method. 

To  verify  the  accuracy  of  the  electrical  resistivity  method, 
extensive  laboratory  tests  were  conducted  and  the  results  were  com- 
pared with  theoretical  and  empirical  work  reported  in  the  literature. 
Most  of  the  earlier  work  studying  the  relationship  between  electrical 
resistivity  and  porosity  was  conducted  on  formation  rocks  or  disper- 
sions.  The  geophysical  techniques  used  in  well  logging  analysis 
involve  electrical  resistivity  measurements  in  order  to  obtain: 

1.  Subsurface  correlation  with  formations  in  other  wells, 
structural  mapping. 

2.  Lithology  of  the  formations. 

3.  Depth  and  thickness  of  productive  zones. 


4.  Formation  porosities. 

5.  Oil  and  gas  saturation. 

6.  Volume  of  oil  in  a  reservoir. 

The  porosity  of  formation  rocks  ranges  from  5  to  30  percent, 
while  the  porosity  of  emulsions  and  dispersions  ranges  from  50  to 
90  percent.   Although  these  ranges  cover  the  porosity  values  of 
sediments  (which  range  from  26  to  90  percent),  extensive  testing  is 
necessary  before  it  can  be  considered  valid  to  apply  the  models  and 
theories  used  for  rocks,  emulsions,  and  dispersions  to  sediments. 

Electrical  resistivity  has  also  been  used  in  order  to  deter- 
mine the  dielectric  constant  of  pure  crystals,  the  volume  concentra- 
tion of  certain  biological  materials  such  as  red  blood  cells,  and 

1  2 

ground  water  detection.  (Meredith  ,  Patten  and  Bennett  ) 

Very  little  electrical  resistivity  work  has  been  conducted  on 
sediments.   The  limited  work  that  has  been  conducted  on  sediments 
generally  has  not  considered  clay  mineral  content  and  type,  particle 
size,  shape,  distribution  and  structure,  and  consolidation  pressure. 

In  this  study  the  following  parameters  were  varied  in  order  to 
measure  their  influence  on  the  prediction  of  porosity  as  determined 
by  electrical  resistivity  measurements: 

1.  Type  of  sediment  (sands,  silt,  kaolinite,  illite,  and 
montmorillonite  clays,  and  marine  sediments), 

2.  Size  of  particles, 

3.  Shape  of  particles, 

4.  Distribution  of  particle  size, 

5.  Structure  of  particles, 


6,  Consolidation  pressure, 

7.  Salinity  of  interstitial  water. 

The  temperature  was  measured  and  a  correction  factor  was  applied  in 
order  to  standardize  all  electrical  resistivity  measurements  to 
25  Centigrade.   By  measuring  the  electrical  resistivity  as  each 
parameter  was  varied,  while  at  the  same  time  determining  porosity  by 
the  standard  method,  the  influence  of  each  parameter  on  the  accuracy 
of  this  method  of  porosity  prediction  could  be  evaluated. 

Since  the  terms  found  in  this  work  are  used  both  in  submarine 
soil  mechanics  and  geophysics,  they  are  defined  below  in  order  to 
eliminate  confusion. 

The  two-phase  mixture  of  soil  and  water  is  considered  a  satu- 
rated sediment,  or  simply  a  sediment.   The  electrical  resistivity 
measurements  are  taken  on  both  the  sediment  and  the  interstitial 
water.   The  resistivity  of  the  sediment  is  divided  by  the  resistivity 
of  the  interstitial  water  to  obtain  a  ratio  called  the  Formation 
Factor.   This  ratio  is  used  as  a  direct  index  of  porosity. 

Electrical  resistivity  (specific  resistance)  is  a  measure  of 
the  ability  of  a  solution,  solid  or  mixture  of  unit  length  and  cross 
section  to  resist  the  flow  of  electric  current.   In  solids,  current 
is  carried  by  free  electrons;  In  solutions,  the  current  is  carried  by 
ions.   As  electrons  and  ions  move  through  solids  and  solutions,  a 
frictional  drag  results.   This  frlctlonal  drag  per  unit  volume  is 
considered  the  electrical  resistivity  of  the  solid  or  solution. 
Electrical  conductivity  (specific  conductance)  is  a  measure  of  the 


ability  of  a  solution,  solid,  or  mixture,  of  unit  length  and  cross 
section  to  carry  an  electrical  current.   Consequently  resistivity 
and  conductivity  are  inversely  related  measures  of  the  same  property, 
Studies  dealing  with  formation  rocks  have  related  resistivity  to 
porosity  while  studies  dealing  with  emulsions  and  dispersions  have 
usually  related  conductivity  to  porosity. 

Listed  below  are  some  of  the  symbols,  assumptions,  and  defini- 
tions used  in  this  paper: 

dispersed  phase  (d)  -  solid  soil  particle, 

continuous  phase  (c)  -  interstitial  water, 

mixture  (m)  -  the  mixture  of  the  continuous  and  dispersed 
phases  -  sediment. 
The  properties  of  these  phases  are  assumed  to  be  isotropic  and  uni- 
form. 

r   ,   -  resistivity  of  the  continuous  phase  dispersed 

C  i  0  •lu 

phase  and  mixture  respectively. 

k   ,   -  conductivity  of  the  continuous  phases  dispersed 
c  •  o  *  m 

phase  and  mixture  respectively. 
In  order  to  make  relationships  dimensionless  the  resistivity 
and  the  conductivity  of  the  continuous  phase  is  considered  unity. 
These  dimensionless  notations  are  denoted  by  Roman  capital  letters: 
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R  -  —  •  1   ,  (1) 

c   r 
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Rd  =  r    -  '  <2> 
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R  -  -  ,  (3) 

m   r   ' 
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k 

m   k 
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Formation  Factor  (FF)  which  is  often  used  in  the  Geophysical 

literature  is  equal  to  R  and  the  reciprocal  of  K  : 
1        m  r  m 


FF  =  R  =  —  .  (7) 

m   K 
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Porosity  is  normally  defined  as  the  volume  of  the  voids  to  the  total 
volume.   Porosity  in  a  saturated  sediment  (two  phase  system)  can  be 
defined  as  the  ratio  of  the  volume  of  water  to  the  total  volume: 


n  /   \      VOldS  /ov 

Porosity  (n)  = ,  (8) 

total 


if  saturated 


water  ,     *,nS 

n  = .  (9) 

total 


REVIEW  OF  THE  LITERATURE 

The  literature  contains  many  theoretical  and  empirical  models 
that  were  developed  in  order  to  predict  porosity  from  electrical 
measurement.   The  following  review  examines  the  models  that  might 
logically  be  applied  to  a  saturated  sediment.   No  one  model  can 
cover  the  wide  range  of  sediment  types  and  mixtures  (sands,  silts 
and  clays).   Factors  such  as  particle  size,  shape,  distribution  and 
structure,  and  consolidation  pressure  influence  model  selection. 
In  addition  conduction  of  the  clay  fraction  changes  as  the  ionic 
concentration  of  the  interstitial  water  varies.   The  first  part  of 
this  review  lists  some  of  the  properties  peculiar  to  sediments  that 
need  to  be  considered  in  models  in  order  to  enable  accurate  porosity 
prediction  from  electrical  measurement. 

Except  for  native  metals,  sulfides,  a  few  oxides,  graphite, 
and  high-grade  coal,  the  resistivity  of  most  minerals  that  make  up 
sediments  is  measured  in  millions  of  ohm-meters.   Although  the  soil 
mineral  is  non-conductive,  the  soil  particle  carries  an  electrical 
charge.   The  net  electrical  charge  may  arise  from  any  one  or  a  combi- 
nation of  the  following  factors: 

1,  Isomorphous  substitution  (important  in  illite  and 
montmorillonite  clays), 

2,  Broken  bond  charges  (important  in  kaolinite  clays), 

3,  Exposed  hydroxyl  groups, 
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4.  Absence  of  cations  in  the  crystal  lattice, 

5.  Presence  of  organic  matter. 

The  cation  exchange  capacity  is  a  measure  of  the  electrical 
charge  carried  by  a  soil  particle.   The  magnitude  of  the  cation  ex- 
change capacity  of  a  sediment  depends  on  the  amount  and  type  of  clay 
fraction  present  and  the  clay  particle  surface  area*   The  clay  frac- 
tion is  usually  defined  as  that  part  of  the  soil  mass  which  has  an 
equivalent  Stokes  diameter  of  less  than  two  microns.   Particles 
smaller  than  one  micron  are  considered  in  the  colloidal  range.   It 
is  in  the  colloidal  range  that  surface  properties  begin  to  dominate 
the  physical  and  chemical  properties  of  materials. 

For  most  clay  particles  the  specific  surface  (the  magnitude 
of  the  surface  area  per  unit  mass)  is  a  good  indicator  of  the  rela- 
tive influence  of  electric  forces  on  the  behavior  of  the  particle. 
Table  1  shows  the  cation  exchange  capacity  and  specific  surfaces  of 
the  common  clay  minerals. 


TABLE  I 

THE  CATION  EXCHANGE  CAPACITY  AND  SPECIFIC  SURFACE 
OF  KAOLIN ITE,  ILLITE,  AND  MONTMORILLONITE  CLAYS 


CLAY  MINERAL 

,meq.  . 

vioo  g; 

2 
(2L) 

g 

CATION  EXCHANGE  CAPACITY*' 

SPECIFIC  SURFACE^ 

Kaolini te 

lllite 

Montmorillonite 

3  -  15 
10-40 
80  -  150 

10 

100 

1000 

*  3 
Wu 


**  4 

Lambe  and  Whitman 


When  individual  clay  particles  are  dropped  into  water  the  ions  on 
the  surface  of  the  particle  tend  to  move  out  from  the  surface  to  be- 
come exchangeable  ions.  Water  molecules  are  attracted  to  these  ions 
and  the  two  together  form  what  is  called  the  double  layer  around  the 
clay  particle.  The  thickness  of  the  double  layer  is  defined  as  the 
distance  from  the  soil  surface  to  the  point  where  the  ion  concentra- 
tion is  in  equilibrium  with  the  interstitial  water. 

Many  factors  influence  the  size  of  the  double  layer.   The 
Gouy-Chapman  Double  Layer  theory  is  often  used  to  explain  the  double 
layer  and  the  factors  affecting  it.   This  theory  assumes  that  there 
is  no  interaction  between  clay  particles  and  that  the  clay  particles 
are  plates.   This  theory  states  that  the  size  of  the  double  layer  is 
directly  proportional  to  the  temperature  and  dielectric  constant  of 
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the  solution  and  inversely  proportional  to  the  concentration  and 
ionic  valence  of  the  electrolyte  in  the  interstitial  water  and  the 
electric  charge  on  the  clay  particle.   Other  factors  influencing 
the  double  layer  are  pll  and  size  of  ions  in  solution. 

If  the  cation  exchange  capacity  and  the  specific  surface  are 
low  and  a  small  double  layer  exists,  almost  all  the  interstitial 
water  would  be  free  water  and  would  not  be  attracted  to  the  mineral 
skeleton.   Conversely  if  the  cation  exchange  capacity  and  the  specif- 
ic surface  are  high  and  a  large  double  layer  exists,  it  would  be 
quite  difficult  to  separate  the  clay  mineral  and  interstitial  water 
phases. 

One  of  the  earliest  works  examining  the  effect  of  charged 
clay  particles  and  the  associated  double  layer  on  electrical  resis- 
tivity was  done  on  reservoir  rocks  by  Winsauer  et  al.   These  in- 
vestigators concluded  that  charged  clay  particles  do  influence  the 
conductivity  of  shaly  sand  and  limestone  formations,  but  that  this 
influence  became  negligible  at  high  ionic  concentrations  of  inter- 
stitial water. 

After  conducting  laboratory  experiments,  Hill  and  Milburn 
made  similar  conclusions  and  developed  empirical  equations  for 
reservoir  rocks  relating  Formation  Factor  to  the  resistivity  of  the 
interstitial  water,  porosity  and  the  cation  exchange  capacity.   In 
their  equations  as  the  electrical  resistivity  of  the  interstitial 
water  decreased,  i.e.,  became  more  saline,  the  effect  of  the  cation 
exchange  capacity  decreased. 

One  of  the  first  works  to  study  the  effect  of  the  cation  ex- 
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change  capacity  on  the  electrical  resistivity  of  a  sediment  was  con- 
ducted by  Berg.    Using  kaolinite  clay  and  sodium  chloride  solutions, 
he  showed  that  the  Formation  Factor  remained  constant  when  the  water 
resistivity  was  less  than  80  ohm-cm.   However  when  water  resistivity 
was  greater  than  80  ohm-cm,  the  Formation  Factor  decreased.   Water 
that  has  an  electrical  resistivity  of  80  ohm-cm  would  have  a  salin- 
ity of  seven  parts  per  thousand.   Sea  water  having  a  salinity  of 
35  parts  per  thousand  would  have  a  resistivity  of  25  ohm-cm.   This 
work  agrees  with  the  work  on  reservoir  rocks.  '   The  effect  of  the 
charge  of  kaolinite  clay  particles  on  the  electrical  resistivity  ap- 
peared to  be  negligible  at  high  ionic  concentrations  of  interstitial 
water. 

o 

Sheeler  et  al.   expanded  the  knowledge  on  sediments  by  study- 
ing several  types  of  clay.   He  found  that  for  kaolinite  clay  the  ef- 
fect of  cation  exchange  capacity  on  the  electrical  conductivity  and 
thus  Formation  Factor  persisted  up  to  a  3  percent  sodium  chloride 
solution,  then  became  negligible.   However  when  using  montmoril lonite 
the  ion  effect  persisted  up  to  an  8  percent  sodium  chloride  solution. 
Sea  water  is  approximately  3.5  percent  sodium  chloride.   This  study 
suggests  that  the  type  of  clay  present  in  a  sediment  such  as 
montmori lloni te  may  influence  the  electrical  resistivity  measurements 
due  to  the  cation  exchange  capacity  even  though  the  salinity  is 

higher  than  sea  water. 

9 
Doyce   reported  that  the  percentage  of  clay  size  material 

varied  inversely  with  the  electrical  resistivity  of  fine  grained 

marine  sediments  obtained  from  the  Bering  Sea.   Boyce  explained  these 
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results  could  be  due  to: 

1.  a  specific  conductivity  increase  resulting  from  a 
double  ionic  layer  of  clay  mineral  in  the  clay  size 
fraction. 

2.  a  secondary  conductivity  increase  caused  by  an  increase 
in  porosity,  since  the  clay  size  fraction  and  porosity 
are  directly  related. 

Fritsch  and  Tauber   studied  kaolinite,  illite,  and 
montmorillonite  clays  at  low  and  high  salt  concentrations.   They 
observed  that  the  electrical  resistivity  measurements  were  over- 
whelmingly influenced  by  the  specific  surface  area  of  the  clay  par- 
ticles when  distilled  water  was  used,  but  at  high  salt  concentrations 
the  resistivity  measurements  were  independent  of  this  influence.   As 
illustrated  in  Table  1,  the  specific  surface  area  is  an  indicator 
of  cation  exchange  capacity.   These  authors  concluded  that  further 
work  was  necessary  in  order  to  establish  the  limits  of  this  influ- 
ence for  different  clays  at  different  concentrations. 

Atkins  and  Smith   studied  sodium  montmorillonite,  calcium 
montmorillonite  and  kaolinite,  which  have  cation  exchange  capacities 
of  74.0,  60.1  and  16.1  meq/100  g  respectively.   They  found  that  the 
Formation  Factor  remained  constant  even  though  the  interstitial 
water  conductivity  ranges  were  1.1  to  8.7,  0.1  to  19.7,  and  0.009 
to  8.5  mho/meter  respectively.   They  found  the  same  result  with 
illite  clay,  but  the  range  of  interstitial  water  conductivity  used 

was  not  given. 
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Kermabon  et  al.    further  studied  the  effect  of  ionic  soiu- 
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tions  on  the  cation  exchange  capacity  of  clay  particles.   In  their 
laboratory  experiments  the  marine  clay  particles  had  a  grain  size 
smaller  than  3.9  microns  and  the  sodium  chloride  solutions  ranged 
from  0.7  to  10  percent.   The  porosity  of  the  clay  sediment  samples 
was  constant.   Throughout  this  wide  range  of  water  salinity  the  re- 
lation between  the  conductivity  of  sediment  and  the  conductivity  of 
the  interstitial  water  remained  constant;  therefore,  any  influence 
the  charge  of  the  clay  particles  had  on  conductivity  was  masked  by 

the  ionic  solutions. 

13 
Smith   found  similar  results  in  laboratory  experiments  with 

"coarse"  clay.   He  used  sodium  chloride  solutions  ranging  from 
0.5  to  10  percent. 

The  results  of  these  studies  indicate  that  in  fresh  water  the 
charge  on  clay  particles  (cation  exchange  capacity)  greatly  influ- 
ences the  conductive  properties  of  the  sediment.   This  charge  and 
consequently  its  effect  on  electrical  resistivity  measurements  in 
fresh  water  depends  on  the  amount  and  type  of  clay  present.   It  is 
apparent  that  in  a  fresh  water  environment  cation  exchange  capacity 
must  be  accounted  for  by  the  model  or  by  empirically  measuring  and 
correcting  for  its  effect  on  resistivity  measurements  if  accurate 
porosity  predictions  are  to  be  made.   However  in  a  marine  environ- 
ment, the  influence  of  the  cation  exchange  capacity  of  clays  on 
electrical  resistivity  measurements  appears  to  be  negligible  in  most 
cases.   The  reason  for  this  difference  between  fresh  and  salt  water 
can  be  deduced  by  visualizing  the  current  flow  in  solution  as  the 
migration  of  free  ions.   In  a  sediment  of  clay  and  fresh  water,  most 
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of  the  ions  in  solution  are  under  the  influence  of  clay  particles 
and  therefore  the  current  flow  through  the  sediment  is  inhibited. 
If  the  interstitial  xvater  is  saline,  there  are  sufficient  ions  to 
neutralize  the  clay  particles  and  allow  the  current  to  flow  freely. 
If  the  clay's  cation  exchange  capacity  is  high,  more  ions  will  be 
required  to  neutralize  its  effect  on  conductivity.   Although  cation 
exchange  capacity  is  usually  neglected  in  the  marine  environment, 
more  testing  is  needed  before  this  influence  on  electrical  resisti- 
vity measurements  for  all  types  of  clays  can  be  assumed  to  be 
negligible. 

The  above  studies  indicate  that  the  cation  exchange  capacity 
of  some  clays  need  not  be  considered  as  an  influence  on  electrical 
resistivity  measurements  in  the  salinity  of  ocean  water.   However 
other  characteristics  of  clay  must  be  considered  if  accurate  poros- 
ity predictions  are  to  be  made  from  electrical  resistivity  measure- 
ments.  The  clay  type,  sand  and  silt  present  in  a  sediment  determine 
the  size,  shape,  and  distribution  of  particles,  and  influence  parti- 
cle structure. 

The  M.I.T.  soil  classification  of  sediment  particle  size 
is  given  below: 

from  0.06  mm  to  2.0  mm  is  considered  sand, 

0.002  mm  to  0.06  mm  is  considered  silt,  and 

a  particle  smaller  than  0.002  is  considered  clay  sized. 
The  particle  distribution  of  a  sediment  is  defined  as  a  measure  of 
spread  such  as  standard  deviation  from  the  mean.   If  the  results  of 
past  work  are  applied  to  sediments,  it  would  appear  that  the  parti- 
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cle  size  and  distribution  might  influence  electrical  resistivity 
in  the  following  ways: 

If  the  porosities,  particle  shapes,  and  standard  deviations 
of  two  distributions  are  the  same  for  samples  with  different  mean 
particle  sizes,  the  electrical  resistivity  will  usually  be  lower 

for  the  larger  particle  size.   This  effect  was  suggested  In  experi- 
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ments  with  glass  beads  by  Meredith  and  by  Wyllie  and  Gregory. 

As  the  particle  size  increased,  the  electrical  resistivity  de- 
creased. 

If  the  porosity  and  mean  particle  shape  and  size  are  the  same 
while  the  distributions  of  two  sediments  are  different,  the  electri- 
cal resistivity  will  be  higher  in  the'sediment  with  the  greater  as- 
sortment of  particles.   This  effect  was  suggested  by  studies  on  rocks 
by  Semenov  "  and  work  conducted  on  glass  beads  by  Meredith. 

The  effect  of  shapes  in  a  sediment  is  complicated.   A  sedi- 
ment may  be  made  up  of  particles  of  different  shapes  ranging  from 
spheres  to  plates.   Usually  particles  in  the  sand  and  silt  size  range 
are  closer  to  a  sphere  in  shape  than  the  plate-like  clay  particles. 
The  different  particle  shapes  affect  the  path  lengths  that  migrating 
ions  take  when  current  flows.   If  two  sediments  have  approximately 
the  sane  particle  size  and  spread  and  equal  porosities  but  differ  in 
particle  shape,  the  electrical  resistivity  would  be  expected  to  be 
higher  for  a  sediment  composed  of  disc-shaped  particles  than  a  sedi- 
ment made  up  of  spherical  shaped  particles.   This  effect  was  sug- 
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gested  by  studies  on  various  shapes  by  Wyllie  and  Gregory,     Discs, 

cubes,  cylinders,  and  triangular  prisms  resulted  in  higher  electri- 


16 

cal  resistivity  values  than  spheres  at  the  same  porosity  values. 
They  also  found  that  Ottawa  sand  (sphericity  -  0.895)  and  glass 
beads  (sphericity  =  1.00)  had  similar  Formation  Factor  versus  poros- 
ity relationships. 

Meredith   found  that  the  conductivity  of  dispersions  at  a 
porosity  of  50  percent  resulted  in  a  10  to  20  percent  error  if  the 
shape  difference  between  spheres  and  discs  was  not  taken  into  con- 
sideration.  The  error  decreased  with  increasing  porosity. 

Semenov    theoretically  developed  equations  for  formation 
rocks  based  on  shape.   Semenov' s  equations  for  rocks  of  spherical 
shaped  grains  and  disc  shaped  grains  (having  the  same  distribution) 
resulted  in  Formation  Factor  values  of  2.7  and  10.1  respectively  at 
a  porosity  of  50  percent. 

Atkins  and  Smith  '  stated  that  at  the  same  porosity  angular 
particles  have  higher  Formation  Factors  than  spherical  shaped  par- 
ticles.  A  shape  factor  was  determined  for  five  clays  and  one  mica. 
Shape  was  inferred  to  be  the  major  factor  causing  differences  be- 
tween the  clays.   No  electron  microscope  photographs  were  taken,  and 
data  such  as  mean  particle  size  and  distribution  were  not  given. 

In  this  paper,  particle  structure  will  refer  to  the  orienta- 
tion of  individual  soil  particles  in  the  sediment  mass.   Particle 
structure  in  sands  and  silts  usually  depends  on  packing.   The  type 
of  packing  is  important  because  it  affects  the  porosity  and  the  path 
length  required  for  migrating  ions  and  consequently  the  electrical 
resistivity  measurement.   The  two  extreme  types  of  packing  of  uni- 
formly sized  spherical  particles  are  loose  or  cubic  packing  and 
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dense  or  rhombohedral  packing  as  shown  in  Figure  1. 


FIGURE  1 


SPHERICAL  PARTICLES  IN  A  DENSELY  AND  A 
LOOSELY  PACKED  CONFIGURATION 


CUBIC  PACKING  RHOMBOHEDRAL  PACKING 

(loose)  (dense) 


Structure  in  clays  is  greatly  influenced  by  those  factors  af- 
fecting the  size  of  the  double  layer.   In  clays  the  net  electrical 
force  between  particles  is  the  result  of  many  repulsive  and  attrac- 
tive forces.   If  the  double  layer  of  clay  particles  is  small,  the  at- 
tractive forces  are  dominant  and  the  clay  particles  will  move  toward 
each  other  and  become  attached  (Flocculant  Structure).   If  the  double 
layer  is  large,  the  repulsive  force  is  more  dominant  and  the  clay 
particles  will  tend  to  move  away  from  each  other  (Dispersed  Struc- 
ture).  The  two  clay  structures  are  shown  in  Figure  2. 

FIGURE  2 

ORIENTATION  OF  CLAY  PARTICLES  IN  A 
FLOCCULANT  AND  A  DISPERSED  STRUCTURE 


FLOCCULANT  DISPERSED 

STRUCTURE  STRUCTURE 
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In  the  flocculant  structure,  large  voids  may  be  enclosed  by 
the  clay  skeleton  and  result  in  sediments  of  very  high  porosities. 
On  the  other  hand  the  dispersed  structure  results  in  lower  porosity 
values.   Aside  from  the  facts  of  porosity  differences,  it  would  ap- 
pear by  visualizing  the  two  structures  that  the  flow  of  current 
through  these  two  structures  would  also  be  different. 

If  the  current  path  becomes  more  tortuous  its  length  is 
increased.   Tortuosity  of  the  current  path  may  be  defined  as  the 
ratio  of  the  average  statistical  length  of  current  flow  lines  to  the 
simple  physical  length  in  a  straight  line  from  one  end  of  a  sediment 
sample  to  the  other. 

By  definition  the  resistance  (R)  between  flat  faces  of  a 
cylindrical  sediment  sample  is: 


r  =  r  I  =  r  —■  ,  (10) 

m  A    c  A   '      > 
e 


where: 

1  is  the  length  of  the  sediment  sample, 
A   is  the  cross-sectional  area  of  the  sediment  sample, 
1   is  the  effective  average  statistical  length  of  the  sedi- 
ment through  the  pore  structure, 

A  is  the  effective  average  statistical  cross-sectional  area 
e 


of  the  pore  structure. 
In  terms  of  Formation  Factor: 
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1 
e 

FT  -  ^  -  — L-  -  -£-,  (11) 


r     A       A   ' 
c     c       e 


where  t  is  the  tortuosity  factor  previously  defined.   Porosity  can 
be  expressed  as  follows.: 


V   .  ,    A  1    A 
voids    e  e    e  /ii\ 

n=  — —  =  _  =  T.  t.  (12) 

total 


If  the  porosity  is  the  same  equations  (1])  and  (12)  can  be  combined 
and  tortuosity  may  be  defined  as  a  function  of  porosity  and  Formation 
Factor: 


or 


t  -,/FF'n  .,  (13) 

t2 

FF  =  —.  (14) 

n 

As  the  structure  of  a  sediment  becomes  more  complex  and  the 
current  paths  become  more  devious,  the  electrical  resistivity  will 
increase.   It  is  at  this  point  that  the  knowledge  of  both  the  elec- 
trical and  hydraulic  properties  of  sediments  and  the  relationship  be- 
tween them  becomes  important  if  physical  properties  are  to  be  proper- 
ly interpreted  by  electrical  resistivity  measurements.   Much  work  has 
been  done  on  the  flow  problem  through  porous  media.   The  general  flow 
formula,  be  the  flow  fluid,  electrical,  thermal  or  chemical,  is: 

F  =  L'X»A.  (15) 
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F  is  the  flow  rate, 

X  is  the  gradient, 

A  is  the  cross-sectional  area, 

L  is  a  coefficient  of  conduction. 
Numerous  equations  directly  inter-relating  hydraulic,  thermal,  and 
electrical  properties  of  porous,  solution-saturated  materials  have 
been  derived.   For  example,  Hutt  and  Berg   developed  prediction 
equations  by  investigating  thermal  and  electrical  conductivities  of 
sandstone  rock  and  ocean  sediment  and  relating  the  two  conductivi- 
ties through  the  common  parameter  of  porosity.   Although  such  a  sys- 
tem seems  reasonable  it  is  important  to  remember  that  the  success  of 
any  electrical  method  for  measuring  thermal  conductivity  in  a  sedi- 
ment would  depend  on  the  range  of  electrical  conductivity  of  the 
soil  particle  encountered  in  a  given  situation.   Therefore  it  is 
imperative  when  inter-relating  flow  problem  that  these  differences 
be  understood. 

Schopper   contends  that  if  quantitative  evaluations  between 
electrical  and  hydraulic  properties  are  to  be  improved,  the  differ- 
ences between  electrical  and  hydraulic  tortuosity  should  be  taken 
into  consideration.   The  essential  difference  is  that  hydraulic 
tortuosity  takes  into  account  a  direction  factor  in  addition  to  the 
normal  tortuous  path  length,  while  the  normal  tortuous  path  length 
is  the  only  factor  involved  in  electrical  tortuosity.   Hydraulic 
permeability  is  a  good  indicator  of  hydraulic  tortuosity,  and  there- 
fore improved  relationships  between  porosity  and  Formation  Factor 
would  take  into  account  hydraulic  permeability. 
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1  O  1  Q  jr\ 

Fatt   and  Wyble   and  Glanville    investigated  the  effects 
of  applied  pressure  on  electrical  resistivity  values  of  reservoir 
rocks  and  sandstone.   These  rocks  of  low  porosity  were  placed  under 
very  high  pressures  ranging  from  1,000  to  12,000  p.s.i.g.   The  ef- 
fect of  pressure  on  electrical  resistivity  seems  to  be  much  greater 
than  the  effect  on  porosity,  since  electrical  resistivity  was 
increased  much  more  than  porosity  decreased.   Similar  work  on  sedi- 
ments appears  to  be  nonexistent. 

Mathematical  Models 

When  selecting  mathematical  models  which  best  fit  a  porous 
medium  such  as  a  saturated  sediment,  it  is  necessary  to  be  aware  of 
the  influence  of  the  charge  on  clay  particles,  clay  type,  salinity 
of  the  continuous  phase,  sediment  particle  size,  shape,  distribution, 
and  structure.   In  order  to  further  refine  models,  awareness  of  other 
factors  such  as  hydraulic  permeability  and  consolidation  pressure 
may  be  useful.   In  addition  the  assumptions  of  theoretical  models 
need  special  consideration,  and  extrapolation  beyond  the  range  of 
experimentation  (in  empirical  models)  is  not  recommended. 

The  electrical  conductivity  of  dispersions  has  been  studied 
extensively.   Some  studies  date  back  over  a  hundred  years.   Only 
those  theories  that  are  most  likely  to  apply  to  sediments  are  con- 
sidered below, 

21 
In  1801  Maxwell  '  conducted  one  of  the  earliest  theoretical 

investigations  and  derived  an  equation  for  the  assemblage  of  spheres 

in  a  continuous  media.   The  assemblage  was  a  suspension  of  spheres 
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that  were  not  in  physical  contact  with  each  other.   Maxwell's 
equation,  in  non-dimensional  conductivity  terms,  follows: 


3  K .  -  2  nK.  4-  2  n 

\   "   A   -  n  +  3 •  <16> 

d 


If  the  spheres  are  assumed  to  be  non-conductive,  i.e.  K.  approaches 
zero,  Maxwell's  equation  reduces  to: 


K  =  tj2-2-  ,  (17) 

m   3  -  n  '  ' 


or  in  resistivity  terms: 


FF  =  R  =  V"^  •  <18> 

m    2  n 


This  equation  can  be  applied  to  a  sediment  mixture  consisting  of  a 
fluid  of  relatively  low  electrical  resistivity  and  a  network  of  soil 
grains,  not  in  contact  with  each  other,  that  are  insulators  of  high 
electrical  resistivity.   For  instance  soil  granules  composed  of  the 
following  common  minerals  shown  in  Table  2  would  have  a  high  resis- 
tivity. 

TABLE  2 

ELECTRICAL  RESISTIVITY  OF  SOME  COMMON  SOIL  MINERALS 

Mineral  Resistivity  (ohm  meter) 

11     12 
Clay  minerals      10   -  10 

Quartz  1012  -  1014 

Calcite  107  -  1012 
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Wh 


en  these  soil  granules  of  high  resistivity  (r  ,)  are  placed  in  a 
fluid  such  as  sea  water  which  has  a  relatively  low  resistivity  (r  ) 
(0.25  ohm  meters),  the  non-dimensional  resistivity  of  the  dispersed 
phase  (soil  granule)  can  be  expressed  as: 


Rd  =  ^ >«,  (19) 

c 


or  the  inverse  relationship: 

Kd  =  IT >  °'  (20) 

c 

22  23 

Wiener   and  Woodside  and  Messmer  "  developed  the  Series  and 

Parallel  Model  Equations  which  are  considered  the  limiting  relation- 
ships for  porosity  determination  by  means  of  conductivity.   A 
complete  development  of  these  equations  is  given  in  Appendix  A. 

In  the  Parallel  Equation,  the  interstitial  water  and  solid 
soil  particle  are  assumed  to  be  parallel  to  current  flow.   The  re- 
sulting equation  in  non-dimensional  form  follows: 


(Parallel)   K  =  K,  +  n(l  -  K,).  (21) 

m    d         d 


The  Series  Equation  evolves  when  the  solid  soil  particles 
and  interstitial  water  are  assumed  to  be  perpendicular  to  the  flow 
of  current: 

(Serles)  Km  .  1+n;L  .  iy  .  (22) 

d 
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22 
Wiener   recognized  that  other  models  would  fall  between  the  Series 

and  Parallel  Equations.   He  proposed  a  linear  combination  of  the  two 

models.   Later  this  linear  combination  was  referred  to  as  the 

23 
Three-Resistor  Model  by  Woodside  and  Messmer.    The  Three-Resistor 

Model  is  given  below: 


K  =  p  (Parallel  Model)  +  (1  -  p)  (Series  Model).   (23) 


In  this  combination  p  represents  the  fraction  of  the  sediment 

described  by  the  Parallel  Model, 
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Woodside  and  Messmer   also  proposed  a  Geometric  Equation 

in  which  the  sediment  is  assumed  to  be  the  geometric  mean  of  the 
conductivities  of  the  components  of  the  sediment: 


K  =  K.(L"n).  (24) 

m    d 


However  if  it  is  assumed  that  the  solid  soil  particle  is  non- 
conductive,  then  K,  approaches  zero.   Consequently  the  Series 
Equation  and  the  Geometric  Equation  have  no  meaning,  and  the 
Parallel  Equation  is  reduced  to: 


K  =  n,  (25) 

m    ' 


and  the  Three-Resistor  Model  is  reduced  to: 


K  =»  pn.  (26) 

m   r 


An  extensive  review  of  the  work  up  until  1959  including 
investigations  by  Rayleigh  and  Runge  was  conducted  by  Meredith. 
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He  concluded  that  effects  of  particle  size,  shape,  distribution  and 
structure  cannot  be  ignored  if  accurate  theoretical  predictions  of 
dispersions  are  to  be  made.   Meredith  developed  the  Distribution 
Model  for  treating  a  concentrated  dispersion  of  spherical  and  non- 
spherical  particles.   The  relationships  in  the  Distribution  Model 
that  are  most  applicable  to  a  sediment  environment  follow: 


nonconducting  .  . 

spheres:         K  =  ,"    .    .1    , r—  »  (27) 

r  m   (5  -  n)  (3  +  n) 


nonconducting 

randomly  oriented  n      ,          .  x 

j             v  9n  (n  +  1)                 /0QN 

rods:  K  =  -r> x  /c  , r  ,  (28) 

m  (4  -  n)  (5  +  n)  ' 


rods  perpendicularly 

oriented  to  field  ,   ,  .  x 

current:  K  =  — \ .  (29) 

m  3  -  n 


Although  the  above  relationships  are  presented  in  non-dimensional 
conductivity  terms,  it  is  simple  to  convert  them  to  Formation  Factor 
by  using  the  relationship: 


FF  =  £— .  (7) 

m 


Formation  Factor-porosity  relationships  proposed  by  Maxwell,  Wiener, 

Meredith,  Woodside  and  Messmer  are  presented  in  Figure  3.   Wherever 

K,  is  used  in  formulation  it  is  assumed  to  be  0.01.   Since  this  value 
d 

is  assumed  to  be  very  small,  0.01  was  arbitrarily  chosen. 

Electrical  resistivity  methods  became  established  as  a  Geo- 


26 


90    100 


POROSITY  (%) 


Picture  3.   Formation  Factor-Porosity  Relationships  Proposed  by 
Maxwell,  Wiener,  Meredith,  woodside  and  Mossmer. 
(In  the  Series  and  Geometric  equations  the  value  of 
K.  =  0.01  is  assumed;  in  all  other  equations  Kd  =  0. ) 
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physical  exploration  technique  in  1927  when  Conrad  Schumberger 
showed  that  electrical  resistivity  measurements  of  an  oil  well 
(then  called  electrical  coring),  could  distinguish  between  produc- 
tive formation  rocks  and  those  that  are  non-productive.   This 
development  stimulated  extensive  research  on  electrical  resistivity 

measurements  as  an  indicator  of  physical  and  chemical  properties. 

24 
In  1941  Archie   defined  formation  resistivity  factor  i.e., 

Formation  Factor,  as: 


r 

FF  =  —  ,  (30) 

r 
c 


where  r   is  the  resistivity  of  the  saturated  sediment  and  r   is  the 
m  '  c 

resistivity  of  interstitial  water.   Archie  also  empirically  deter- 
mined an  equation  relating  FF  to  the  formation  porosity: 

FF  =  n"m,  (31) 

where  n  is  the  porosity  of  the  formation,  and  m  is  the  slope  of  the 
line  in  a  log-log  plot  of  Formation  Factor  versus  porosity  diagram. 
From  his  experimental  work,  Archie  determined  m  to  be  in  the  range 
of  1.8  to  2.0  for  consolidated  sandstones,  and  about  1.3  for  loosely 
consolidated  sands, 

Winsauer  et  al.   investigated  FF  to  n  plots  and  suggested  the 
relation  should  be: 


FF  =  an""1  ,  (32) 
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where  (a)  is  not  always  unity,  but  is  often  a  factor  accounting  for 
the  amount  of  cementation  between  mineral  grains.   Therefore  (a)  is 
referred  to  as  the  cementation  factor.   From  equation  (32)  the 
"Humble  Formula"  was  developed  and  is  probably  the  most  widely  used 
equation  in  the  well  logging  industry: 


FF  =  0.62n"2,15  .  (33) 


A  comprehensive  review  of  geophysical  well  logging  through  1961  has 
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been  compiled  by  Dakhnov.    He  summarized  the  factors  that  affect 

electrical  resistivity  of  a  porous  media  by  a  general  formula.   If 

the  equation  is  applied  to  a  saturated  sediment,  it  becomes: 


rm  =  fl(c)  f2(n)  f3(S)  VT)  f5(Q)  f6(rd}  f7(rc)f  °4) 


where  r   is  the  resistivity  of  the  sediment, 

m  '  ' 

c  the  amount  of  clay/silt  in  the  sediment, 

n  the  porosity  of  the  sediment, 

S  the  partial  saturation  of  the  sediment, 

T  the  temperature  of  the  sediment, 

Q  the  cation  exchange  capacity  of  the  soil  mineral, 

r,  the  resistivity  of  the  soil  mineral, 

r   the  resistivity  of  the  interstitial  water. 
c 

Since  some  of  the  variables  are  sometimes  omitted,  the  complexity  of 
the  above  equation  is  often  reduced.   For  example  if  the  soil  miner- 
al (r.)  is  nonconductive,  this  term  is  omitted.   In  a  saturated 
sediment  (S)  is  a  constant  equal  to  one.   If  the  environment  is 
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saline  the  interstitial  water  reduces  the  effect  of  cation  exchange 
capacity  (Q)  and  cation  exchange  capacity  is  usually  ignored.   How- 
ever if  the  interstitial  water  has  a  low  ionic  concentration  this 
term  must  be  considered.   The  temperature  (T)  can  usually  be  measured 
and  a  correction  to  electrical  resistivity  can  be  made  based  on  the 
following  formula: 


r--  =  rr 


[1  +  t   (T  -  25°)],  (35) 


where  r_,.  is  the  resistivity  at  25  Centigrade, 

r    the  resistivity  measured  at  temperature  T, 
t   the  temperature  coefficient  of  resistivity,  which  usually 
is  about  2.5  percent  per  degree  Centigrade, 
In  most  situations  then,  the  resistivity  of  a  saturated  sediment 
(at  a  set  temperature)  will  be  determined  by  three  parameters  if  the 
environment  is  saline:   r  ,  c,  and  n;  and  by  four  parameters  if  the 
ionic  solution  of  the  interstitial  water  is  low:   r  ,  c,  Q  and  n. 

Dakhnov  also  obtained  a  theoretical  relation  between  Formation 
Factor  and  porosity  of  uniform  nonconducting  spheres  in  a  conducting 
fluid: 

FF  .  1+  0.25(1  -n)^\  (36) 

1  -  (1  -  n)Z/J 

Semenov   suggested  a  theoretical  equation  for  cemented  rocks 
consisting  of  grains  which  have  the  forms  of  ellipsoids  of  revolu- 
tion.  The  formula  for  the  Formation  Factor  along  the  axis  of 


30 


revolution  is: 


In  n 


_  f  +  v  .,  In  (1  -  v  ) 


In  this  equation  f  is  defined  as: 

3 


f  ■  -j 1,  (38) 

(1  -  e  )(arthe  -  e) 


where 


•/ 


h  2 
1  -  (?)  ,  (39) 


a 


e  is  the  eccentricity  of  an  ellipse,  a  and  b  being  semimajor  and 
semiminor  axis  of  the  ellipsoid.   The  extreme  cases  of  this  equa- 
tion are  when  f  =  0  (nonconducting  grains  of  a  flat  disc  shape) 
and  when  f  =  2  (nonconducting  grain  would  be  spherical).   This 
theoretical  work,  although  developed  for  the  electrical  resistivity 
of  rocks,  combined  two  additional  important  factors,  that  of  parti- 
cle distribution  which  is  termed  the  sorting  coefficient  (v  ),  and 
a  particle  shape  factor  (f).   These  two  factors  would  also  be 
important  in  a  sediment  environment. 

Formation  Factor,  porosity  relationships  proposed  by  Archie, 
Dakhnov  and  Semenov  are  present  in  Figure  4.   The  Humble  Formula 
is  also  presented  in  Figure  4. 
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Figure  4.   Formation  Factor-Porosity  Relationships  Proponed  by 
Archie,  Dakhnov,  Semenov  and  the  Humble  Formula.' 
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More  recent  works  which  are  not  reviewed  previously  in  the 
literature  are  discussed  below.   These  recent  works  are  either  di- 
rectly applicable  or  actually  conducted  on  sediments. 

Atkins  and  Smith   have  experimentally  found  m,  which  they 
call  shape  factor  exponents,  for  various  mineral  particles: 

Mineral  Particle        m 

Na  Montmorillonite  3.28 

Ca  Montmorillonite  2.70 

lllite  2.11 

Kaolinite  1.87 

Sand  1.60 

They  also  concluded  that  (a)  the  cementation  factor  in  formula  32 
was  not  unity  but  was  a  factor  of  proportionality  which  accounts  for 
amount  of  clay  to  sand  ratio.   Only  if  a  system  consisted  of  one 
type  of  mineral  particle  would  (a)  be  unity. 

It  is  interesting  to  note  that  if  m  could  be  found  by  experi- 
mentation, then  for  that  type  of  clay  only  the  determination  of  the 
Formation  Factor  would  be  necessary  and  the  porosity  could  be  esti- 
mated.  Although  this  study  was  conducted  in  order  to  understand 
how  the  clay  mineral  fraction  influences  reservoir  rocks,  the  results 
could  be  very  useful  in  the  prediction  of  porosity  from  Formation 

Factors  in  sediments. 

9 

Boyce  took  43  electrical  resistivity  measurements  on  sedi- 
ments taken  from  cores  obtained  from  seven  locations  in  the  Bering 
Sea,  some  of  which  were  600  nautical  miles  from  each  other.   Using 
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a  consolidometer  he  squeezed  the  interstitial  water  out  of  these 
sediments  and  measured  the  electrical  resistivity  of  this  water  in 
order  to  obtain  Formation  Factors.   The  porosity  of  these  measure- 
ments ranged  from  58.3  to  87.4  percent.   By  the  method  of  least 
squares  for  an  exponential  function  the  following  formula  was 
derived: 

FF  =  1.30n"1,45.  (40) 

Boyce  found  that  this  formula  had  an  error  of  +  15  percent.   He 
stated  that  this  rather  large  error  could  be  caused  by  different 
types  and  percentages  of  electrical  conductors,  clay  minerals,  and 
pore  textures.   The  large  error  suggests  that  600  miles  is  too  large 

an  area  to  attempt  to  define  with  one  formula. 

12 

Kermabon  et  al.   studied  the  electrical  resistivity  of  deep 

sediments  from  21  long  cores  taken  in  a  4  by  10  square  mile  area  in 
the  Tyrrhenian  Sea,   The  porosity  of  these  sediments  ranged  from 
50  to  87  percent.   By  assuming  the  resistivity  of  the  interstitial 
water  was  the  same  as  the  bottom  sea  water,  Formation  Factors  were 
determined.   More  than  2500  measurements  of  porosity  and  Formation 
Factors  were  compared.   A  third  degree  polynomial  curve  was  found 
to  best  fit  their  data  relating  porosity  and  Formation  Factor 


n(7.)  =  -5.9021  FF3  +  40.0416  FF2  -  105.3889  FF  +  171.2504  ,  (41) 


or, 


m   1  +0.7193  (1-n1-4615) 

rr  1.4615        '  ^    ' 

n 
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According  to  Kermabon  et  al.  the  above  equation  fit  the  data 

obtained  with  only  a  +  2.5  percent  porosity  prediction  error. 

13 

Smith   studied  electrical  resistivity  measurements  of  North 

Atlantic  deep-sea  cores  and  in  situ  resistivity  measurements  on  the 
continental  shelf  around  Wales,  England.   The  resistivity  of  the 
interstitial  water  was  either  estimated  from  bottom  water  salinity 
or  actual  core  water  salinity  measurements  were  taken  in  the  labora- 
tory.  A  best  fit  straight  line  through  99  measurements  resulted  in 
the  following  equation: 


FF  =  1.35n"1,2.  (43) 


Smith  noted  that  the  data  clearly  grouped  into  two  classes:   for 
silts  and  clays  of  porosities  greater  than  60  percent, 


FF  =  n"2,  (44) 


and  for  sands  and  coarse  silts  of  porosities  less  than  60  percent 


FF  =  n"1,5.  (45) 


He  concluded  that  there  was  no  need  to  use  a  polynomial  formula  of 

12 
the  type  described  by  Kermabon  et  al.,   but  that  once  a  Formation 

Factor  is  measured  an  approximate  value  of  porosity  could  be  ob- 
tained from  a  general  empirical  equation,  and  that  after  selective 
sampling  a  better  value  for  porosity  could  be  obtained  by  a  more 
refined  equation. 

Porosity  versus  Formation  Factor  equations  are  plotted  in 
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Figure  5  from  Boyce,  Kcrniabon  et  al.f  and  Smith.   Figure  6  shows 
Atkins  and  Smith's  formulations  using  equation  (32)  for  sodium 
montmorillonite  (m  =  3.28),  illite  (m  =  2.11),  kaolinite  (1.87) 

and  sand  (m  =  1.60). 

26 
Pautot   was  able  to  correlate  several  sedimentary  layers  by 

taking  electrical  resistivity  measurements  on  sediment  cores  ob- 
tained some  distance  apart.   He  did  this  work  both  in  a  lake  and  in 
a  marine  environment.   Although  no  formulation  between  electrical 
resistivity  and  sediment  properties  was  given,  the  ability  to  match 
horizontal  layers  may  indicate  a  similarity  between  physical  and 

chemical  properties  with  electrical  properties  of  sediments. 

27  28 

Bouma  et  al.   and  Chmelik  et  al.   developed  laboratory  and 

in  situ  electrical  resistivity  equipment  for  marine  sediments. 
Their  objective  was  to  make  correlations  on  lithology  and  geotech- 
nical  properties.   Although  no  formulations  were  given,  they  at- 
tempted to  correlate  electrical  properties  with  the  chemical,  phys- 
ical and  engineering  characteristics  of  sediments.   Properties  such 
as  pH,  Eh,  water  content,  density,  carbonate  content,  grain  size 
analysis,  X-ray  radiographs,  photographs,  cone  penetrameter  and 
vane  shear  measurements  were  compared  with  electrical  resistivity 
measurements.   Certain  trends  were  apparent  from  their  data.   Elec- 
trical resistivity  was  indirectly  proportional  to  the  water  content 

and  percentage  of  clay  size  present. 

29 
Ball   studied  the  literature  in  order  to  determine  the 

feasibility  of  using  electrical  conductivity  techniques  for  sediment 

porosity  prediction.   He  concluded  that  the  following  points  need 
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Figure  5.   Formation  Factor-Porosity  Relationships  Proposed  by 
Boyco,  Kermabon,  et.  al. ,  and  Smith. 
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Figure  6.   Formation  Factor-Porosity  Relationships  Proposed  by 
Atkins  and  Smith  for  Sodium  Montmorilloniter  Tlliter 
Kaolinito  and  Sand. 
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further  work  both  in  the  laboratory  and  in  the  marine  environment. 

1.  What  theoretical  or  empirical  equation  best  describes 
the  electrical  conductivity-porosity  relationship? 

2.  What  is  the  effect  of  clay  on  electrical  conductivity 
measurements? 

3.  How  can  the  conductivity  of  interstitial  water  in  marine 
sediments  be  conveniently  determined? 

4.  What  is  the  best  method  for  restricting  the  volume  over 
which  a  measurement  would  be  made? 

The  first  two  questions  are  considered  in  this  section  of  the  paper, 
while  the  last  two  questions  are  investigated  in  the  second  section. 
The  main  motivation  behind  this  literature  review  and  the  following 
laboratory  experiments  is  to  justify  the  use  of  electrical  resisti- 
vity techniques  as  an  in  situ  measuring  method  to  predict  porosity 
in  a  marine  environment.   Consequently  sediment  type,  interstitial 
water  salinity,  and  formulations  applicable  to  the  marine  sediment 
environment  were  studied  more  extensively  than  those  applicable  to 
a  fresh  water  environment. 


EXPERIMENTAL   PROCEDURE 

Materials 

The  following  types  of  materials  were  obtained  for  testing: 

1.  Kaolinite  clay  (Albion  Kaolin  from  Babcock  and  Wilcox 
Refractories  Division), 

2.  Illite  clay  (Grundite  from  A.  P.  Green  Refractories 
Company), 

3.  Montmorillonite  clay  (Black  Hills  Bentonite  from 
International  Minerals  and  Chemical  Corporation), 

4.  Providence  silt, 

5.  Ottawa  sand  (well-rounded  grain  shape), 

6.  Two  types  of  glacial  sands  (angular  grain  shapes), 

7.  Four  types  of  marine  sediments. 

One  type  of  marine  sediment  was  obtained  by  a  Ewing  gravity 
corer  and  the  Deep  Ocean  Sediment  Probe  (DOSP)  in  April  1970  from 

shallow  water  near  Ponce,  Puerto  Rico.   The  properties  of  this  sedi- 

30 
ment  were  extensively  investigated  by  Lewis.    The  other  three 

types  of  marine  sediments  were  taken  by  an  Electrical  Resistivity 
Measuring  System,   These  samples  were  taken  during  the  Fall  of  1971 
from  Narragansett  Bay. 

The  properties  of  the  marine  sediments  and  the  other  labora- 
tory materials  are  presented  in  Table  3.   Co-ordinates  of  marine 
sediment  cores  measured  in  this  investigation  are  given  in  Table  4. 
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TABLE  3 
MINERALOGY,  CLASSIFICATION,  GRAIN  SIZE,  SHAPE,  SPECIFIC  GRAVITY, 


NAME  OF  SEDIMENT 

MINERALOGY1 

SEDIMENT  CLASSIFICATION 

LIQUID 
LIMIT 

7.  Sand 

7.  Silt 

7.  Clay 

Kaolinite 

Kaolinite 
Quartz 

2 

68 

30 

37.4 

Illite 

Illite 
Quartz 

7 

53 

40 

44.8 

Montmorillonite 

Montmorillonite 

K-Feldspar 

Quartz 

10 

12 

78 

135.0 

Providence  Silt 

Quartz 
Plagioclase 

4 

92 

4 

28.5 

Ottawa  Sand  (a) 

Quartz 

100 

- 

- 

- 

Ottawa  Sand  (b) 

Quartz 

100 

- 

- 

- 

Glacial  Sand(l-a) 

Quartz 

100 

- 

- 

- 

Glacial  Sand(l-b) 

Quartz 

70 

30 

- 

- 

Glacial  Sand  (2) 

Quartz 

98 

2 

- 

- 

Puerto  Rican 
Sediment 

Aragonite 
Calcite 
Quartz 
Plagioclase 

70  -  50 

18  -  24 

7-18 

38.0 

to 
72.2 

Narragansett  Day 
(Station  A) 

Quartz 

K-Feldspar 

Mica 

41  -  33 

57  -  51 

7  -  14 

31.0 

to 
36.0 

Narragansett  Bay 
(Station  B) 

Quartz 

K-Feldspar 

Mica 

58  -  43 

34  -  53 

4  -  9 

27.0 

to 
28.1 

Narragansett  Bay 
(Station  C) 

Quartz 

97 
_ 

3 

- 

- 

1.   From  X-ray  Diffraction  analysis  (see  Appendix  C). 
2*  fi  -  -    log-  (diameter  in  millimeters), 

jtf__  =  diameter  for  which  507.  of  the  sample  is  finer  grained. 
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LIQUID  LIMIT,  PLASTIC  LIMIT,  pH,  AMD  ORGANIC  CONTENT  OF  SEDIMENTS 


GRAIN  SIZE 

SHAPE4 

SPECIFIC 
GRAVITY 

PLASTIC 
LIMIT 

pH 

ORGANIC 

CONTENT 

(7.) 

<502 

•/ 

8.38 

1.81 

platy 

2.60 

26.3 

4.1 

- 

8.25 

2.99 

platy 

2.79 

26.5 

2.0 

- 

- 

- 

platy 

2.73 

54.1 

7.2 

- 

5.97 

1.21 

angu- 
lar 

2.75 

23.0 

5.6 

- 

.52 

.17 

well- 
rounded 

2.65 

- 

- 

- 

1.26 

.30 

well- 
rounded 

2.65 

- 

- 

- 

2.32 

.34 

angu- 
lar 

2.65 

- 

- 

- 

3.47 

.50 

angu- 
lar 

2.65 

- 

- 

- 

1.00 

1.31 

sub- 
angular 

2.65 

- 

- 

- 

4.4 

to 

6.2 

1.80 

to 
4.40 

- 

2.78 

to 
2.86 

to 
53.8 

7.3 

to 

7.8 

1 

A. 2 

to 

4.9 

2.20 

to 
2.70 

mixed 
platy 
&  angu- 
lar 

2.70 

27.0 

to 
33.9 

4.5 

to 

6.4 

5-10 

3.3 

to 

4.1 

1.90 

to 
2.20 

mixed 
angular 
&  platy 

2.70 

to 
26.1 

6.2 

to 

6.6 

4-5 

2.18 

1.10 

angular 

2.65 

- 

7.3 

- 

3.  6  .  (Dispersion  measure)  =  1/2  (&q, 

4.  Obtained  by  electron  microscope 
and  conventional  microscope. 


<*16> 
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TABLE  4 
LATITUDE   AND    LONGITUDE   OF   MARINE   SEDIMENTS 

I.      Puerto  Rico  Cruise,   April  1970 

1.  17°57.4'    N      -     66°37.4'   W 

2.  17057.2f   N      -     te°21.T   W 

3.  17°56.5'    N      -      66°36.6'    W 

4.  17°56.2'    N      -     66°35.0'   W 

5.  17054.4f    N      -     66°35.5'   W 

II.  Narragansett  Bay,  Fall  1971 

Station  A  -  41°34.0'  N  -  7l°25.4*  W 

Station  B  -  41°29.5'  N  -  71°24.7'  W 

Station  C  -  41°22.7'  N  -  71°30.6I  W 


43 


The  details  of  the  tests  used  to  obtain  these  properties  are  given 
in  Appendix  C.   These  include  X-ray  analysis,  grain  size  analysis, 
hydrometer,  specific  gravity,  liquid  and  plastic  limit  determina- 
tion, pH,  and  organic  content. 

The  illite  clay,  kaolinite  clay  and  Providence  silt  were 
made  as  nearly  homoionic  to  sodium  ion  as  possible  by  washing 
batches  of  clay  in  concentrated  (IN)  solutions  of  sodium  chloride. 
The  excess  salt  was  then  removed  by  leaching  clays  and  silt  with 
distilled  water  and  the  conductivity  of  the  solution  was  checked 
after  each  dilution  by  means  of  a  Beckman  Conductivity  Bridge 
(Model  RC-16B2).   When  the  solution  conductivity  was  that  of  normal 
tap  water,  the  clays  and  silt  were  subsequently  dried  in  an  oven  at 
105  C,  lightly  pulverized,  and  stored  in  sealed  jars.   This  prepa- 
ration of  samples  was  not  conducted  on  Black^Hills  Bentonite  since 
it  is  a  sodium  montmorillonite  clay  and  therefore  is  already  homo- 
ionic  to  the  sodium  ion.   The  sands  were  sieved  to  insure  proper 
grain  size  and  then  washed  in  tap  water  to  remove  as  much  ionic 
contamination  as  possible.   The  sands  were  then  dried  and  stored. 

Laboratory  Apparatus 

The  laboratory  testing  apparatus  was  developed  in  the  follow- 
ing stages: 

1.  a  simple  apparatus  to  determine  the  feasibility  of  the 
method. 

2.  a  modified  design  to  obtain  variation  of  key  parameters. 

3.  redesigned  equipment  to  facilitate  testing  and  allow  an 
increased  number  of  tests. 
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Simple  apparatus 

The  first  resistivity  cell  was  constructed  from  a  seven  inch 
length  of  acrylic  tubing  (I.D.  1,38  in.,  O.D.  1.75  in.)   (see 
Figure  7).   The  aluminum  top  and  bottom  electrodes  were  machined  down 
to  fit  the  inner  diameter  of  the  tubing  with  0-ring  seals  to  allow 
for  a  tight  seal.   A  small  hole  was  drilled  in  the  top  electrode  to 
provide  an  escape  for  excess  water  while  the  cell  was  being  closed. 
The  four  side  platinum  electrodes  were  installed  (see  Figure  8)  by 
drilling  circumaxially  two  small  holes,  V8  inch  apart,  and  looping 
0.025  in.  platinum  wire  through  the  holes,  thus  providing  four  ap- 
proximate point  electrodes  spaced  3/4  inch  apart  vertically  along 
the  interior  of  the  cell  wall.   The  holes  were  then  made  air  tight 
by  using  a  GACO  sealant.   The  platinum  wire  was  then  soldered  to  a 
standard  copper  wire  in  order  to  connect  the  cell  to  the  instruments. 
For  testing  actual  marine  sediments  the  top  and  bottom  electrodes 
were  redesigned  in  order  to  allow  sediment  samples  of  various  lengths 
to  be  used  in  the  cell,   (See  Figure  9), 

Another  simple  testing  apparatus  (see  Figure  10)  was  made  by 
placing  top  and  bottom  current  electrodes  on  a  marine  sediment  sample 
located  in  the  original  core  liner  used  to  obtain  and  store  the 
sample.   Top  and  bottom  electrodes  were  machined  to  the  same  diameter 
as  the  core  liner,   A  hole  (3/16  in,  in  diameter)  was  drilled  along 
the  perimeter  of  the  top  electrode.  This  hole  allowed  an  insulated 
V8  inch  diameter  stainless  steel  tube  to  pass  vertically  through  the 
top  electrode  into  the  sediment.   Located  at  the  end  of  the  stain- 
less steel  tubing  were  two  vertically  separated,  platinum  potential 
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Figure  7.  The  Simple  Cell  Apparatus. 
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Figure  8.  The  Four  Side  Platinum  Electrodes  of  the 
Simple  Cell  Apparatus. 
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Figure  9.   simple  Cell  Apparatus  Used  to  Test  Samples 
of  Various  Lengths. 
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Figure  10.   Testing  Apparatus  and  Electronic  Equipment  Used 
to  Test  Marine  Sediments  in  the  Laboratory. 
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electrodes.   The  wiring  from  the  platinum  potential  electrodes 
passed  through  the  tubing  to  the  instrumentation. 

Modified  apparatus 

A  modified  design  was  required  in  order  to  test  various 
porosities  of  the  fine  grained  sediments  such  as  silts  and  clays. 
This  design  also  made  it  possible  to  see  if  consolidation  pressure 
altered  electrical  resistivity  measurements  to  such  an  extent  as  to 
mislead  the  porosity  prediction. 

Sedimentation  tubes  were  fabricated  with  a  piston  loading 
rod  which  applied  pressure  to  the  sediment.   This  pressure  could  be 
applied  incremently  by  adding  weights  to  the  loading  rod.   Each 
tube  was  constructed  from  a  nine  inch  length  of  acrylic  tubing 
(I.D.  ■  1.42  in.,  O.D,  =  1.75  in.).   The  tubes  were  bored  to  insure  a 
concentric  inner  diameter,  and  allow  the  piston  head  and  rod  assembly 
to  slide  freely  inside  the  tube.   Small  holes  were  drilled  into  the 
piston  head  to  allow  fluid  drainage.   The  piston  head  was  made  of  the 
same  material  as  the  tubes  to  prevent  any  binding  due  to  temperature 
changes.   A  bottom  cap  of  the  same  design  and  dimensions  as  the 
piston  head  was  placed. in  the  bottom  of  the  tubes  to  prevent  sediment 
loss  from  the  bottom  when  the  pressure  was  applied.   In  addition 
micro-pore  filters  were  used  above  the  bottom  cap  and  below  the  piston 
head  to  prevent  the  sediment  from  oozing  out  of  the  fluid  drainage 
holes.   A  PVC  (poly-vinyl-chloride)  plastic  piston  rod  was  screwed 
into  the  piston  head  and  supported  a  PVC  loading  plate.   Five  side 
platinum  electrodes  were  installed  in  the  same  manner  described  for 
the  first  cell.   Ribbon  wire  was  used  to  connect  the  four  cells  to 
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the  instrumentation  (see  Figures  11  and  12).   By  connecting  individual 
cells  and  by  using  a  simple  switching  device,  different  electrode 
combinations  could  be  used.   The  four  cells  conveniently  fit  into  a 
beaker . chamber  and  were  held  upright  by  a  plastic  support  frame. 
Four  of  these  chambers  were  fabricated  at  the  URI  Engineering  Instru- 
ment Shop.   Details  of  the  design  are  given  in  Appendix  B, 

The  electronic  equipment  (see  Figure  13)  used  with  the  simple 
apparatus  and  modified  apparatus  was  obtained  from  New  London  Labora- 
tory Naval  Underwater  Systems  Center,  New  London,  Connecticut.   The 
1  ma,  1000  Hz  signal  sent  to  the  current  electrodes  was  supplied  by 
a  constant  current  generator  and  power  supply.   The  voltage  differ- 
ence obtained  from  the  potential  electrodes  was  sent  to  an  AC/DC 
converter.   The  potential  difference  between  electrodes  was  measured 
from  a  digital  frequency  counter  and  the  same  signal  monitored  on  an 
oscilloscope.   The  electronic  instrumentation  was  calibrated  and  the 
accuracy  of  the  system  was  determined  to  be  within  +  1  percent. 

Redesigned  apparatus 

The  laboratory  apparatus  was  redesigned  in  order  to  reduce  the 

time  required  to  consolidate  the  fine  grained  sediments.   A  modified 

31 
version  of  a  pressure  vessel  designed  and  fabricated  by  Maus   was 

incorporated.   The  pressure  vessel  was  constructed  from  a  3}j  inch 

length  of  aluminum  pipe  which  was  capped  on  both  ends  with  cast  iron 

pipe  caps.   The  air  pressure  is  regulated  and  enters  the  vessel  through 

a  top  fitting.   The  porous  disc  sits  on  the  bottom  end  cap.   This  end 

cap  also  acts  as  a  reservoir  for  interstitial  water  draining  out  of 

the  sediment  when  the  pressure  is  applied.   Pressures  up  to  19  psig 
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Figure  11.  One  Cell  of  the  Modified  Test  Apparatus. 
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Figure  12.  One  Chamber  of  the  Modified  Test  Apparatus. 
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Figure  13.   Electronic  Equipment  Used  with  the  Simple 
and  Modified  Testing  Apparatus. 
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can  be  applied  without  bubbling  air  through  the  porous  disc.   The 
cell  was  constructed  of  thin  walled  plastic  tubing  (I.D.  =  1.38  in., 
O.D.  =  1.50in.)  and  5.5  inches  long.   Four  platinum  wire 
(dia.  =  0.025  in.)  electrode  rings  were  cemented  in  grooved  indenta- 
tions on  the  inside  tube  wall  with  epoxy.   The  vertical  separation 
between  rings  was  0.375  in.   The  cell  was  cemented  to  the  porous 
disc  with  epoxy.   The  platinum  ring  electrodes  and  thermistor  (used 
to  measure  temperature  changes)  were  connected  to  wiring  that  fed 
through  an  air  tight  seal  located  at  the  top  of  the  pressure  vessel. 
This  wiring  connected  the  cell  to  the  instrumentation.   Figure  14 
illustrates  the  cell,  pressure  vessel,  regulator  and  gauge  assembly. 
Appendix  B  contains  other  details-  of  design. 

New  electronic  instrumentation  allowed  continuous  electrical 
conductivity  measurements  during  the  consolidation  process  and 
recorded  the  results.   This  electronic  instrumentation  consisted  of: 
(1)  conductivity  signal  conditioner,  (2)  impedance  matching  amplifier, 
and  (3)  slow  speed  single  channel  recorder.   The  conductivity  signal 
conditioner  was  manufactured  by  Honeywell  Inc.  (Model  description 
552022-1002-100-DOO)  as  part  of  a  water  quality  measuring  system 
(see  Figure  10).   Since  the  signal  conditioner  is  basically  a  solid 
state  null  balancing  A.C.  amplifier  only  the  sensor  (electrodes)  had 
to  be  modified  to  allow  this  system  to  be  used  on  sediments.   The 
signal  conditioner  supplies  a  60  Hz  0.6  ma  current  to  drive  the 
sensor  array,  and  compares  the  resultant  output  of  the  sensor  voltage 
electrodes  with  a  stable  reference  voltage.   A  differential  voltage  is 
used  to  regulate  the  drive  current  and  it  produces  an  output  linearly 
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Figure  14.   Cell,  Pressure  Vessels,  Regulator  and  Guage 
Assembly  Used  .with  the  Redesigned  Equipment. 
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proportional  to  the  conductivity  of  the  measured  medium.   This  out- 
put voltage  can  be  read  on  the  signal  conditioner  panel  meter,  a 
digital  voltmeter,  or  sent  to  an  impedance  matching  amplifier  and 
then  to  a  chart  recorder  where  a  continuous  record  can  be  maintained. 

Since  sediment  conductivity  varies  with  operating  temperature, 
compensation  is  required  to  refer  the  measurement  to  standard  condi- 
tions (usually  25  C).   The  thermistor,  located  inside  the  pressure 
vessel  and  attached  to  the  plastic  cell,  directly  alters  sensor  cur- 
rent, thereby  compensating  for  temperature  changes  during  testing. 

According  to  Honeywell  operating  specifications,  the  calibrated 
accuracy  of  the  signal  conditioner  and  sensor  is  +  0.50  percent  for 
the  full  scale  range.   Calibration  curves,  schematic  diagrams  and 
detailed  operating  principles  are  presented  in  Appendix  B. 

Methods  of  Testing 

Simple  apparatus 

In  the  simple  apparatus  first  described,  an  electrical  cur- 
rent was  allowed  to  flow  between  the  top  and  bottom  electrodes  of  the 
resistivity  test  cell,  and  a  flow  pattern  developed  that  was  analo- 
gous to  Darcy's  Law  for  one  dimensional  fluid  flow.   In  this  study 
when  a  voltage  potential  was  measured  between  two  platinum  potential 
electrodes,  the  resistivity  of  the  sample  was  calculated  as  follows: 


r  = 

m 


^f  =  K  f  ,  (46) 


where: 
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r  =  the  resistivity  of  the  sample, 

A  =  cross  sectionai  area  of  the  ceil, 

a  ■  the  voltage  electrode  spacing, 

^E  =  voltage  difference  between  voltage  electrodes, 

I  =  the  current, 

K  =  the  geometric  cell  constant. 
If  the  current  was  allowed  to  flow  between  the  top  and  bottom 
platinum  electrodes,  a  different  flow  pattern  developed.   The  sec- 
ond part  of  Equation  46  is  still  applicable  if  the  voltage  potential 
is  measured  between  the  two  inner  platinum  electrodes.   The  geometric 
cell  constant  (K)  can  no  longer  be  determined  simply  by  dividing  the 
cross  sectional  area  by  the  distance  between  potential  electrodes 
but  is  usually  determined  empirically  by  using  solutions  of  known 
resistivities. 

The  cell  constant  was  determined  by  making  measurements  with 
the  cell  filled  with  standard  solution  of  sodium  or  potassium  chloride, 
for  which  the  resistivity  may  be  found  in  chemical  tables.   In  this 
study,  cell  constants  for  all  the  geometrical  electrode  configurations 
were  obtained  using  0.01N  KC1,  0.1N  NaCl,  0.5N  NaCl,  and  IN  NaCl  solu- 
tions.  All  cell  constants  of  simple  geometric  configuration  were 
compared  with  measured  values  and  less  than  1  percent  deviations  were 
encountered  between  the  two  values. 

Since  variation  of  interstitial  water  was  desired  in  the  labora- 
tory both  to  simulate  ocean  as  well  as  estuary  conditions  and  to 
observe  the  effect  solution  concentration  has  on  the  electrical  resis- 
tivity of  the  sediment,  sodium  chloride  solutions  of  3.5,  7.0,  17.5 
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and  35.0  parts  per  thousand  were  used  in  the  early  stages  of  testing. 
In  later  testing  standard  sodium  chloride  solutions  of  0.1,  0.2,  0.3, 
0.5  and  0.6  normality  were  used  as  the  interstitial  fluid  of  the 
sediment  matrix. 

In  early  laboratory  testing  using  the  simple  cell  the  electri- 
cal resistivity  of  sands  was  tested  using  the  following  procedure: 

The  interstitial  salt  solution  was  first  poured  into  the  cell 
and  the  sand  was  added  slowly.   The  sand  settled  out  rapidly  in  a 
loose  packed  configuration.   The  temperature  and  length  of  the  sand- 
water  system  was  recorded  and  the  top  electrode  of  the  cell  was  then 
secured.   The  resistance  measurements  were  taken  using  at  least  three 
separate  configurations  of  electrodes.   Ten  measurements  of  each  con- 
figuration were  taken  and  the  average  value  was  multiplied  by  its  ap- 
propriate geometric  cell  constant  to  give  the  sediment  resistivity. 
Electrode  configuration  resistivity  measurements  were  compared  with 
each  other  and  if  variations  greater  than  3  percent  existed,  the 
test  was  re-run.   On  completion  of  the  first  series  of  tests  the 
cell  was  tapped  approximately  twenty  times  and  the  sand  densified 
into  a  closely  packed  configuration.   Sediment  resistivity  measure- 
ments were  again  taken  and  the  length  of  the  sand-water  system  was 
recorded.   On  completion  of  these  tests,  the  top  electrode  was  re- 
moved and  the  temperature  again  recorded.   Temperature  variations 
were  seldom  more  than  a  degree  Centigrade,  since  the  sand  and  solu- 
tion were  usually  at  room  temperature  prior  to  mixing.   The  sand 
matrix  was  then  dried  at  105  C  for  24  hours  and  the  weight  of  sand 


was  obtained  for  porosity  measurements. 
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Clays  and  silt  were  tested  in  a  similar  manner  in  the 
simple  cell  apparatus,  with  the  exception  that  the  clays  and  silt 
were  allowed  to  equilibrate  24  hours  in  the  respective  salt  solu- 
tion under  a  partial  vacuum  to  remove  entrapped  air.   Slurries  were 
then  added  to  the  cell  and  allowed  to  settle  under  a  partial  vacuum 
in  order  to  avoid  formation  of  air  bubbles.   The  settling  period 
varied  depending  on  the  individual  sediment.   Approximately  six  or 
seven  slurry  replenishments  with  total  time  periods  of  24  to  36  hours 
for  settling  were  required  to  make  a  cell  ready  for  testing.   Due  to 
the  long  settling  periods  the  clay  measurements  were  only  made  for 
7.0,  17.5  and  35.0  parts  per  thousand  sodium  chloride  solution.   In 
addition  densif ication  of  the  highly  porous  clays  and  silt  were  not 
performed  since  tapping  the  cell  caused  the  resuspension  of  the  sedi- 
ment.  The  procedure  previously  explained  for  sand  samples  was  used 
to  measure  electrical  resistivity  and  porosity  of  clays  and  silt. 
The  total  time  to  run  sand  tests  was  approximately  twenty  minutes. 
Although  the  actual  testing  time  for  clays  and  silt  was  ten  minutes, 
a  test  cell  was  occupied  for  an  average  period  of  24  hours  due  to 
the  long  settling  periods. 

A  natural  marine  sediment  was  studied  in  the  early  stages  of 
this  investigation  in  order  to  (a)  determine  a  sediment  resistivity 
profile  versus  porosity  and  water  content,  (b)  determine  the  Forma- 
tion Factor  by  measuring  or  estimating  the  resistivity  of  interstitial 
water,  and  (c)  gain  experience  for  more  extensive  future  testing.   As 
previously  mentioned,  the  first  marine  sediment  was  obtained  by  both 
a  gravity  corer  and  the  DOSP.   Certain  sections  of  these  cores  had 
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previously  been  segmented  into  four  inch  sections  for  storage. 
Sixteen  of  these  short  segmented  sections  were  tested  in  the 
acrylic  resistivity  cell,   1  3/8  inch  diameter  samples  were  taken. 
The  electro-osmosis  sample  cutter  was  used  in  order  to  reduce  sample 
disturbance  and  core  shortening.   The  samples  were  transferred  to  a 
plastic  tube  of  the  same  dimensions  and  sealed  with  rubber  stoppers. 
This  tube  was  then  placed  in  a  constant  temperature  bath  (25  C)  for 
24  hours  prior  to  testing. 

In  order  to  determine  water  content  a  sample  was  taken  from 
the  excess  sediment.   Approximately  the  same  amount  of  sample  was 
placed  in  a  test  tube  and  diluted  with  100  ml  of  distilled  water. 
The  tube  was  sealed  and  the  mixture  was  shaken  vigorously.   After 
the  slurry  in  the  tube  settled,  a  conductivity  measurement  of  the 
solution  was  taken  using  the  Beckman  conductivity  meter.   If  it  is 
assumed  that  the  number  of  ions  in  solution  remained  constant  and 
conductivity  changes  are  linear  for  small  dilution  changes,  it  is 
possible  to  approximate  the  original  conductivity  of  the  interstitial 
water.   In  addition,  when  a  small  quantity  of  interstitial  water  was 
available,  an  American  Optical  Corp.  salinity  T/C  refractometer 
(Model  10402)  was  used  to  measure  salinity.   This  permitted  the 
determination  of  the  resistivity  of  the  interstitial  water. 

After  24  hours  the  marine  sediment  sample  was  removed  from 
the  temperature  bath  and  placed  in  the  testing  cell.   Three  resisti- 
vity measurements  were  taken  for  each  sample.   These  measurements 
corresponded  to  a  segment  of  the  sample  which  was  dissected  upon 
removal  of  the  sample  from  the  testing  cell.   The  individual  porosity 
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was  then  obtained  for  each  sample. 

Those  marine  samples  not  dissected  into  small  segments  were 
analyzed  directly  in  their  original  plastic  corer  liner.   The  stain- 
less steel  tube  containing  the  potential  electrodes  was  guided 
through  the  hole  on  the  top  current  electrode  and  vertically  pushed 
down  through  a  core  sample.   Electrical  resistivity  versus  depth  for 
the  entire  length  of  the  core  was  obtained.   Since  the  tubing  passes 
along  the  plastic  corer  liner  interior  wall,  negligible  sample  dis- 
turbance resulted  in  the  central  area  of  the  core  samples.   This  al- 
lowed electrical  resistivity  measurements  to  be  obtained  without 
destroying  the  sediment  sample  for  other  engineering  soil  tests. 

Modified  apparatus 

The  procedure  used  to  test  the  effect  of  consolidation  pres- 
sure on  clays  and  silts  with  the  modified  apparatus  was  as  follows: 
A  layer  of  Ottawa  sand  was  placed  in  the  four  liter  beaker  chamber. 
This  layer  allowed  the  load  applied  by  the  loaded  cells  to  be  equal- 
ly distributed  over  the  base  of  the  beaker.   The  support  frame  con- 
taining the  four  cells  was  placed  on  the  layer  of  sand,   A  micro- 
pore filter  was  placed  in  the  bottom  of  each  cell.  The  chamber  was 
then  filled  with  a  solution  of  known  normality,   A  slurry  was  pre- 
pared in  the  manner  previously  described,  with  a  solution  of  the 
same  normality  as  the  chamber  solution.   Approximately  50  ml  of  the 
slurry  was  introduced  to  the  cells.   After  the  slurry  settled,  a  new 
charge  was  added.   This  process  was  continued  until  all  four  cells 
in  the  chamber  were  filled.   A  micro-pore  filter  was  placed  on  top 
of  the  sediment  in  each  cell  and  the  piston  and  loading  rod  was 
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placed  on  the  sediment  giving  an  initial  pressure  of  about 

2 
0.01  kg/cm  . 

The  initial  vertical  length  between  the  loading  plate  and  the 
support  frame  was  measured.   The  weight  of  the  loading  plate  de- 
creased the  length.   After  two  or  three  days  the  movement  stopped. 
The  change  in  the  length  indicated  sediment  consolidation.   The  load 
was  then  doubled  and  the  procedure  repeated.   The  loads  used  re- 
sulted in  pressures  of  0.018,  0.036,  0.073,  0.110,  0.146,  0.292, 

2 

0.366,  0.438  and  0.511  kg/cm  .   Electrical  resistivity  measurements 

and  changes  in  lengths  were  taken  daily. 

Four  to  six  weeks  were  required  to  complete  a  chamber  test. 
During  this  time  period  the  test  apparatus  was  placed  under  a  clear 
plastic  hood  in  order  to  prevent  evaporation  which  would  change  the 
salinity  of  the  chamber  water  (see  Figure  13).   A  foam  rubber  mat, 
saturated  with  water,  was  placed  under  the  test  chamber.   A  humidity 
dial  was  placed  inside  the  hood.   A  thermometer  was  located  in  the 
beaker  chamber  in  order  to  allow  temperature  corrections  to  be  made. 
At  the  end  of  the  tests  the  cells  were  lifted  from  the  chamber  and 
the  sediment  height  in  each  cell  was  recorded.   All  of  the  sediment 
was  then  extruded  from  the  cell  cylinder  and  analyzed  for  solid  soil 
weight,  water  content  and  porosity. 

Redesigned  apparatus 

For  the  redesigned  equipment  a  calibration  curve  was  obtained 
by  filling  the  cell  with  standard  solutions  of  sodium  chloride  of 
known  conductivity  and  measuring  the  signal  conditioner  output  volt- 
age.  Sodium  chloride  of  solutions  of  0.1,  0.15,  0,2,  0.5  and  0.6 
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normality  were  used  to  obtain  calibration  curves.   Flow  patterns, 
cell  constants  and  calibration  curves  are  given  in  Appendix  B. 

The  method  used  to  test  clays,  silt  and  marine  sediments  in 
the  redesigned  test  apparatus  was  as  follows:   a  slurry  was  prepared 
in  the  manner  previously  described.   Water  from  the  slurry  was  ex- 
tracted and  placed  in  the  cell  to  determine  an  interstitial  water 
conductivity.   The  water  was  removed  from  the  cell,  the  slurry  was 
poured  into  the  cell  and  its  height  was  recorded.   The  slurry  was 
analyzed  for  water  content.   The  vessel  was  then  assembled  and  pres- 
sure was  applied.   All  sediments  were  tested  at  four  different  pres- 
sures:  0.25,  1.00,  4.00,  and  16.00  psi.   Each  pressurized  test  ran 
approximately  12  hours  or  until  consolidation  at  the  pressure  tested 
had  ceased.   The  pressure  vessel  was  disassembled.   The  sample  height 
and  the  volume  of  water  in  the  bottom  reservoir  was  measured. 

The  sample  was  extruded  and  analyzed  for  water  content  and 
porosity.  The  electrical  conductivity  of  the  sediment  sample  was 
continuously  recorded  during  the  entire  test  period.  Sample  calcu- 
lations for  porosity  and  water  content  are  presented  in  Appendix  C, 


RESULTS  AND  DISCUSSION 

The  conductivity  of  kaolinite  clay,  illite  clay,  Providence 
silt,  Ottawa  sand  (a)  and  angular  glacial  sand  (1-a)  were  plotted 

against  the  conductivity  of  the  interstitial  water  (see  Figure  15). 

12 
The  results  from  Kermabon  et  al.  '  for  a  clay  having  a  porosity  of 

50  percent  is  also  presented  in  this  figure.   For  each  of  the  sedi- 
ments a  linear  relation  was  plotted  by  the  method  of  least  squares. 
The  slope  of  these  lines  is  FF,  the  Formation  Factor.   The  porosity 
for  each  sediment  is  given  on  Figure  15.   The  variation  of  porosity 
for  each  type  of  sediment  was  within  +  2  percent,  and  this  is  probably 
the  cause  of  deviation  from  the  least  square  line  in  the  larger  grained 
sediments.   In  the  clays  this  deviation  could  also  be  influenced  by 
the  cation  exchange  capacity.  The  effect  that  ions  on  the  clay  sur- 
face have  on  the  Formation  Factor  is  illustrated  by  the  results 
presented  in  Figure  16.   In  this  figure  the  Formation  Factor  for  the 
three  clays  used  in  laboratory  experiments  (kaolinite,  illite,  and 
montmorillonite)  were  plotted  against  increasing  interstitial  water 
conductivity.   Though  porosity  was  held  constant  the  curves  show 
that  the  Formation  Factors  of  sodium  montmorillonite  and  illite  clays 
increased  with  increasing  interstitial  water  conductivity  from 
1.1  to  3.0  mho/meter.  From  this  point  on  the  Formation  Factor  be- 
came relatively  constant.   This  result  differs  from  the  result  found 
by  Atkins  and  Smith;  in  their  work  the  cation  exchange  capacity  of 
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Figure  15.   Electrical  Conductivity  of  the  Sediment  Versus 
the  Conductivity  of  the  Interstitial  Water  for 
Tllito  clay,  Kaolinite  clay,  Providence  silt, 
Ottawa  sand,  and  Glacial  sand. 
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Figure  16.   Formation  Factor  Versus  the  Conductivity  of  the 
Interstitial  Water  for  Kaolinite,  Illite  and 
Montmorillonite  Clay. 
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sodium  montmorillonite  was  of  no  influence  on  resistivity  at  an 
interstitial  water  salinity  of  1.1  mho/meter.   Although  the 
kaolinite  clay  data  presented  in  Figure  16  had  more  scatter,  it  is 
apparent  that  the  Formation  Factor  for  this  clay  remained  fairly 
constant  over  the  range  of  interstitial  water  conductivities  used 
in  these  tests  (1.1  to  5.6  mho/meter).   These  results  agree  with  the 
conclusion  made  by  Berg  that  the  effect  of  cation  exchange  capacity 
on  Formation  Factor  was  negligible  for  kaolinite  clay  above  inter- 
stitial water  conductivities  of  1.25  mho/meter. 

Sheeler  et  al.   stated  that  the  effect  of  cation  exchange 
capacity  on  Formation  Factor  persisted  up  to  interstitial  water  con- 
ductivities of  approximately  6  mho/meter  for  kaolinite  clay,  and 
11.8  mho/meter  for  montmorillonite  clay.  However  examination  of  the 
data  on  kaolinite  clay  presented  by  Sheeler  et  al.  indicates  that 
above  an  interstitial  water  conductivity  of  1.6  mho/meter  the  rate  of 
increase  of  the  Formation  Factor  is  so  small  that  it  could  be  consid- 
ered constant.   In  the  case  of  the  montmorillonite  clay  the  low 
porosities  could  account  for  the  increase  in  Formation  Factors  with 
interstitial  water  conductivities  up  to  11.8  mho/meter.  The  poro- 
sities of  clays  used  by  Sheeler  et  al.  ranged  from  48  to  67  percent, 
while  the  porosities  used  in  this  investigation  ranged  from 

87  to  94  percent. 

9 

Although  Doyce  offered  cation  exchange  capacity  as  a  pos- 
sible explanation  of  the  variation  in  his  data  (in  which  the  lowest 
interstitial  water  conductivity  encountered  in  the  cores  was  ap- 
proximately 5  mho/meter),  the  results  of  the  present  investigation 
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indicate  that  above  3  mho/meter  the  effect  of  cation  exchange  capa- 
city on  Formation  Factor  is  negligible  for  sediments.   Five  mho/meter 
is  equivalent  to  a  salinity  of  32.6  parts  per  thousand  at  25  degrees 
Centigrade,  while  3  mho/meter  is  equivalent  to  18.5  parts  per  thou- 
sand.  Although  no  cation  exchange  capacity  values  were  given  by 

12         13 

Kermabon  et  al.   and  Smith,   after  conducting  laboratory  experiments 

on  marine  clays  having  porosities  of  approximately  50  percent,  both 
concluded  that  there  was  no  appreciable  change  in  Formation  Factor 
when  using  interstitial  water  conductivities  ranging  from  1.4  to 
13  mho/meter. 

Another  interesting  result  of  the  present  investigation  is 
that  below  1.5  mho/meter  the  Formation  Factor  for  montmorillonite 
clay  was  less  than  one.  This  means  the  sediment  has  become  more 
conductive  than  the  interstitial  water.   This  was  also  noted  by 

o 

Sheeler  et  al.  when  working  with  montmorillonite  clay  of  lower  po- 
rosities.  This  fact  would  be  extremely  important  when  working 
with  montmorillonite  clay  with  low  interstitial  water  conductivities, 
since  a  correction  factor  would  have  to  be  developed  before  Formation 
Factor  could  be  considered  an  accurate  index  of  porosity. 

The  Formation  Factor  and  the  porosity  data  of  a  specific 
sediment  are  given  in  Appendix  C.   A  relationship  using  Equation  32 
was  developed.   If  Equation  32  is  re-written  the  following  relation- 
ship is  obtained: 

Log  FF  =  log  a  -  m  log  n.  (32) 
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From  the  above  relationship  it  is  evident  that  (m)  is  actually  the 
slope  of  the  line  and  (a)  is  the  value  of  the  Formation  Factor  at  a 
porosity  of  100  percent  or  the  intersection  of  the  y-axis  of  a 
log-log  plot.   Using  a  linear  regression  technique  these  coefficients 
(m  and  a)  were  determined  by  the  least  squares  method.   An  example 
of  a  complete  data  scatter  plot  relating  Formation  Factor  to  poros- 
ity for  kaolinite  clay  is  presented  in  Figure  17.   The  dashed  lines 
in  Figure  17  are  the  +  2  percent  error  lines  located  from  the  line 
of  best  fit. 

Data  for  all  the  sediments  tested  in  this  investigation  are 
plotted  as  Formation  Factor  versus  porosity  in  Figure  18.   The  lines 
of  best  fit  are  drawn  for  each  sediment  tested:   the  sands,  silts, 
clays  and  marine  sediments.   The  length  of  the  line  depends  on  the 
maximum  and  minimum  values  tested.   This  figure  represents  a  graphic- 
al method  of  predicting  porosity  from  Formation  Factors  obtained  by 
electrical  resistivity  measurements.  The  ranges  of  porosity  and 
Formation  Factor  obtained  are  presented  in  Table  5. 


TABLE  5 

THE  POROSITY  RANGE  AND  FORMATION  FACTOR  RANGE  OF  THE 
SANDS,  SILT,  CLAYS,  AND  MARINE  SEDIMENTS  TESTED 


SEDIMENT  TYPE 

POROSITY  RANGE  (7.) 

FORMATION  FACTOR 
RANGE 

Sands 

26  -  48 

3.1  -  5.9 

Silt 

41  -  62 

2.6  -  5.3 

Clays 

41  -  93 

1.1  -  3.1 

Marine  sediments 

.44  -  83 

1.6  -  5.7 
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Figure  17.      Formation  Factor-Porosity  Relationship 
for  Kaolinite  Clay. 
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Figure  18.      Formation  Factor-Porosity   Relationship 
for  all  Sediments  Tested. 
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The  following  data  for  each  sediment  tested  is  presented  in 
Table  6:   number  of  data  points  used  in  the  statistical  analysis, 
slope,  intercept,  standard  error  of  estimate,  percentage  of  data 
points  that  fall  within  +  2  percent  error  lines,  and  the  maximum 
and  minimum  Formation  Factor  and  porosity  values.   The  techniques 
used  in  this  statistical  analysis  are  discussed  in  Appendix  C. 

Because  of  the  ranges  of  Formation  Factor  and  porosity  (see 
Figure  18)  the  discussion  will  be  facilitated  by  dividing  the 
comparisons  into  the  following  three  groups: 

a)  sands  and  silt, 

b)  clays  and  silt, 

c)  clays,  silt  and  marine  sediments. 

Sands  and  Silt 

Formation  Factor  versus  porosity  is  plotted  for  the  five  sands 
and  Providence  silt  in  Figure  19.   Formation  Factor  for  the  sands  and 
silt  in  this  figure  ranged  from  3.12  to  5.90,  while  porosity  ranged 
from  26  to  48  percent.   The  Formation  Factors  for  Ottawa  sands  and 
glacial  sands  (1-a)  and  (1-b)  were  very  similar  at  the  different  poros- 
ities, while  the  more  naturally  occurring  sediments  (Providence  silt, 
glacial  sand  (2)  and  Narragansett  Bay  sand  (Station  C))  were  each 
uniquely  different.   Narragansett  Bay  sand,  although  actually  a  marine 
sediment,  was  considered  under  sands  and  silts  since  97  percent  of  the 
particles  were  in  the  sand  size  range. 

The  lines  of  best  fit  for  glacial  sand  (2)  and  Providence  silt 
at  porosities  of  41  ,to  46  percent  visually  suggest  that  electrical 
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Figure   19.      Formation  Factor-Porosity    Relationships   for-  Five 
Sands  and   Providence  Silt. 


76 


resistivity  decreases  with  increasing  particle  size.   Providence  silt 
which  has  a  mean  diameter  of  0.016  mm  had  a  higher  Formation  Factor 
than  glacial  sand  (2)  which  has  mean  diameter  of  0.50  mm.   At  a 
porosity  of  43  percent,  the  Formation  Factors  were  4.9  and  3.7 
respectively.   The  distributions  of  these  two  sediments  were  similar. 
As  discussed  later,  the  shape  of  the  silt  was  similar  to  the  sand, 
that  is  angular.   In  addition  to  particle  size,  structure  was  also 
variable  in  these  two  sediments.   At  43  percent  porosity  Providence 
silt  is  densely  packed  and  glacial  sand  (2)  is  loosely  packed.   How- 
ever the  results  discussed  later  show  that  structure  of  sands  does 
not  influence  electrical  resistivity  measurements  independent  of 
porosity.   Consequently  it  appears  that  particle  size  did  account 
for  the  difference  in  the  Formation  Factor  of  glacial  sand  (2)  and 
Providence  silt. 

These  results  are  in  agreement  with  Meredith,  who  reported 
that  mixtures  of  particles  that  were  smaller  in  size  by  10  diameters 
resulted  in  slightly  higher  Formation  Factors  than  mixtures  of  larger 
particles.   The  diameters  of  glass  beads  used  in  Meredith's  work 

ranged  from  6  mm  to  0.02  mm  while  the  size  of  sand  was  not  given, 

14 
The  data  of  Wyllie  and  Gregory   agreed  that  the  effect  of  different 

sized  glass  spheres  on  Formation  Factor  was  noticeable  only  when  the 
difference  in  particle  size  was  large.   The  diameters  of  the  glass 
spheres  used  in  their  study  ranged  from  3  mm  to  0.033  mm. 

An  additional  size  test  was  conducted  in  order  to  see  if  small 
differences  in  particle  size  of  a  sediment  could  be  detected.   Elec- 
trical resistivity  measurements  of  Ottawa  sand  (a)  with  a  diameter  of 
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0,70  mm  and  Ottawa  sand  (b)  with  a  diameter  of  0.42  mm  resulted  in 
Formation  Factors  of  3.6  and  3.5  respectively.   The  porosities, 
particle  shapes,  and  standard  deviations  of  these  two  sands  were 
approximately  the  same.   The  work  of  Meredith  on  glass  beads  agrees 
with  this  work  on  sediments  that  the  effect  of  particle  size  on 
conductivity  is  negligible  when  the  size  difference  is  small. 

The  effect  of  distribution  on  electrical  resistivity  can  be 
seen  in  Figure  19  by  comparing  Narragansett  Bay  sand  with  glacial 
sand  (1-a).   The  Formation  Factor  of  Narragansett  Bay  sand  which  had 
a  larger  distribution  was  higher  than  the  Formation  Factor  of  glacial 
sand  (1-a).   For  example  at  a  porosity  of  40  percent  the  Formation 
Factors  were  A. 36  and  3.85  respectively.  The  standard  deviation  of 
the  particle  size  for  Narragansett  Bay  sand  was  1.10  while  glacial 
sand  (1-a)  was  0.34.   Particle  size  and  shape  were  approximately  the 
same  (see  Table  3).  These  results  agree  with  the  work  of  Semenov 
on  sedimentary  rocks,  and  the  work  of  Meredith  on  glass  beads. 
Meredith  found  that  Formation  Factor  was  higher  for  a  mixture  of 
spheres  6  mm  and  0.75  mm  in  diameter  than  for  spheres  6  mm  in  dia- 
meter. 

At  lower  porosities  the  effect  of  distribution  was  found  to 
be  even  more  pronounced.   At  a  porosity  of  32  percent  the  Formation 
Factor  of  Narragansett  Bay  sand  was  5.70  and  glacial  sand  (1-a)  was 
4.77.   Therefore  the  difference  between  the  two  was  0.51  at 
40  percent  porosity  and  increased  to  0.93  at  32  percent  porosity. 
The  increased  assortment  of  particle  sizes  causes  the  current  to  take 
a  more  tortuous  path  and  thus  the  sediment  is  more  resistive.  As 
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the  sediment  becomes  densified  (lower  porosity),  the  variety  of 
particle  sizes  has  an  even  greater  effect  on  the  tortuosity. 

The  Narragansett  Bay  sand  had  a  small  quantity  of  platy  mica 
particles  which  could  be  thought  to  account  for  some  of  the  dif- 
ference between  Narragansett  Bay  sand  and  glacial  sand  (1-a).   How- 
ever as  discussed  later  in  this  paper,  particle  shape  differences  in 
sediments  do  not  seem  to  significantly  influence  electrical  resisti- 
vity measurements. 

Figure  20  a  to  f  are  photographs  of  Ottawa  sand,  glacial 
sands  1-a,  1-b,  2,  Narragansett  Bay  sand  and  Providence  silt  respec- 
tively. An  evaluation  of  angular  and  spherical  shapes  as  an  influ- 
ence on  electrical  resistivity  can  be  made  by  comparing  Ottawa  sand 
with  glacial  sand  (1-a).   Although  Ottawa  sand  is  spherical  in  shape 
and  the  glacial  sand  (1-a)  is  angular  in  shape,  the  Formation  Factor 
versus  porosity  relationship  is  very  similar.  At  a  porosity  of 
35  percent  Ottawa  sand  had  a  Formation  Factor  of  4.20  while  glacial 
sand  (1-a)  had  4.35.   Size,  distribution  and  structure  were  approxi- 
mately the  same  for  these  two  sands.  On  Figure  19  95  percent  of  the 
data  points  lie  within  +  2  percent  of  the  line  of  best  fit.   Conse- 
quently the  difference  between  the  line  of  best  fit  for  Ottawa  sand 
and  glacial  sand  (1-a)  is  not  significant.   These  results  indicate 
that  electrical  resistivity  is  quite  similar  for  angular  and  spheri- 
cal shapes. 

14 
Wyllie  and  Gregory   obtained  a  Formation  Factor  of  4.2  for 

glass  spheres  and  spherical  Ottawa  sand,  4.4  for  angular  beach  and 

creek  sands,  and  4.7  for  lucite  triangular  prisms  at  35  percent 
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Figure  20-a.      Ottawa  Sand. 
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Figure  20-b.     Glacial  Sand    (1-a) . 
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Figure  20-c.      Glacial  Sand    (1-b) . 
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Figure  20-d.      Glacial  Sand    (2) 
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Figure  20-e.      Narragansett  Bay  Sand. 
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Figure  20-f.  .Providence  Silt. 
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porosity.   The  sphericity  was  1,0,  0.8,  and  0,6  respectively.   The 
Formation  Factor  in  their  work  seems  to  increase  as  sphericity  de- 
creases.  However  the  particle  size  was  not  held  constant.   The 
particle  size  of  the  glass  spheres  was  3  mm  in  diameter  and  the  tri- 
angular prisms  were  0,25  mm  in  diameter.   As  discussed  previously 
this  difference  in  particle  size  is  large  enough  to  influence  the 
resistivity.   The  Formation  Factor  would  be  higher  for  the  triangu- 
lar prisms  since  the  diameter  was  12  times  smaller  than  the  glass 
spheres. 

The  influence  of  structure  on  the  electrical  resistivity  of 
sands  was  tested  by  making  resistivity  measurements  before  and  after 
densif ication.   The  electrical  resistivity  and  porosity  were  changed 
by  densif ication  in  the  same  proportion.   The  Formation  Factors  of 
the  five  sands  at  40  percent  and  32  percent  porosities,  the  percent- 
age change  in  the  Formation  Factor  and  the  measure  of  spread  are 
given  in  Table  7,   For  a  20  percent  decrease  in  porosity  the  Forma- 
tion Factor  increased  19  to  24  percent.  As  shown  in  Table  7,  the 
Formation  Factors  of  the  two  sand  types  with  a  large  distribution 
changed  a  little  more  than  the  sands  with  a  small  spread  in  particle 
size. 

Tortuosity  factors  as  defined  by  Equation  13  are  given  in 
Table  8  for  loosely  and  densely  packed  sands.   Tortuosity  increases 
as  density  increases  for  those  sands  that  have  a  large  variety  of 
particle  sizes.   Densif ication  does  not  seem  to  affect  the  tortuosity 
of  those  sands  with  a  small  spread  in  particle  size. 


TABLE  7 

THE  FORMATION  FACTORS  AT  40  PERCENT  AND  32  PERCENT  POROSITY, 
THE  PERCENTAGE  CHANGE  IN  FORMATION  FACrOR  AND  THE 
MEASURE  OF  SPREAD  FOR  FIVE  SANDS. 
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Name 

(1) 

FF 
40% 

FF 
32% 

%  FF  (2) 
increased 

(3) 

H 

Ottawa  sand 

3.73 

4.60 

19 

0.25 

Glacial  sand  (1-a) 

3.85 

4.77 

19 

0.34 

Glacial  sand  (1-b) 

3.73 

4.60 

19 

0.50 

Glacial  sand  (2) 

3.96 

5.00 

21 

1.31 

Narragansett  Bay 
sand  (Station  C) 

4.35 

5.70 

24 

1.10 

(1)  FF  _,  -  Formation  Factor  at  40%  porosity. 

FF    —  FF 

(2)  %  FF  increased  =»   32%  "   40%  (100). 

FF32% 

(3)  6.  (Dispersion  measure)  =  £  (0Q4  -  0,6),  where  0  .  and  0   are 
cumulative  percentages  of  the  sample  which  is  coarser  grained  than  84 
and  16  percentiles  respectively.  The  0  scale  is  defined  as  -leg, 
(diameter  in  millimeters). 
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TABLE   8 


COEFFICIENTS  OF  PERMEABILITY,  FORMATION  FACTORS,  POROSITY,  TORTUOSITY, 
AND  THE  PRODUCT  OF  THE  FORMATION  FACTOR  AND  THE  COEFFICIENT  OF  PERMEA- 
BILITY  FOR  FIVE   LOOSELY   AND  DENSELY    PACKED   SANDS. 

Name  of  Sediment     CooCf iciont  Formation     Porosity     Tortuosity     KpFF 

of  Perinea-  Factor  (n)     (%)         y/FF;n 

bility   K_,  (FF) 
(cm/sec ) 


Ottawa  Sand 
(loose) 

0.239 

Ottawa  Sand 
(dense) 

0.154 

Glacial  Sand 
(1-a) (loose) 

0.0891 

Glacial  Sand 
(1-a) (dense) 

0.0502 

Glacial  Sand 
(1-b) (loose) 

0.0224 

Glacial  Sand 
(1-b) (dense) 

0.00699 

Glacial  Sand 
(2)  (loose) 

0.115 

Glacial  Sand 
(2)  (dense) 

0.0413 

Narragansett  Bay 
Sand  (loose) 

0.0879 

Narragansett  Bay 
Sand  (dense) 

0.0248 

3.55 


5.30 


3.43 


5.00 


3.32 


4.80 


3.50 


5.85 


4.30 


5.75 


42.0 


27.5 


45.0 


30.2 


44.0 


31.2 


44.0 


27.2 


40.5 


31.8 


1.22 


1.20 


1.24 


1.23 


1.21 


1.22 


1.24 


1.26 


1.32 


1.35 


0.848 


0.816 


0.306 


0.251 


0.0744 


0.0336 


0.403 


0.252 


0.378 


0.143 
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Coefficients  of  permeability  (Kp)  and  the  corresponding 
Formation  Factor  and  porosity  are  also  presented  in  Table  8.   The 
coefficient  of  permeability  as  obtained  by  a  variable  head  perme- 
ability test,  was  reduced  as  the  sands  densified  and  Formation 
Factors  increased.   The  coefficient  of  permeability  decreased  as 
the  mean  particle  size  decreased.   It  decreased  more  for  densified 
sands  with  a  large  distribution  of  particle  sizes.   As  the  porosity 
of  a  sand  decreased,  the  coefficient  of  permeability  decreased 
logarithmically  and  the  Formation  Factor  increased  logarithmically. 
Due  to  the  cancelling  effect  of  these  two  logarithmic  variations, 
the  product  of  Formation  Factor  times  coefficient  of  permeability  of 
densely  and  loosely  packed  sand  can  be  plotted  on  a  linear  scale. 
This  linear  assumption,  along  with  a  Formation  Factor  versus  poros- 
ity plot  allow  the  coefficient  of  permeability  to  be  obtained  by 
measuring  only  the  Formation  Factor.   Thus  by  knowing  any  one  of  the 
three  variables  (FF,  n,  Kp),  the  other  two  variables  can  be  predicted 
from  standard  curves.   The  product  of  the  Formation  Factor  and  the 
coefficient  of  permeability  (FFKp)  is  presented  for  all  the  sands  on 
Table  8.   The  relationship  between  porosity  and  FFKp  for  loosely  and 
densely  packed  sands  is  plotted  in  Figure  21.   By  assuming  a  poros- 
ity of  36  percent  and  using  Figures  19  and  21,  the  FF  FFKp  and  the 
Kp  for  the  five  sands  used  in  these  tests  are  given  in  Table  9. 

The  sands  that  are  more  common  in  nature  have  the  greater  slopes 
in  Figure  21.   For  these  sands  a  small  change  in  porosity  results  in  a 
larger  change  in  FFKp  than  for  the  sorted  sands.   Consequently  these 
graphs  predicting  porosity  and  permeability  from  electrical  resisti- 
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Porosity  Versus  the  Product  of  the  Formation  Factor  and 
the  Coefficient  of  Permeability  (FFKp)  for  Five  Sands. 
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TABLE  9 

FORMATION  FACTORS,  COEFFICIENTS  OF  PERMEABILITY  AND  THE  PRODUCT  OF 
FORMATION  FACTOR  AND  THE  COEFFICIENT  OF  PERMEABILITY  AT  36  PERCENT 
POROSITY  FOR  FIVE  SANDS. 


Name 

FF36% 

FFKp  36% 
(cm/sec) 

5>  36% 
(cm/sec) 

Ottawa  Sand 

4.12 

.834 

.202 

Glacial  Sand  (1-a) 

4.26 

.274 

.0643 

Glacial  Sand  (1-b) 

4.12 

.048 

.0117 

Glacial  Sand  (2) 

4.40 

.321 

.0730 

Narragansett  Bay  Sand 

4.95 

.258 

.0521 
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vity  measurements  would  be  very  useful  in  a  natural  environment. 

Clays  and  Silt 

A  Formation  Factor-porosity  diagram  with  curves  of  best  fit 
for  kaolinite,  illite,  and  sodium  montmorillonite  clays  and 
Providence  silt  is  shown  in  Figure  22.   This  figure  includes  over 
200  measurements  made  with  all  three  versions  of  the  laboratory  test 
equipment.   The  statistical  data  for  this  figure  is  presented  in 
Table  6. 

The  results  of  this  investigation  on  clays  and  silts  suggest 
that  particle  shape  does  not  influence  electrical  resistivity  inde- 
pendent of  porosity.   Kaolinite  clay  (0.6N  NaCl)  and  Providence  silt 
have  almost  identical  Formation  Factors  over  a  porosity  range  of 
51  to  61  percent  although  the  former  is  platy  shaped  and  the  latter 
is  angular  in  shape.   These  shapes,  enlarged  2000  and  10,000  times 
by  an  electron  microscope,  can  be  seen  in  Figures  23-a,  23-b,  20-f 
and  24.   The  particle  distribution  for  the  two  sediments  was  similar; 
the  measure  of  spread  for  kaolinite  was  1.81  while  Providence  silt 
was  1.21.   The  particle  size  of  Providence  silt  was  approximately 
five  times  larger  than  kaolinite.   It  was  suggested  in  the  discussion 
on  results  of  sands  that  a  particle  size  difference  of  about  10  times 
is  necessary  before  size  difference  is  detectable  by  electrical  re- 
sistivity measurements.   In  the  results  and  discussion  on  cation  ex- 
change capacity  it  is  indicated  that  cation  exchange  capacity  would 
not  be  a  factor  influencing  the  electrical  resistivity  measurements 
of  these  two  sediments  since  the  salinity  was  high  (0.6  normal  sodium 
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Figure  22.  Formation  Factor-Porosity  Relationships  for 
Kaolinite,  Illite  and  Kontmorillonite  Clays 
and  Providence  Silt. 
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Figure  23.   Photographs  of  (a)  Kaolinite  Clay  Enlarged  2000 
a,b   Times  by  Electron  Microscope,  (b)  Kaolinite  Clay 
Enlarged  10,000  Times. 
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Figure  24.   Photograph  of  Providence  Silt  Enlarged  10,000 
Times  by  Electron  Microscope. 


95 


chloride).   Since  both  would  have  a  flocculant  structure  at  this 
salinity,  structure  would  not  be  a  variable.   Since  shape  is  the 
only  major  variable  between  the  two,  and  since  the  Formation  Factors 
are  nearly  identical  in  this  porosity  range,  it  appears  that  parti- 
cle shapes  such  as  platy  or  angular  do  not  influence  electrical  re- 
sistivity.  The  results  for  illite  clay  further  substantiate  this 
conclusion.   Although  illite  clay  is  of  a  different  mineral  content 
from  kaolinite,  the  Formation  Factors  were  again  almost  identical. 
This  result  indicates  that  this  difference  in  mineral  content  has 
little  if  any  effect  on  electrical  resistivity  measurements.   Illite 
clay  has  a  platy  shape  that  can  be  seen  in  Figure  25-a  and  25-b. 
The  mean  particle  size  of  illite  clay  is  about  the  same  as  for 
kaolinite  clay  (five  times  smaller  than  Providence  silt).   The  Forma- 
tion Factor  of  illite  clay  and  Providence  silt  is  nearly  identical 
though  shape  differs  and  appears  to  be  the  only  major  variable.   The 
opposite  is  true  of  kaolinite  clay  and  montmorillonite  clay.   Although 
these  two  have  a  very  similar  shape,  the  Formation  Factor  versus 
porosity  relationsbip  for  these  two  is  very  different  (Figure  26-a, 
and  26-b).   For  example  at  a  porosity  of  83  percent,  Formation  Fac- 
tor of  kaolinite  is  1.45  while  montmorillonite  is  2.00.   As  discussed 
later  this  difference  is  probably  due  to  particle  size  and  distribu- 
tion.  The  results  for  illite  and  montmorillonite  clays  reinforce  the 
conclusion  that  the  influence  of  particle  shape  on  electrical  resis- 
tivity measurements  of  sediments  appears  to  be  small. 

The  standard  deviation  of  particle  size  for  illite  clay  was 
2,99  and  might  appear  to  be  a  variable  when  compared  with  Providence 
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Figure  25.   Photographs  of  (a)  Illite  Clay  Enlarged  2000  Times 
a,b   by  Electron  Microscope,  (b)  Illite  Clay  Enlarged 
10,000  Times. 
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Figure  26.   Photographs  of  (a)  Montmorillonite  Clay  Enlarged 
a,b   2000  Times  by  Electron  Microscope,  (b)  Montmoril- 
lonite Clay  Enlarged  10,000  Times. 
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silt.   Although  this  standard  deviation  is  twice  that  of  Providence 
silt,  this  difference  is  probably  not  significant.   As  a  descriptive 
measure  of  spread,  standard  deviation  is  less  accurate  when  the 
sediment  grains  are  very  fine.   This  is  due  to  the  assumption,  in 
the  hydrometer  method  of  analysis,  that  all  particles  are  spherical 
and  do  not  aggregate  while  settling.   This  assumption  can  be  used 
with  reasonable  accuracy  for  most  particles  above  0.005  mm  in  size. 
However  below  0.005  mm  particles  vary  widely  in  shape  and  aggregate 
more  readily.   The  formula  for  the  standard  deviation  of  the  parti- 
cle size  6  .   is  based  on  the  84th  percentile  diameter  (see  Table  3), 
0 

Eighty  four  percent  of  the  Providence  silt  particles  were  coarser 
than  0.006  mm  while  84  percent  of  the  illite  clay  particles  were 
coarser  than  0.0006  mm.   It  is  apparent  that  error  could  have  been 
introduced  by  aggregation  and  non-spherical  particles  in  the  illite 
clay  analysis. 

Descriptive  measures  such  as  phi  median  diamter  (^c0)  or 
phi  deviation  measure  (6  .)  determined  by  a  standard  hydrometer  meth- 
od of  analysis  cannot  be  used  for  montmorillonite  clay.   Approximate- 
ly 80  percent  of  the  particle  sizes  are  below  0.005  mm  in  this  clay. 
Other  methods  of  analysis  such  as  optical  equipment  or  centrifugal 
methods  are  required.   However  time  and  cost  make  these  methods  im- 
practical.  Spread  or  distribution  can  be  estimated  from  the  per- 
centages of  sand,  silt  and  clay  that  make  up  the  sediment,   Mont- 
morillonite clay  which  has  a  large  percentage  of  clay  sized  parti- 
cles could  be  expected  to  have  a  large  spread  in  particle  sizes,  A 
large  distribution  would  be  expected  to  increase  Formation  Factor, 
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Montmorillonite  was  more  resistive  than  illite  and  kaolinite,  al- 
though size  probably  influenced  this  result  also.   As  discussed 
earlier  large  distributions  increased  Formation  Factor  for  sands 
when  mean  particle  size  was  constant.   Just  as  was  found  in  the 
results  on  sands  and  silts,  particle  size  appears  to  influence  elec- 
trical resistivity  measurement  of  clays  although  porosity  is  constant. 
As  stated  earlier  higher  Formation  Factors  were  obtained  for  mont- 
morillonite than  for  kaolinite.   One  of  the  two  variables  (in  addi- 
tion to  distribution)  is  particle  size.   The  mean  particle  diameter 

of  montmorillonite  is  reported  to  be  «  0.001  mm  while  kaolinite  is 

32 
0.003  mm.   According  to  Osthaus   the  particle  diameters  for  mont- 
morillonite clays   ranged  from  0.0016  to  0.00005  mm.   These  measure- 
ments were  made  using  a  super-centrifuge  operating  at  various  speeds 
up  to  50,000  rpm.   This  very  large  size  difference  in  particle  size 
most  likely  accounts  for  some  of  the  difference  in  the  Formation 
Factor  of  montmorillonite  and  kaolinite.   At  82  percent  porosity, 
the  Formation  Factor  of  these  two  clays  was  2.10  and  1.48  respec- 
tively.  It  is  doubtful  that  the  mineral  difference  of  these  two 
clays  influenced  the  Formation  Factors,  since  both  are  non-conductive 
and  the  interstitial  water  was  0.6  normal  sodium  chloride. 

Smaller  differences  in  mean  particle  size  (i.e.   Providence 
silt  (0.016  mm)  and  kaolinite  (0.003  mm))  did  not  appear  to  influ- 
ence the  electrical  resistivity  measurements.  As  seen  in  Figure  22 
the  Formation  Factors  were  nearly  identical  over  the  porosity  range 
of  51-61  percent  for  these  sediments.   Because  of  the  larger  sized 
particles,  Providence  silt  was  not  in  a  stable  condition  above 
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61  percent  porosity.   If  the  particle  size  is  smaller  (as  in 
kaolinite  clay)  the  sediment  exists  in  a  stable  condition  at  poros- 
ities up  to  84  percent.   Electrical  resistivity  measurements  were 
conducted  up  to  a  91  percent  porosity  by  producing  a  slurry.   The 
data  of  Atkins  and  Smith  '  for  kaolinite  and  sodium  montmorillonite 
clay  is  also  presented  in  Figure  22.   The  Formation  Factor-porosity 
data  agrees  well  with  the  results  obtained  in  this  investigation. 
However  in  their  work  the  difference  between  kaolinite  and  montmorillo- 
nite was  attributed  to  shape  and  values  for  a  given  shape  factor  (m) 
were  determined.   As  previously  discussed  the  results  of  the  present 
study  on  clay  sediments  suggest  that  the  difference  between  the  two 
is  due  to  particle  size  aind  distribution. 

Meredith,  after  studying  the  electrical  conductivity  of 
dispersions  (spheres  and  discs),  concluded  that  errors  as  high  as 
10  to  20  percent  were  due  to  shape  differences.   Semenov  found  that 
disc  shaped  particles  in  rock  resulted  in  higher  Formation  Factors 
than  rocks  of  spherical  particles  at  equal  porosities.   The  results 
of  the  preseiat  study  indicate  that  shape  is  not  a  factor  influencing 
the  electrical  resistivity  of  sediments.   The  earlier  discussion  on 
sands  and  Hilts  showed  tlhat  differences  between  spherical  and  angu- 
lar particles  were  not  detected  by  electrical  resistivity  measure- 
ments.  The  results  on  clays  and  silts  indicate  that  nearly  identi- 
cal electrical  resistivity  measurements  resulted  from  sediments  with 
angular  or  platy  particle;  shapes  as  the  only  appat'ent  variable. 

The  effect  of  structure  on  electrical  resistivity  can  be  seen 
in  Fig,ure  22  if  kaolinite  (0.6N  NaCl)  is  compared  with  kaolinite 
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using  0.2  molar  sodium  tetra  phosphate  solution  (Quadrafos)  as  the 
interstitial  water  of  the  sediment.   The  structure  naturally  en- 
countered in  a  saline  environment,  a  flocculant  structure  (see 
Figure  2),  resulted  in  a  lower  Formation  Factor  than  when  Quadrafos, 
a  common  def locculating  agent,  was  used  to  obtain  a  dispersed  struc- 
ture.  The  difference  between  the  two  increased  as  porosity  de- 
creased.  For  example  at  a  porosity  of  80  percent  the  Formation  Fac- 
tor for  kaolinite  (0.6N  NaCl)  was  1.33  while  kaolinite  (Quadrafos) 
was  1.55.   At  a  porosity  of  52  percent  the  Formation  Factors  were 
2.30  and  3.60  respectively. 

In  the  dispersed  structure  the  clay  particles  are  in  a 
preferred  orientation,  that  is  one  particle  tends  to  be  oriented  to 
the  next  particle  in  the  same  manner.   However  in  the  flocculant 

structure  the  clay  particles  are  in  a  randomly  oriented  pattern  to- 

33 
ward  each  other.   According  to  Brindley   preferred  orientation  will 

bring  out  or  increase  the  intensity  of  X-ray  diffraction  patterns  of 

basal  spacing  (d-spacing)  reflections.  The  d-spacing  is  the  spacing 

between  crystal  lattice  planes. 

A  side  test  was  conducted  to  check  if  the  particle  orienta- 
tion assumed  previously  could  be  validated  by  X-ray  diffraction 
techniques.   Two  X-ray  diffraction  slides  of  equal  thickness  were 
prepared  under  identical  conditions  for  kaolinite  clay  with 
0.6N  sodium  chloride  as  the  interstitial  water  and  kaolinite  clay 
with  0.02M  sodium  tetra  phosphate. 

Using  the  General  Electric  diffractometer  (model  XRD-5S)  in 
the  Metallurgical  Research  Laboratory  at  the  University  of  Rhode 
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Island,  each  sample  was  scanned  from  8  to  26  degrees.   The  proce- 
dure used  is  presented  in  Appendix  C.   This  scan  covered  the  two 
major  d-spacing  values  of  kaolinite  clay;  7.2  and  3,6  Angstroms. 
The  X-ray  dif fractograms  of  the  two  samples  are  presented  in 
Figure  27.   The  vertical  axis  is  a  logarithmic  indication  of  inten- 
sity corresponding  to  an  approximate  quantitative  index  of  the 
mineral  present.   Although  many  factors  affect  intensity  of  the  dif- 
fracted X-ray  beams,  orientation  of  particle  and  the  reinforcement 
of  intensity  caused  by  preferred  orientation  or  the  reduction  of  in- 
tensity caused  by  random  orientation  of  the  same  mineral  particles 
would  be  a  major  factor  determining  the  magnitude  of  intensity.  The 
difference  in  intensity  between  the  two  samples  is  illustrated  in 
Figure  27.  An  average  intensity  count  was  made  at  the  two  major 
d-spacing  values.  The  results  are  shown  in  Table  10. 

TABLE  10 

THE  AVERAGE  INTENSITY  COUNTS  OF  THE  TWO  MAJOR  D-SPACING 
VALUES  OF  TWO  STRUCTURES  OF  KAOLINITE  CLAY 


Sample  name 

7.2  d-spacing 
Intensity  (counts) 

3.6  d-spacing 
Intensity  (counts) 

Kaolinite  (0.6N  NaCl) 
Kaolinite  (Quadrafos) 

8758 
47333 

! 

8884 
47842 

The  results  show  that  the  kaolinite  clay  (Quadrafos)  reinforced  the 
intensity  when  compared  with  kaolinite  clay  (0.6N  NaCl).  These  re- 
sults indicate  that  kaolinite  clay  (Quadrafos)  was  of  a  dispersed 
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Figure  27.      X-Ray  Diffractografns  of  Two  Structures  of  Kaolinite 
Clay. 
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structure  while  kaolinite  clay  (0.6N  NaCl)  was  of  a  flocculant 
structure. 

The  salinity  of  marine  and  estuarine  environments  results  in 
clay  deposits  formed  by  f locculation.   Extremely  high  lateral  or 
vertical  pressure  or  perhaps  even  disturbance  could  cause  the  floc- 
culant structure  to  collapse  to  a  dispersed  structure.   These  results 
suggest  that  the  electrical  resistivity  measurements  would  change  as 
the  structure  changed  from  flocculated  to  dispersed. 

the  effect  of  consolidation  pressure  on  the  conductivity  of 
clay  and  silt  sediments  can  be  seen  in  Figure  28.   These  curves  il- 
lustrate that  the  electrical  conductivity  is  inversely  proportional 
to  the  applied  pressure  in  the  range  from  0.25  to  16.00  psi.   In  this 
range  the  Formation  Factor  increased  in  proportion  to  porosity  de- 
crease.  Consequently  consolidation  pressure  seems  to  have  no  effect 
on  electrical  resistivity  of  sediments  independent  of  porosity. 

Figure  28  presents  four  curves  showing  illite  clay  consoli- 
dating at  pressures  of  0.25,  1.00,  4.00,  and  16.00  psi  and  the  re- 
sulting conductivity  values  of  2.7,  2.5,  2.0,  and  1.4  mho/meter 
respectively.   Similar  results  obtained  for  Providence  silt,  kaolinite 
clay  and  montmorillonite  clay,  can  be  found  in  Appendix  C,   Porosity 
was  determined  at  the  beginning  and  end  of  the  consolidation  pres- 
sure tests* 

As  discussed  in  the  Review  of  the  Literature  several  studies 

18  19  20 
on  formation  rocks   '   '   have  reported  that  Formation  Factor  in- 
creases much  faster  than  the  decrease  in  porosity.   The  pressures 
used  (1,000  to  12,000  psig)  were  much  greater  than  the  16  psi  used 
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Figure  28.  Conductivity  Versus  Time  Curves  of  lllite  Clay  at 
Consolidation  Pressures  of  0.25,  1.00,  4.00  and 
16.00  psig. 
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in  this  study.   It  would  seem  possible  that  sediments  would  be  more 
easily  influenced  by  overburden  pressure  than  rocks.   However  the 
greater  increase  of  Formation  Factor  over  porosity  decrease  was  not 
found  on  sediments  at  the  low  consolidation  pressures  used  in  this 
study. 

A  typical  conductivity  versus  time  recording  obtained  by 
the  redesigned  equipment  is  presented   in  Figure  29  for  kaolinite 
clay  (0.6N  NaCl).   Most  tests  were  run  for  12  hours.   At  the  times 
designated  as  t.  and  t.  in  the  figure,  a  Formation  Factor  was  ob- 
tained and  compared  with  porosity  values  determined  on  the  initial 
slurry  and  after  consolidation.   These  two  values  were  compared  with 
Formation  Factor-porosity  values  for  kaolinite  (0.6N  NaCl)  from 
Figure  22.   At  t.  and  t«  the  Formation  Factors  were  1.64  and  3.27 
respectively.   These  values  fall  within  the  +  2  percent  porosity 
error  lines  for  kaolinite  (0.6N  NaCl)  in  Figure  22.   The  values  for 
conductivity  of  the  sediment  sample  were  changed  to  Formation  Factor 
by  dividing  the  conductivity  of  the  interstitial  water  by  the  con- 
ductivity of  the  sediment.   By  comparing  these  Formation  Factors 
with  those  of  Figure  22,  values  of  porosity  were  obtained  and  used 
on  the  vertical  scale  of  Figure  29.   Since  porosity  and  void 
ratio  (e)  are  uniquely  related  by: 


•  -  &  •  <"> 


the  void  ratio  was  included  as  an  addition  on  the  vertical  axis. 
For  saturated  sediments  a  relationship  exists  between  the  void  ratio 
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Figure  29.   Conductivity,  Formation  Factor,  Porosity,  Void 

Patio  and  Water  Content  Versus  Time  of  Kaolinite 
Clay  at  a  Consolidation  Pressure  of  16  psig. 
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and  water  content  (w)  that  can  be  obtained  by  knowing  the  specific 
gravity  of  the  sediment: 


w  =  -~  .  (48) 


The  specific  gravity  (G)  of  kaolinite  clay  is  2.60.  Water  content 
was  calculated  and  also  included  on  the  vertical  axis  of  Figure  29. 
Therefore  Figure  29  is  a  continuous  record  of  the  change  in  Forma- 
tion Factor,  porosity,  void  ratio  and  water  content  of  a  kaolinite 
clay  sediment  as  it  is  consolidated  by  pressure  of  16.00  psig  over  a 
twelve  hour  period.  The  total  test  time,  consolidation  pressure, 
and  for  times  t.  and  t~  the  Formation  Factors,  porosities,  void 
ratios,  and  water  contents  for  all  sediments  tested  by  the  redesigned 
equipment  are  presented  in  Appendix  C. 

The  continuous  conductivity  readings  by  the  redesigned  equip- 
ment can  be  extremely  useful.   If  accurate  Formation  Factor-porosity 
relationships  are  determined  for  a  sediment  as  they  were  in  this 
study,  then  from  the  continuous  conductivity  readings  (Formation  Fac- 
tor) the  porosity,  void  ratio,  and  water  content  can  be  obtained  at 
any  point  during  the  monitoring.  Using  the  void  ratios  obtained  at 
the  end  of  a  test,  a  void  ratio  versus  the  log  of  the  consolidation 
pressure  is  plotted  in  Figure  30.   In  classical  soil  mechanics  this 
type  of  curve  is  very  useful  in  predicting  consolidation.   If  a 
Formation  Factor-porosity  curve  of  a  sediment  is  sufficiently  accu- 
rate, then  the  redesigned  equipment  could  be  used  as  a  consolido- 
meter.  Curves  such  as  Figure  30  could  be  obtained  by  simply  contin- 
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Figure  30.      Formation  Factor,    Porosity  and  Void   Ratio  Versus 
the  Log  of  the  Consolidation  Pressure  of   Illite 
Clay. 
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ually  recording  the  electrical  conductivity  of  the  sample  and 
incrementally  increasing  the  pressure  after  set  periods  of  time 
without  removing  the  sample.   Then  by  using  the  Formation  Factor- 
porosity  relationship  the  porosity  and  consequently  void  ratio 
could  be  obtained  and  plotted  versus  the  corresponding  pressure. 
Currently  tests  are  being  conducted  with  a  more  advanced  teflon 
cell  in  the  redesigned  equipment.   This  cell  reduces  side  wall 
friction  and  allows  more  extensive  and  elaborate  testing  to  be 
conducted. 

Since  the  water  driven  out  of  the  sample  during  consolida- 
tion was  collected  in  the  reservoir  of  the  redesigned  equipment,  it 
was  possible  to  obtain  an  approximate  value  of  the  coefficient  of 
permeability  (K  ),   Thus  by  using  Darcy's  Law» 


where 


Kp-f.  i.A.  (49) 


K   is  the  coefficient  of  permeability  (cm/sec), 

q   the  quantity  of  water  collected  in  the  reservoir  during 

3 
the  test  (cm  ), 

t   the  total  time  of  the  test  (sec), 

L  the  average  length  of  the  sample  (obtained  by  measuring 
the  initial  length  of  sample  and  adding  it  to  the  final 
sample  length,  then  dividing  by  two)  (cm), 

Ah  the  pressure  used  in  the  test  (cm), 

2 

A   the  cross-sectional  area  of  the  sample  (cm  ), 


Ill 


An  example  of  the  results  of  this  type  of  calculation  is  presented 
in  Figure  31.   In  this  figure  the  void  ratio  at  the  end  of  a  test  is 
plotted  against  the  coefficient  of  permeability  for  the  illite  clay. 
From  this  curve  it  is  possible  to  obtain  other  coefficients  of  perme- 
ability from  Formation  Factors,  if  the  Formation  Factor-porosity 
relationship  has  been  developed. 

Marine  Sediments 

For  the  four  different  marine  sediments  tested 
90  to  100  percent  of  the  Formation  Factor-porosity  data  points  fell 
within  +  2  percent  of  the  line  of  best  fit.   This  indicates  that 
electrical  resistivity  measurements  can  accurately  predict  porosity 
of  marine  sediments.   The  values  that  fell  outside  the  +  2  percent 
error  lines  were  unusual  samples  usually  containing  pieces  of  coral 
or  shells. 

Two  vertical  profiles  of  a  marine  sediment  core  sample  ob- 
tained by  using  the  simple  laboratory  test  apparatus  are  presented  in 
Figures  32-a  and  32-b.   Porosity,  water  content,  and  sediment  resis- 
tivity (corrected  to  a  temperature  of  25  C)  are  shown.   The  points 
on  the  curves  represent  discrete  measurements.   The  horizontal  dashed 
line  indicates  where  the  core  was  cut  to  facilitate  handling  and 
storage. 

The  electrical  resistivity  of  this  sediment  was  inversely  pro- 
portional to  the  water  content  and  porosity.   An  exception  is  illus- 
trated in  Figure  32-a.  A  large  shell  approximately  two  inches  in 

diameter  was  located  at  a  depth  of  20  inches.  At  this  location  the 
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Figure  31.  Formation  Factor,  Porosity  and  Void  Ratio 
Versus  the  Coefficient  of  Permeability  of 
Illite  Clay. 
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Figure  32.  Two  Vertical  Depth  Profiles  of  Porosity,  Water 
a,b   Content  and  Electrical  Resistivity  of  a  Marine 
Sediment. 
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resistivity  of  the  sediment  increased  sharply,  but  the  water  content 
and  porosity  remained  fairly  constant.   These  sharp  changes  in  re- 
sistivity are  indicators  of  structure  and  fabric  changes  along  the 
core  sample.  One  1\   inch  diameter  core  showed  sharp  increases  in 
resistivity.  Wet  sieving  a  five  inch  section  of  this  core  through  a 
40  mesh  sieve  (the  section  with  the  highest  electrical  resistivity) 
produced  over  320  grams  of  shells,  coral  and  other  calcium  carbonate 
material. 

Formation  Factor-porosity  data  obtained  by  using  both  the 
simple  and  redesigned  equipment  for  the  four  marine  sediments  and 
the  lines  of  best  fit  are  presented  in  Figure  33.  Since  some  of 
this  data  was  obtained  using  slurries,  the  vertical  dashed  lines  in 
Figure  33  indicate  the  porosity  below  which  the  sediments  exist  in 
a  stable  condition. 

Discussion  of  the  factors  which  influence  the  Formation 
Factor-porosity  relationship  for  marine  sediments  is  difficult  be- 
cause with  the  exception  of  Narragansett  Bay  sand  (Station  C),  the 
marine  sediments  tested  were  mixtures.  The  mixtures  of  particles 
found  in  a  marine  sediment  can  be  seen  in  Figure  34.  This  photo- 
graph (after  a  575X  electron  microscope  enlargement)  of  the  sediment 
obtained  at  Narragansett  Bay  (Station  B)  shows  a  large  angular  sand 
grain  and  a  variety  of  plate,  rod  and  angular  shaped  particles.  The 
purpose  of  this  picture  is  not  to  determine  particle  shape  which  has 
been  found  to  have  little  effect  on  electrical  resistivity  measure- 
ments of  sediments,  but  to  illustrate  the  wide  variety  of  particle 
sizes  and  shapes  found  in  a  marine  sediment.  Measures  such  as  ^_n 
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Figure  33.      Formation  Factor-Porosity   Relationships  of  Four 
Marine  Sediments, 
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Figure  34.   Photograph  of  a  Marino  Sediment  from  Station  B 
Enlarged  575  Times  by  Electron  Microscope. 
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and  6  .  can  be  obtained  for  marine  sediments  but  the  values  for 
0 

different  sections  of  a  marine  core  sample  vary  so  widely  that 
little  can  be  concluded  from  these  measurements.   Consequently 
sediment  classification,  that  is  percentage  of  sand,  silt,  and  clay, 
is  more  useful  than  particle  size  and  shape  when  observing  the  ef- 
fect marine  sediment  particles  have  on  electrical  resistivity  meas- 
urements. For  example  it  was  observed  that  the  Formation  Factor  of 
sediments  with  a  higher  sand  content  change  less  with  a  change  in 
porosity  than  of  sediments  with  a  lower  sand  content.   The  reduc- 
tion in  negative  slope  of  the  Formation  Factor-porosity  curve  for 
Narragansett  Bay  sand  (Station  C)  can  be  seen  in  Figure  18,   If  the 
line  of  best  fit  is  extrapolated  to  the  vertical  axis  at  a  porosity 
of  100  percent,  those  sediments  with  a  higher  sand  content  resulted 
in  higher  intercept  values.  Table  11  illustrates  these  points  by 
showing  the  Formation  Factor  intercept  values  and  the  negative 
slopes  along  with  the  percentages  of  sand  in  the  marine  sediments 
tested. 


TABLE  11 

FORMATION  FACTOR  INTERCEPT  VALUES,  NEGATIVE  SLOPES, 
AND  PERCENTAGE  SAND  OF  FOUR  MARINE  SEDIMENTS 


Sediment 

(Intercepts) 

FF  at  1007.  porosity 

Slopes 

Sand  7.  range 

Narragansett  Bay 

(Station  C) 
Puerto  Rican 

sediment 
Narragansett  Bay 

(Station  B) 
Narragansett  Bay 

(Station  A) 



1.47 
1.36 
1.15 

1.00 

• 

-1.20 
-1.77 
-1.75 

-2.22 

i 

97  -  100 
70  -  50 
58  -  43  . 
33  -  41 

» 
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The  marine  sediments  which  contained  a  higher  percentage  of 
clay  were  stable  at  a  higher  porosity  than  the  sediments  with  less 
clay.  This  point  is  illustrated  in  Table  12,  which  shows  the  per- 
centage of  clay  and  the  highest  porosity  at  which  the  marine  sedi- 
ments existed  in  a  stable  condition. 

TABLE  12 

PERCENTAGE  CLAY  AND  THE  HIGHEST  POROSITY  AT  WHICH 
THE  FOUR  MARINE  SEDIMENTS  WERE  STABLE 


Name 

Clay  (7.) 

Highest 
Porosity  (7.) 

Puerto  Rican  marine  sediment 
Narragansett  Bay  (Station  A) 
Narragansett  Bay  (Station  B) 
Narragansett  Bay  (Station  C) 

7-18 
7-14 
4-9 
0 

74 
74 
67 
41 

The  sediments  that  can  exist  at  a  higher  porosity  result  in  a  larger 
range  of  Formation  Factors  since  the  range  of  porosity  at  which  elec- 
trical resistivity  measurements  can  be  made  is  greater. 

Data  such  as  coefficient  of  permeability,  porosity,  void  ratio 
and  water  content  versus  time  for  various  consolidation  pressures  is 
presented  in  Appendix  C. 


Model  Equations 

The  coefficients  of  the  empirical  model  equation  (Equation  32) 
were  determined  for  each  sediment.  A  similar  method  is  used  by  well 
logging  engineers  to  determine  empirical  model  equations  for  different 
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rock  formations.   Table  13  presents  the  empirical  models  for  each 
sediment  and  the  range  of  Formation  Factors  and  porosity  determined 
in  this  investigation. 

Since  94  percent  of  the  data  points  used  to  obtain  the  line 
of  best  fit  fell  within  +  2  percent  porosity  error  lines  (the  other 
six  percent  fell  within  +  4  percent),  these  equations  can  be  consid- 
ered good  empirical  prediction  equations  for  each  individual  sedi- 

12 
ment.   Kermabon  et  al.  "  obtained  a  porosity  error  line  of 

+  2.5  percent  with  2500  data  points  used  to  determine  a  prediction 

equation  from  21  cores  taken  from  a  40  square  mile  area  in  the 

Tyrrhenian  Sea.   The  prediction  equations  developed  by 

12       9         13 
Kermabon  et  al.,   Boyce,   and  Smith  '  (presented  in  the  Review  of 

the  Literature  section  as  Equations  41,  40,  and  43  respectively), 
were  tested  with  this  data  and  predicted  the  porosity  of  Providence 
silt,  kaolinite  clay  (0.6N  NaCl),  illite  clay  and  Narragansett  Bay 
sediment  (Station  A  and  B)  within  +  5  percent  in  the 
50  to  70  percent  porosity  range.   Although  the  "Humble  Formula", 
Equation  33,  is  the  most  widely  used  prediction  equation  for  rock 
formations,  it  is  not  at  all  applicable  to  sediments* 

The  two  factors  which  influenced  Formation  Factor-porosity 
equations  the  most  in  the  saline  environment  of  these  laboratory 
experiments,  were  the  mean  size  of  particles  (if  the  mean  size  dif- 
fered by  at  least  ten  diameters)  and  the  distribution  or  spread  of 
particle  sizes.  The  data  in  Figure  18  is  grouped  into  two  distinct 
areas:   in  one  area  were  the  sands,  and  in  the  other  area  were  the 
silts  and  clays.  As  discussed  previously,  classification  into 
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TABLE  13 

THE  COEFFICIENTS  OF  THE  MODEL  EQUATION,  POROSITY  RANGE 
AND  FORMATION  FACTOR  RANGE  OF  THIRTEEN  SEDIMENTS. 


Name  of  Sediment 


Coefficients  of  the 

Model  Equation*    Porosity   Formation 


FF  =  an 


-m 


Range  (%)  Factor  Range 


Ottawa  sand 
Glacial  sand  (1-a) 
Glacial  sand  (1-b) 
Glacial  sand  (2) 
Narragansett  Bay  (Station  C) 
Kaolinite  clay  (0.6N  NaCl) 
Illite  clay 
Montmorillonite  clay 
Kaolinite  clay  (Quadrafos) 
Providence  silt 
Narragansett  Bay  (Station  A) 
Narragansett  Bay  (Station  B) 
Puerto  Rico  marine  sediment 


FF  =  1.55  n 


-0.97 


FF  =1.57  n 


-0.96 


FF  =  1.35  n 


-1.10 


FF  =  1.45  n 


-1.10 


FF  =  1.47  n 


-1.20 


FF  «=  n 


-1.97 


FF  =  1.08  n 


-1.82 


FF  =  n 


-4.11 


FF  =  1.05  n 


-1.24 


FF  =  1.11  n 


-1.75 


FF  =  n 


-2.22 


FF  =  1.15  n 


-1.75 


FF  =  1.36  n 


-1.77 


26 
30 
31 
31 
32 
51 
55 
81 
41 
41 
47 
44 
46 


43 
47 
47 
46 
40 
92 
81 
96 
94 
62 
82 
74 
74 


3.5 
3.2 
3.1 
3.4 
4.3 
1.2 
1.6 
1.4 
1.1 
2.5 
1.5 
2.0 
2.1 


5.7 
5.0 
4.7 
5.4 
5.7 
3.8 
3.3 
2.1 
3.1 
5.2 
5.0 
4.9 
5.0 


*On  double  logarithmic  paper  (m)  of  the  model  equation  is  tho 
slope  of  the  line  of  best  fit  and  (a)  is  the  intercept  of  the  equation's 
line  with  the  vertical  axis  at  -100  percent  porosity. 
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percentages  of  sand,  silt  and  clay  generally  indicate  mean  particle 
size  and  spread.   A  prediction  equation  such  as 

FF  =  1.6/n,  (50) 

would  approximate  the  porosity  of  the  sands  of  this  study  to  within 
+  5  percent  porosity.  A  prediction  equation  such  as 

FF  =  1.1/n2,  (51) 

would  best  approximate  the  porosity  of  silts  and  clays  to  within 
+  5  percent  porosity  with  the  exception  of  montmorillonite,  Mont- 
morillonite, which  has  a  high  clay  content  (78  percent)  resulted  in 
increased  slope  value  (m).  This  probably  can  be  explained  by  the 
fact  that  this  sediment  has  a  much  smaller  mean  particle  size  and 
greater  distribution  than  found  in  sediments  containing  less  than 
70  percent  clay.  A  prediction  equation  such  as: 

FF  =  1/n3,  (52) 

or 

FF  =  1/n4,  (53) 

would  approximate  montmorillonite  and  similar  clays  better.   Of 
the  four  marine  sediments  tested  in  this  study  the  porosity  of 
three  was  predicted  within  +  5  percent  by  the  equation  which  best 
approximated  clays  and  silts,  while  the  porosity  of  Narragansett 
Bay  sand  (Station  C)  was  predicted  within  +  5  percent  by  the  equa- 
tion which  best  approximated  sands.   It  appears  that  porosity  can 
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be  predicted  within  +  5  percent  by  obtaining  an  electrical  resisti- 
vity measurement  and  knowing  the  major  type  of  sediment  present: 
(a)  sand  (b)  silt  and  clay  (c)  predominantly  clay.   Table  13  sug- 
gests that  a  more  accurate  prediction  is  possible  by  determining 
the  electrical  resistivity  and  the  porosity  of  enough  samples  to 
obtain  a  Formation  Factor-porosity  curve  (at  least  ten  samples  with 
a  wide  range  of  porosity),  and  developing  an  equation  that  approxi- 
mates that  specific  sediment.   This  is  easily  accomplished  using  the 
redesigned  equipment.   This  equation  could  then  be  used  to  predict 
the  porosity  (within  +  2  percent)  of  the  area  by  measuring  only  the 
electrical  resistivity  of  the  sediment  and  its  interstitial  water. 
As  discussed  in  the  Review  of  the  Literature  section,  most 
of  the  theoretical  model  equations  were  developed  for  dispersions 
or  porous  mediums  such  as  rocks.  However  the  following  observations 
compare  the  results  of  this  study  with  the  theoretical  work  in  the 
literature. 

1,   All  data  obtained  in  this  investigation  lie  within  the 
Parallel  and  Series  Equation  (Equations  21  and  22). 

2*   The  equations  that  considered  particle  shape  such  as  Maxwell, 
Meredith,  and  Semenov  (Equations  18,  27  and  37  respectively) 
agree  with  the  Formation  Factor-porosity  relationships  of  this 
study;  however  the  results  on  particle  shape  in  this  investi- 
gation show  that  in  sediments  shape  is  not  an  influential  factor. 

3.   The  non-dimensional  conductivity  term  (K.)  is  unnecessary  in  a 
saline  environment  (above  18.5  parts  per  thousand  at  25  C)  since 
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this  term  approaches  zero.   Equations  that  consequently  have 

no  meaning  in  a  saline  environment  are:   Equations  22  and  24. 

However  in  a  fresh  water  environment  the  term  K,  /   0  and  these 

d 

equations  may  be  applicable,  especially  if  a  clay  with  a  high 
cation  exchange  capacity  such  as  sodium  montmorillonite  is 
present. 

4.   Kaolinite  clay  with  0.02M  sodium  tetra  phosphate  solution  as 

interstitial  water  followed  the  Maxwell  Equation  (Equation  18) 
for  non-conducting  spheres  though  the  particle  shape  of  this 
clay  was  platy  (see  Figure  23-a  and  23-b).   The  structure  of 
this  sediment  was  dispersed  and  therefore  the  grains  were  not 
in  direct  contact  with  each  other.  As  stated  earlier  shape  does 
not  seem  to  influence  the  resistivity  of  sediments.   These 
results  suggest  that  the  Maxwell  Equation  is  a  good  prediction 
equation  for  sediments  of  dispersed  structure  regardless  of 
their  shape.  The  dispersed  structure  of  the  kaolinite  clay  is 
close  to  one  of  the  boundary  equations  (Parallel  Equation), 
This  is  as  expected  since  the  dispersed  orientation  is  an 
extreme  case. 

5.   In  the  general  theoretical  formula  (Equation  34)  for  resistivi- 
ty of  a  porous  medium  such  as  a  sediment,  the  term  f(c)  is  a 
factor  that  depends  on  the  amount  of  clay/#silt  in  the  sediment. 
The  results  of  this  investigation  indicate  that  the  term  f(c) 
is  more  dependent  on  the  mean  particle  size  and  distribution 
or  spread  than  on  the  clay/silt  shape  or  mineral  content. 


CONCLUSIONS 

1,  Formation  Factor-porosity  relationships  were  determined  for 
three  clays  (kaolinite,  illite  and  montmorillonite),  Providence 
silt,  four  sands  and  four  marine  sediments.   Formation  Factors 
ranged  from  1,1  to  5.9  while  porosities  ranged  from 

26  to  93  percent.   The  results  show  that  the  porosity  of  these 
sediments  correlated  well  with  electrical  resistivity  measure- 
ments.  In  fact  94  percent  of  the  data  points  fell  within 
+  2  percent  of  the  line  of  best  fit,  while  the  other  six  percent 
fell  within  +  4  percent  of  the  line  of  best  fit. 

2.  Mathematical  model  equations  and  a  graphical  presentation  of 
the  data  showed  that  the  data  grouped  into  two  classes, 

a)  sands  and  b)  silts  and  clays.   In  the  mathematical  equation 
(FF  =  an  ),  the  intercept  value  (a)  was  higher  for  sands  (1.6) 
than  for  clays  (1).  The  slope  value  (m)  was  approximately 
minus  one  for  sands,  minus  two  for  silts  and  clays,  and  minus 
three  or  four  for  sediments  having  at  least  70  percent  clay 
sized  particles.  This  description  of  (a)  and  (m)  excludes 
kaolinite  clay  particles  in  a  dispersed  structure.  While  the 
intercept  value  was  similar  to  the  other  clays  (1.05),  the 
slope  (m)  was  closer  to  that  of  sands  than  clays  (-1.24). 

Kaolinite  clay  particles  in  a  dispersed  structure  followed  the 

21 
theoretical  equation  of  Maxwell  "  for  non-conducting  spheres. 
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3.  The  cation  exchange  capacity  influenced  the  electrical  conducti- 
vity measurements  of  sodium  montmorillonite  and  illite  clays 
when  the  conductivity  of  the  interstitial  water  was  1.1  to 

3.0  mho/meter.  Above  3.0  mho/meter  this  influence  was  negli- 
gible.  The  cation  exchange  capacity  of  kaolinite  was  negligi- 
ble over  the  range  of  interstitial  water  salinities  tested 
(1.1  to  5.6  mho/meter). 

4.  The  mineral  difference  between  kaolinite  and  illite  clays  did 
not  influence  the  electrical  resistivity  measurements  in  this 
study. 

5.  As  the  particle  size  of  sediments  decreased,  electrical  resis- 
tivity increased.  When  the  mean  particle  size  of  a  sediment 
was  30  times  larger  than  that  of  a  similar  sediment  at  the  same 
porosity,  the  Formation  Factors  were  considerably  different. 
However  small  differences  in  particle  size  (about  5  times) 
were  not  detectable  by  electrical  resistivity  measurements, 

6.  The  distribution  of  particle  size  was  found  to  affect  electri- 
cal resistivity  measurements;  the  Formation  Factor  was  found  to 
increase  as  the  spread  in  particle  size  increased. 

7.  The  particle  shape  of  the  sand,  silt  and  clay  sediments  tested 
did  not  appear  to  influence  the  electrical  resistivity  measure- 
ments in  this  study.  Spherical  shapes  were  compared  with  angu- 
lar shaped  particles  and  angular  shapes  were  compared  with 
platy  shaped  particles, 

8.  The  structure  of  sands  did  not  influence  the  electrical  resisti- 
vity measurements  independent  of  porosity.   Densely  packed  sand 
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was  compared  with  loosely  packed  sand. 

9.  Consolidation  pressure  (for  normally  consolidated  sediments) 
up  to  16  psig  did  not  influence  the  electrical  resistivity 
measurements  of  the  clays  and  silts  tested  independent  of 
porosity. 

10.  The  structure  of  two  kaolinite  clay  samples  did  affect  electri- 
cal resistivity  measurements,  especially  at  lower  porosities. 
The  flocculated  structure  had  a  much  higher  Formation  Factor 
than  the  dispersed  structure  at  52  percent  porosity. 

11.  The  marine  sediments  varied  so  much  in  particle  size,  shape 
and  distribution  within  core  sections  that  percentage  sand, 
silt  and  clay  is  a  more  useful  classification.  For  the  marine 
sediments  with  higher  sand  content  the  change  in  electrical  re- 
sistivity was  less  with  a  change  in  porosity  than  the  sediments 
with  less  sand, 

12.  Unusual  marine  sediment  samples  which  fell  outside  the 

+  2  percent  error  lines  of  the  Formation  Factor-porosity  line 
of  best  fit  were  found  to  contain  large  amounts  of  shell  or 
coral. 

13.  The  electrical  resistivity  measurements  made  with  the  re- 
designed laboratory  equipment  in  this  study  allowed  perme- 
ability, tortuosity,  and  void  ratio-log  consolidation  pressure 
to  be  determined.   Though  further  study  is  required  to  validate 
this  as  a  method,  electrical  resistivity  measurements  by  the  re* 
designed  equipment  could  provide  a  fast,  easy  and  inexpensive 
method  of  testing  in  the  field  of  submarine  soil  mechanics. 
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14.   The  redesigned  laboratory  equipment  which  consisted  of  a  cell, 
an  electrode  array  and  a  signal  conditioner  is  safe,  simple, 
easy  to  operate  and  economical  ($500.00).  Also  required  for 
these  tests  are  a  standard  volt  meter  and  recording  equipment 
(a  simple  strip  chart  recorder  is  sufficient). 


II 

AN  ELECTRICAL  RESISTIVITY  MEASURING  SYSTEM 

FOR  THE  DETERMINATION  OF  POROSITY  OF 

MARINE  SEDIMENTS  IN  SITU 


INTRODUCTION 

An  Electrical  Resistivity  Measuring  System  was  designed, 
fabricated,  and  tested  as  a  method  of  predicting  porosity  of  in  situ 
marine  sediments.  The  system  consists  of  electrode  arrays,  an  elec- 
tronic measuring  circuit,  and  a  FM  multiplex  telemetry  data  link. 
These  three  components  will  be  added  to  the  Deep  Ocean  Sediment 
Probe  (DOSP). 

The  DOSP  (see  Figure  35)  was  designed  and  built  by  the 

University  of  Rhode  Island  under  contract  to  the  Naval  Underwater 

30 
Systems  Center  (NUSC),  Newport-New  London  (Lewis   ).  The  original 

mission  of  the  DOSP  was  to  obtain  the  velocity  and  attenuation  of 
compressional  sound  waves  through  the  ocean  sediment.  This  was 
accomplished  by  inserting  four  probe  legs,  located  at  the  corners  of 
the  base  structure,  incrementally  into  the  sediment  to  a  distance  of 
five  feet.  In  one  of  the  legs  an  acoustic  sparker  source  provides  a 
signal  which  is  received  by  hydrophones  located  in  the  other  legs. 
In  addition,  a  thin  walled  corer,  centrally  located  equidistant  from 
the  probe  legs,  is  used  to  collect  a  sediment  sample  for  laboratory 
correlation  with  in  situ  acoustic  and  temperature  data  obtained  at 
various  probe  penetrations.  With  a  current  design  depth  of 
5,000  feet,  the  DOSP  provides  an  excellent  stable  ocean  floor  plat- 
form for  in  situ  sediment  measurement  sensors.   In  addition  to  the 

Electrical  Resistivity  Measuring  System,  other  systems  currently  be- 

34 
ing  added  to  the  DOSP  are  a  Nuclear  Gage  (Rose  and  Roney  )  for 

129  ' 


130 


!  PROBE  DRIVE 

-  .r"?--jyn 


^  l  'Ebb*. 


*— c  .  pumpVY 


v 

r 


V/ 


_J  BASEPLATE 
I  (42  IN.  DIAMETER 


•,  I  HYDROPHONE-  \.*7 


jtej**"**. '     L 


\ «» 


Figure  35.   Photograph  of  the  Deep  Ocean  Sediment  Probe 
(DOSP) . 
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sediment  density  measurements,  a  geochemical  probe,  and  a  multi- 
coring  unit  to  obtain  longer  cores. 

The  ability  of  the  DOSP  to  obtain  samples  with  minimum  dis- 
turbance for  correlation  with  acoustical  data  is  a  major  step  for- 
ward.  However,  certain  variables  such  as:   (1)  a  measure  of  dis- 
turbance while  penetrating,  (2)  temperature  changes  and  possible 
release  of  gases  in  solution  while  returning  the  sample  to  the  sur- 
face, and  (3)  sample  handling  on  deck,  enroute  to  the  laboratory, 
and  in  the  laboratory,  are  not  monitored  and  therefore  still  leave 
an  area  for  improvement.   The  Electrical  Resistivity  Measuring  Sys- 
tem will  allow  continuous  resistivity  readings  to  be  taken  while  the 
corer  barrel  is  penetrating  the  ocean  bottom  (both  inside  and  out- 
side the  corer  barrel)  and  also  while  the  sample  is  being  returned  to 

• 
the  surface.   It  will  be  possible  to  estimate  the  degree  of  disturb- 
ance during  these  operations  and  in  addition  to  predict  sediment 
porosity  and  water  content.  Resistivity  readings  could  also  be  taken 
after  sectioning  (just  before  storage)  and  later  in  the  laboratory 
prior  to  testing  analysis,  and  thus  allow  a  measure  of  disturbance 
during  these  handling  procedures. 


REVIEW  OF  LITERATURE 

Most  electrical  resistivity  field  equipment  was  developed  by 

25 
geophysicists  for  oil  or  other  mineral  prospecting.   Dakhnov  "  re- 
viewed almost  all  the  geophysical  well  logging  equipment  used  be- 
fore 1962.   A  detailed  description  of  all  the  electronics  and 
electrode  arrays  used  is  included.   When  and  where  each  electrode 
array  is  best  used  while  surveying  is  also  discussed.   Books  and 
pamphlets  of  more  recent  equipment  used  in  geophysical  well  logging 
are  available  from  Schlumberger  Ltd.,  P.  0.  Box  2175,  Houston,  Texas, 

77001.   Electrical  resistivity  measuring  equipment  and  interpreta- 

35 
tion  of  electrical  resistivity  logs  are  included.   Fisher  "  suggests 

that  the  unique  hardware  which  has  been  developed  for  the  radical 

changes  in  pressure  and  temperature  that  occur  in  oil  well  bore 

holes  could  have  useful  application  in  the  ocean  as  well. 

The  application  of  two  different  geophysical  methods  in  order 

to  determine  the  physical  properties  of  land  soils  has  been  described 

36 
by  Borowczyk  and  Krolikowski.    These  authors  combined  electrical 

resistivity  and  radioactive  isotopic  methods  in  one  probe  called  the 
"Universal  Probe".   This  probe  consisted  of  a  steel  tube  with  elec- 
trical resistivity  electrodes  located  in  the  tube  head.   The  probe 
was  introduced  into  the  soil  by  either  a  vibrating  hammer  or  drill 
rig  -  pile  driver  which  used  a  tripod  for  vertical  support.   While 
sinking  the  probe  into  the  ground  at  intervals  of  10  or  20  cm.,  elec- 
trical resistivity  measurements  were  taken.   After  these  measurements 
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radioactive  isotopic  transmitting  and  detecting  equipment  were  low- 
ered down  the  steel  tube  (inner  diameter  3.5  cm.)  and  the  intensity 
of  the  scattered  gamma  radiation  was  determined.   The  authors  con- 
cluded that  the  Universal  Probe  allowed  accurate  soil  water  contents 
and  porosities  to  be  predicted  by  nuclear  measurements,  soil  bed  and 
thickness  identification  by  electrical  resistivity  measurements,  and 
approximate  soil  strength  by  the  penetration  resistance  measurements 
made  during  probe  insertion.  These  conclusions  were  arrived  at  by 
correlating  the  measurements  made  by  the  probe  with  known  strati- 
graphy and  soil  properties  obtained  by  borings. 

The  electrical  resistivity  surveying  techniques  that  were 

37  38 

first  developed  by  Wenner  '  were  used  in  a  land  survey  by  Perret, 

With  a  simple  apparatus  consisting  of  batteries,  a  double  commutator, 
a  null  indicator  potentiometer,  and  the  necessary  cables,  reels  and 
steel  electrodes,  he  was  able  to  determine  the  depth  of  a  marine  clay 
which  underlies  the  lower  Mississippi  Valley.  At  four  check  points 
the  electrical  resistivity  data  agreed  with  boring  data  within 
3.5  percent.   By  using  electrical  resistivity  surveying  techniques, 
it  was  possible  to  use  boring  intervals  of  5,000  feet  instead  of  the 
usual  1,000  feet.  This  resulted  in  an  80  percent  reduction  in  time 
and  a  59  percent  reduction  in  cost  when  compared  with  surveying  the 
site  only  by  boring  techniques. 

The  standard  geophysical  laboratory  or  field  electrical  re- 
sistivity equipment  used  on  land  can  be  used  to  measure  the  resisti- 
vity of  marine  sediments.   The  first  resistivity  measurements  made 
on  marine  sediments  were  conducted  on  board  ship  after  cores  had  been 
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taken.   Cell  modification  was  required  to  handle  sediments  instead 
of  formation  rocks.   More  extensive  equipment  modifications  are 
necessary  in  order  to  measure  resistivity  of  marine  sediments 
in  situ.  Since  bore  holes  are  not  provided,  a  penetration  type 
probe  is  required.   Deep  penetration  is  not  usually  possible  since 
power  sufficient  to  drive  a  probe  deep  into  the  bottom  is  not  avail- 
able on  the  deep  ocean  floor.   Since  a  stable  platform  is  usually 
unavailable,  distance  of  probe  penetration  is  difficult  to  measure 
accurately.   Either  the  electronics  and  power  must  be  encased  in  a 
water  tight  pressure  vessel  or  long  lengths  of  electrical  cable  must 

be  provided, 

9 

Boyce  measured  electrical  resistivity  of  Bering  Sea  cores  on 

board  ship  and  in  the  laboratory.  The  original  plastic  core  liner 
served  as  the  cell,  A  copper  ring  was  securely  glued  on  a  circular 
jacket.  Four  silver  alloy  needles  were  soldered  to  the  copper  ring. 
The  needle  electrodes  were  geometrically  spaced  90  degrees  apart 
around  the  ring.  When  the  circular  jacket  was  clamped  to  the  core 
liner,  the  needle  electrodes  penetrated  the  core  liner  and  made  con- 
tact with  the  sediment.  Four  of  the  circular  jackets  were  vertically 
spaced  along  the  core  liner.   This  configuration  is  similar  to  four 
ring  electrodes  in  which  the  current  electrodes  are  the  outer  two 
rings  and  the  inner  rings  are  the  potential  electrodes.  The  elec- 
tronics consisted  of  a  FRL  Resistivity  meter  (Type  3320,  manufactured 
by  Scientific  Services  Laboratories,  Inc.,  Dallas,  Texas),   This 
instrument  was  basically  a  potentiometric  nulling  circuit.   The  de- 
vice was  nulled  manually  and  the  resistance  of  the  sediment  was  di- 
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rectly  read  on  the  meter.   The  resistance  reading  was  then  correc- 
ted for  temperature  and  multiplied  by  a  cell  constant  to  give  a  re- 
sistivity value  corrected  to  25  Centigrade.   Interstitial  water  was 
squeezed  out  in  a  press  and  salinity  was  measured  by  an  induction 
salinometer.   Although  sediment  resistivities  were  obtained  both  on 
board  ship  and  in  the  laboratory,  the  resistivity  of  the  interstitial 
water  obtained  in  the  laboratory  was  used  to  obtain  Formation  Factor 

for  both  the  sediment  at  sea  and  in  the  laboratory. 

26 
The  apparatus  used  by  Pautot   to  measure  electrical  resisti- 
vity of  sediment  cores  taken  from  both  the  sea  and  a  lake  environ- 
ment was  a  quadripole  electrode  arrangement  which  was  moved  along 
the  vertical  axis  of  the  core.  The  quadripole  electrodes  were  of 
the  Wenner  arrangement,  A  Tellohm  resistance  meter  measured  the  re- 
sistance directly.  The  resistance  reading  was  multiplied  by  a  geo- 
metric constant  to  obtain  the  resistivity. 

39 

Williams,   in  order  to  determine  bottom  types  (sand,  silt 

and  clay)  designed,  built  and  tested  a  device  that  measured  the 
in  situ  Formation  Factor  at  the  water  sediment  interface.  The  For- 
mation Factor  measuring  electrodes  consisted  of  a  mutually  shared 
copper  ring  electrode  placed  on  the  perimeter  of  a  non-conducting 
plastic  disc.  Two  copper  point  electrodes  were  placed  on  each  side 
in  the  center  of  the  disc.  The  point  electrodes  were  insulated 
from  each  other.   The  disc  was  placed  at  the  water-sediment  inter- 
face to  locally  isolate  the  water  and  sediment  for  separate  resis- 
tance measurements.   Since  the  upper  and  lower  electrode  geometry 
was  identical,  the  ratio  of  the  resistance  at  the  terminals  was 
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equal  to  the  Formation  Factor.   The  electronics  of  the  system 
consisted  of  a  ratio  bridge  circuit  that  was  manually  nulled  to 
directly  read  the  Formation  Factor.   These  Formation  Factor  measure- 
ments are  based  on  the  assumption  that  the  resistivity  of  the  water 
immediately  above  the  sediment  is  the  same  as  the  resistivity  of 
the  interstitial  water  in  the  adjacent  surface  sediments.   The 
temperature  of  the  sediment  and  water  are  also  assumed  to  be  the 
same. 

Since  electrical  resistivity  measurements  on  cores  are  time 

27  28 

consuming  and  expensive,  Bouma  et  al.   and  Chmelik  et  al.  ,  modi- 
fied an  electrical  well  logging  system  to  be  used  with  a  specially 
designed  in  situ  probe.  The  in  situ  probe  consisted  of  a  heavy 
stainless  steel  tube  12  feet  in  length  with  stabilizing  fins  welded 
to  the  upper  end.  The  electrodes  consisted  of  brass  rings  mounted 
to  the  nose  of  the  stainless  steel  tube.  Rings  of  PVC  insulated  the 
electrodes  from  each  other  and  from  the  stainless  steel  rod.  The 
electrode  wiring  which  was  encased  in  a  water  tight  cabling  was 
brought  up  through  the  stainless  steel  tube  to  the  well  logging 
electronics  located  on  board  the  research  vessel.  The  electrodes 
were  designed  to  operate  in  the  normal  mode  and  in  a  guard  mode. 
The  normal  mode  is  an  array  in  which  the  distance  between  the  measur- 
ing electrode  is  large  compared  to  the  distance  between  the  current 
source  and  the  closest  measuring  electrode.   The  guard  mode  is  an 
array  in  which  a  thin  center  current  electrode  is  placed  between  two 
guard  current  electrodes.  This  configuration  causes  the  current  flow 
from  the  center  electrode  to  be  limited  to  a  horizontal  plane  rather 


137 


than  spreading  out.   The  measuring  electrodes  are  usually  located 
between  the  guard  electrodes  and  the  center  current  electrode. 

The  in  situ  probe  was  lowered  on  a  steel  cable  and  allowed  to 
fast  fall  the  last  ten  feet  by  releasing  the  brake.  A  spoked  wheel 
(18  inches  in  diameter)  which  contained  magnets  in  the  wheel  hub  was 
loosely  fitted  to  the  stainless  steel  probe  tube.  As  the  probe 
penetrated  the  bottom  the  wheel  remained  at  the  water  sediment  inter- 
face.  Magnetic  proximity  switches  were  placed  at  one  foot  intervals 
on  the  probe.   Depth  of  probe  penetration  was  obtained  by  an  over- 
riding spike  signal  caused  by  the  closing  of  the  proximity  switches 
as  they  passed  through  the  magnetic  field  located  in  the  wheel  hub. 
As  the  probe  was  extracted  an  excitation  current  and  the  return  sig- 
nal were  sent  over  a  separate  cable  to  obtain  the  resistivity 
measurements. 

Electrical  resistivity  measurements  were  made  on  the  cores 
using  the  same  electrical  well  logging  system,  but  the  electrodes 
were  mounted  to  a  sliding  carriage  which  was  in  synchronous  movement 
with  a  recorder.  Cores  were  vertically  cut  in  half  and  one  section 
was  placed  in  a  rack.  The  electrodes,  consisting  of  brass  pins,  were 
lowered  into  the  sediment.  The  equipment  was  energized  and  the  elec- 
trodes raked  the  sediment  as  the  carriage  moved  down  the  length  of 
the  core  section,  A  continuous  electrical  resistivity  reading  versus 
core  length  was  recorded. 

It  is  difficult  to  interpret  the  results  obtained  with  this 
equipment  since  temperature  is  not  accounted  for  and  interstitial 
water  salinity  and  bottom  water  salinity  were  not  given.  Consequent- 
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ly  Formation  Factors  could  not  be  obtained  and  the  work  could  not 

be  compared  with  other  research  work. 

27 

New  equipment  proposed  by  Bouma  et  al.   eliminates  the  need 

of  the  wheel  depth  indexing  device.  This  is  accomplished  by  having 
a  series  of  point  electrodes  arranged  in  a  spiral  on  a  fiberglass 
barrel  so  that  there  is  a  1.1  inch  vertical  spacing  between  elec- 
trodes. After  the  probe  penetrates  the  bottom  only  sequential 
switching  would  be  required  to  obtain  an  electrical  resistivity  versus 

depth  measurement. 

12 

Kermabon  et  al.  "  developed  both  laboratory  and  in  situ  equip- 
ment. The  resistivity  cell  used  in  the  laboratory  experiments  had  an 
electrode  array  consisting  of  four  electrode  rings.  A  16  Hz  square 
wave  alternating  current  was  fed  through  the  two  outer  ring  electrodes 
and  the  potential  difference  was  read  across  the  two  inner  ring  elec- 
trodes* The  same  16  Hz  square  wave  current  was  used  in  the  in  situ 
resistivity  probe.  The  electrode  array  was  a  four  point  electrode 
system  mounted  on  a  plastic  pad  located  at  the  end  of  the  probe.  The 
outer  two  point  electrodes  were  the  current  electrodes  while  the  in- 
ner two  electrodes  measure  the  potential  difference. 

The  electrical  resistivity  probe  weighed  100  kilograms  in 
water  and  was  13  meters  in  length.  Since  the  probe  was  designed  for 
deep  sea  work,  all  the  electronics  were  encased  in  a  large  pressure 
housing.  The  electronics  consisted  of  a  battery  power  supply,  an 
electrical  resistivity  circuit,  a  depth  monitoring  circuit  and  a  mag- 
netic tape  recording  unit.  The  operating  principles  of  the  in  situ 
probe  were  as  follows: 
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1.  The  probe  was  lowered  to  10  meters  from  the  sea  floor. 

2.  The  probe  was  released  and  allowed  to  free  fall  and  make  bottom 
penetration. 

3.  Upon  penetration  the  reference  platform  remained  on  the  surface 
of  the  sea  floor  while  the  probe  continued  its  downward  thrust. 

4.  After  the  completed  penetration  the  probe  was  extracted;  elec- 
trical resistivity  measurements  were  made  with  the  four  point 
electrode  array  and  the  penetration  was  measured  via  a  mechan- 
ical pulley  hookup  to  the  reference  platform. 

5.  The  electrical  resistivity  versus  penetration  data  was  stored 

on  tape  for  later  playback  on  a  X-Y  recorder  in  the  laboratory. 

12 
To  obtain  Formation  Factors,  Kermabon  et  al,  "  assumed  the  electri- 
cal resistivity  of  interstitial  water  was  the  same  as  bottom  water 
resistivity.   Bottom  water  resistivity  values  were  obtained  by 
measuring  bottom  water  salinity  and  temperature.   It  also  had  to  be 
assumed  that  the  temperature  of  the  bottom  water  was  the  same  as 
that  of  the  sediment  to  a  depth  of  13  meters.   This  equipment  has 
been  transferred  from  SACLANT  ASW  Research  Centre  Laboratory  at 
La  Specie,  Italy,  to  the  Marine  Science  Laboratories  University 

College  of  North  Wales  for  further  testing. 

13 
Smith  '  reported  on  resistivity  measurements  made  in  the 

laboratory  and  in  situ  by  using  commercially  available  electrical 
prospecting  equipment  and  the  Schlumberger  4-electrode  arrangement. 
To  obtain  Formation  Factors  from  the  in  situ  measurements,  the  re- 
sistivity of  the  bottom  water  was  measured  before  penetrating  the 
sediment,  then  this  resistivity  value  was  assumed  to  be  the  same  as 
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the  interstitial  water.   In  the  laboratory  the  resistivity  of  the 
interstitial  water  was  calculated  from  knowledge  of  its  salinity. 

If  porosity  of  marine  sediments  is  to  be  accurately  predicted 
by  an  Electrical  Resistivity  Measuring  System,  a  convenient  way  of 
measuring  the  resistivity  of  the  interstitial  water  must  be  found. 
One  method  involves  the  use  of  a  consolidometer  to  squeeze  out  the 
interstitial  water  in  order  to  determine  salinity  which  is  then 
easily  converted  to  resistivity.   The  squeezing  process  is  laborious, 
time-consuming,  and  requires  large  sediment  samples.   Since  samples 
must  be  analyzed  in  the  laboratory  at  a  later  date,  desiccation  may 
occur  and  introduce  error.  Also  it  is  possible  that  squeezing  causes 
error  by  filtering  the  water,  particularly  if  clay  is  present. 

A  simpler  method  is  to  measure  the  bottom  water  salinity  and 
assume  that  this  is  the  same  as  the  salinity  of  the  interstitial 

water.  The  accuracy  of  each  of  these  methods  has  been  considered. 

9 
Boyce  squeezed  the  interstitial  water  from  cores  and  found  that 

this . salinity  was  not  the  same  as  the  salinity  of  the  bottom  water. 
He  found  that  salinity  of  the  interstitial  water  varied  with  loca- 
tion and  depth.  However  since  the  measurements  were  made  three  to 
five  months  after  the  cores  were  taken,  Boyce  stated  some  desicca- 
tion probably  occurred.  The  maximum  difference  between  the  inter- 
stitial water  salinity  from  that  of  overlying  bottom  water  was 
2.3  parts  per  thousand  while  the  typical  difference  was  1,2  parts 

per  thousand, 

40 
Siever  et  al,   squeezed  the  water  from  twenty-two  cores  . 

taken  in  six  different  areas  and  concluded  that  the  variations  within 
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a  core  can  be  significant.   However  his  data  indicates  that  In  most 
cases  the  average  value  of  interstitial  water  salinity  for  a  core 

was  very  close  to  the  bottom  water  salinity  for  the  area, 

41 
Smirnov  '  studied  a  large  amount  of  empirical  data  which  was 

collected  in  the  1950 's  and  1960's  (including  the  data  of 

40 
Siever  et  al.   )  and  concluded  that  the  salinity  of  the  interstitial 

water  of  marine  sediments  does  correspond  to  the  salinity  of  the 

bottom  sea  water.   He  had  previously  proved  this  point  theoretical- 

42 
ly.   Data  by  Kullenberg  '  agrees  that  the  difference  in  salinity  be- 
tween the  two  is  small.  After  studying  Pacific  red  clay,  Atlantic 
red  clay,  calcareous  ooze,  claying  mud,  varved  clay,  and  lacustrin 
mud,  the  largest  difference  between  the  salinity  of  the  interstitial 

water  and  that  of  bottom  water  was  0.8  parts  per  thousand.  The  data 

12 

of  Kermabon  et  al.  '  also  suggests  the  error  introduced  by  assuming 

these  salinities  to  be  the  same  is  very  small;  for  2500  resistivity 
measurements  with  depths  down  to  13  meters  the  Formation  Factor- 
porosity  line  of  best  fit  had  only  a  2.5  percent  error. 

In  order  to  predict  the  porosity  of  marine  sediments  accurate- 
ly by  an  Electrical  Resistivity  Measuring  System,  it  is  important  to 
use  an  electrode  array  that  restricts  the  volume  of  the  sediment 
sampled.   Past  work  with  in  situ  measurements  has  been  done  with 

electrode  arrays  that  encompass  a  large  volume  of  sediment. 

27 

Bouma  et  al.   found  that  less  detail  was  recorded  with  the  in  situ 

probe  than  on  deck  measurements  with  laboratory  equipment.   The 
reduced  detail  was  attributed  to  the  in  situ  probe's  electrode  array 
which  sampled  a  larger  area  than  desired.  Consequently  equipment 
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was  planned  that  was  expected  to  reduce  the  area  sampled. 

12  37 

Kermabon  et  al.  '  used  the  four-electrode  system  by  Wenner 

in  the  nose  of  the  resistivity  probe.   A  detailed  derivation  of  the 
Wenner  electrode  configuration  is  presented  in  Appendix  A.  A  cur- 
rent was  passed  between  the  outer  two  electrodes,  and  the  differ- 
ence in  potential  between  the  inner  two  electrodes  was  measured. 
Theoretically,  if  the  medium  surrounding  the  electrodes  was  of 
uniform,  isotropic,  homogeneous  resistivity,  the  calculation  of  the 
true  resistivity  could  be  obtained  from  the  knowledge  of  the  current, 
the  potential  difference,  the  geometric  arrangement,  and  the  separa- 
tion of  the  electrodes.  Therefore  if  the  electrodes  were  equidistant 
and  the  distance  between  electrodes  was  large  in  comparison  with  elec- 
trode penetration  into  sediment,  the  true  resistivity  would  be: 

where  r_  is  the  true  resistivity  assuming  a  homogeneous,  isotropic, 
uniform  media; 
a,  the  electrode  spacing; 

&V,  the  potential  difference  between  the  voltage  electrodes; 
I,   the  current; 

K,   the  geometric  coefficient  for  the  particular  electrode 
array. 
However,  the  medium  around  the  electrode  array  is  not  usually 
homogeneous,  isotropic,  and  uniform;  therefore  the  value  computed  in 
equation  54  becomes  the  apparent  resistivity  (rA)»  The  apparent  re- 
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sistivity  may  be  larger  than,  less  than,  or  equal  to  the  true  resis- 
tivity.  The  difference  between  the  apparent  resistivity  and  the 
true  resistivity  will  depend  greatly  on  the  volume  of  sediment 
encompassed  by  a  specific  electrode  array  while  the  resistivity 

measurement  was  being  obtained, 

29 

Using  the  skin  depth  conduction  theory,  Ball   showed  the 

volume  sampled  by  a  Wenner  electrode  configuration  was  directly  pro- 
portional to  the  distance  between  current  electrodes  and  inversely 
proportional  to  the  frequency  of  the  current,  and  the  conductivity 
of  the  medium.   Then,  using  a  frequency  of  1000  Hz  and  a  current 
electrode  spacing  of  3  cm,,  Ball  determined  the  volume  of  material 
through  which  a  1000  Hz  A,C,  current  passes  as  a  function  of  the 
conductivity  of  the  material.  Higher  frequencies  are  undesirable 
due  to  A.C.  impedance  effects.  Using  this  electrode  configuration 
and  a  normal  range  of  marine  sediment  conductivities 
(0.015  to  0.035  mho/cm.)  the  volume  of  material  that  the  1000  Hz 
A.C.  current  could  encompass  would  range  from  0.6  to  1.2  cubic  meters. 
However  a  volume  much  smaller  than  0.6  cubic  meters  would  be  desir- 
able if  accurate  porosity  predictions  are  to  be  made.   It  is  pos- 
sible to  limit  the  volume  by  insulating  it  from  its  surrounding 
environment  (inside  the  corer  barrel)  or  by  focusing  an  electrode 
array  (outside  the  corer  barrel)  in  such  a  way  as  to  direct  the  cur- 
rent patterns  to  a  confined  volume. 


IN  SITU  MEASURING  SYSTEM 

The  Electrical  Resistivity  Measuring  System  was  built  in  the 
Ocean/Civil  Engineering  Laboratory  of  the  University  of  Rhode  Island 
for  the  DOSP,   It  consists  of  three  components:  electrode  arrays, 
an  electronic  measuring  circuit  and  a  FM  multiplex  telemetry  data 
link. 

Electrode  Arrays 

The  Electrical  Resistivity  Measuring  System  for  the  DOSP  has 
two  types  of  electrode  arrays:   the  electrode  array  used  inside  the 
corer  tube  and  the  electrode  array  used  outside  the  corer  tube.  The 
electrode  array  used  inside  the  corer  tube  (see  Figure  36)  consists 
of  four  brass  ring  electrodes  vertically  spaced  along  the  inner  wall 
of  the  corer  tube.  The  corer  tube  of  the  DOSP  is  a  continuous  five 
foot  section  of  transparent  acrylic  tubing  (O.D.  3  in.,  I.D.  2  3/4  in,). 
An  additional  14  inches  of  acrylic  tubing  (O.D.  3  1/4  in.,  I.D.  3  in.) 
is  fused  to  the  lower  end  of  the  corer  tube.  This  allows  the  adapta- 
tion of  a  corer  cutter-catcher  assembly  and  provides  additional  struc- 
tural strength  in  the  vicinity  of  the  inner  electrode  array.  The 
placement  of  the  inner  corer  tube  electrodes  was  difficult  because, 
the  inner  corer  tube  electrodes  had  to  be  placed  flush  against  the 
inner  corer  wall.   Electrode  extension  Into  the  sediment  causes 
sample  disturbance  when  the  sediment  sample  is  taken  while  recessing 
the  electrodes  reduces  the  electrical  contact  with  the  medium  sampled. 
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Figure  36.  Nose  of  the  DOSP  Corer  Tube  with  the  Inner  Corer 
Ring  Electrode  Array  and  the  Corer  Cutter-Catcher 
Assembly. 
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In  addition  the  inner  ring  electrodes  require  extra  acrylic  tubing 
material  for  structural  support.   The  apparent  resistivity  (r  )  of 
a  marine  sediment  sample  through  the  corer  tube  is  calculated  by 
the  following  formula: 

rA  "  ^aT^  "  K  I   '  {55) 

where : 

d  is  the  inner  diameter  of  the  corer  tube; 

a,   the  separation  between  the  inner  two  potential  electrode 

rings; 
AV,  the  potential  difference  between  the  voltage  electrodes; 
I,   the  current, 
K,   the  geometric  coefficient  for  the  particular  electrode 

A2 

i  i  *.       Ttd 

array  is  equal  to  -r —  • 

Equipotential  surfaces  for  a  vertical  cross  section  of  the  inner 
corer  tube  electrode  array  are  shown  for  a  homogeneous  medium  in 
Figure  37.  From  the  equipotential  surfaces  an  estimate  of  the  volume 
contributing  to  80  percent  of  the  total  resistance  involved  in  a 
resistivity  measurement  is  580  cubic  centimeters. 

The  electrode  array  used  outside  the  corer  tube  is  based  on 

the  principle  used  in  the  oil  well  logging  device  called  the 

A3 
Microlateralog.  The  Microlateralog  was  introduced  by  Doll   of  the 

Schlumberger  Well  Surveying  Corporation,  The  device  was  designed  to 

indicate  the  resistivity  of  a  small  volume  of  material  immediately 

in  front  of  the  electrode  array.  The  Microlateralog  (see  Figure  38) 
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Figure  37.     Vertical  Cross  Section  of  the   Inner  Corer  Ring 
Electrode  Array  and  Equipotential  Lines. 
(Shaded  area  is  approximately  80  percent  of  the 
total  ohmic  drop  of  potential) . 


148 


Core  .Tube 


Figure  38.   Microlateralog  Pad  Electrode  Array  and  the  80 
a,b   Percent  Equipotential  Line.   (Shaded  area  is 
approximately  80  percent  of  the  total  ohmic 
drop  of  potential) . 
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consists  of  a  center  current  electrode  and  three  circular  ring 
electrodes.   The  largest  or  outer  ring  is  the  other  current  elec- 
trode and  the  inner  two  rings  are  the  voltage  electrodes.   The 
electrodes  are  fabricated  out  of  non-corrosive  monel  and  embedded 
in  an  insulating  supporting  material.   The  small  pad  can  be  easily 
mounted  either  outside  the  corer  barrel  just  above  the  corer 
cutter-catcher  assembly,  or  on  one  of  the  DOSP  probe  legs.   The 
vertical  cross  section  showing  the  direction  and  shape  of  the  cur- 
rent lines  for  a  homogeneous  medium  are  also  shown  in  Figure  38. 
The  dark  equipotential  surface  represented  in  Figure  38  is  such 
that  80  percent  of  the  potential  drop  in  the  current  pattern  occurs 
between  the  electrode  pad  and  this  surface.  Consequently,  the 
volume  of  material  confined  within  the  pattern  extends  to  a  distance 
of  approximately  3.5  inches  from  the  electrode  pad  and  contributes 
80  percent  of  the  total  resistance  involved  in  the  resistivity 
measurement.  This  volume  for  a  homogeneous  medium  is  approximately 
380  cubic  centimeters. 

Since  both  the  inner  and  outer  corer  barrel  electrode  arrays 
require  wire  connection  along  the  corer  tube  to  the  measuring  instru- 
mentation, a  flat,  9  conductor,  polyester  insulated,  transmission 
cable,  manufactured  by  the  Burndy  Tape  Cable  Div. ,  15  Linden  Park, 
Rochester,  N.  Y. ,  is  used  to  transmit  the  current  from  the  signal 
conditioner  to  the  current  electrodes  and  to  send  the  potential  dif- 
ference from  the  voltage  electrodes  back  to  the  signal  conditioner 
(see  Figure  39).  Four  of  the  conductors  are  signal  cables  and  five 
are  ground  cables;  each  signal  cable  is  separated  by  a  ground  cable. 


/ 


150 


■v     ~^/ 


I  I 


ft 
\ 

.1 


V. 

J 


- 


"  . 


I 


" 


*»'  .H.1 jlM^^T-^ 


■     ......     »     .— ■■ 

r .    3 


11       ■!!      IIW^t« 


}         I       «fcTfc|         \ 


" 


hi   .. 


Figure  39.   Photograph  of  the  Inner  Corer  Ring  Electrode  Array 
with  the  Flat  Transmission  Cable  Hook  Up. 
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The  flat  transmission  cable  thickness  is  0.0085  inches.   The  cable 
was  fastened  to  the  corer  tube  with  three  successive  thin  filmed 
coatings  of  epqxy  cement.   With  this  configuration  the  connecting 
wiring  to  the  electrode  arrays  does  not  hinder  corer  penetration 
nor  will  the  wiring  break  due  to  abrasive  materials.  When  the 
Microlateralog  pad  was  placed  on  the  corer  tube,  the  same  type  of 
flat  transmission  cable  was  used  for  electrode  hookup. 

In  addition  a  miniature  guard  mode  electrode  array,  described 
in  the  Literature  Review  section,  was  fabricated  for  laboratory 
testing.   Seven  thin  platinum  wire  electrode  rings  were  vertically 
spaced  on  a  plastic  probe  five  inches  long  and  1/4  inch  in  diameter. 
Thin  wire  connecting  the  platinum  ring  electrodes  to  the  electronic 
circuit  was  embedded  in  the  plastic  probe  head.   The  miniature  guard 
mode  electrode  array  probe  was  not  tested  at  sea  since  the  plastic 
material  used  in  fabrication  was  not  considered  sturdy  enough  to 
withstand  sea  tests. 

A  photograph  of  a  Microlateralog  pad  and  the  plastic  probe 
with  a  guard  mode  electrode  array  is  shown  in  Figure  40.   A  photo- 
graph of  the  inner  corer  ring  electrode  array  is  shown  in  Figure  41. 
Any  of  the  three  electrode  arrays  can  be  used  with  the  electronic 
circuit. 

The  acrylic  corer  tube  with  its  inner  brass  electrode  rings 
was  fabricated  by  Plastic  Fabricators  Inc.,  East  Greenwich,  R,  I. 
The  Microlateralog  pad  and  guard  mode  probe  used  to  measure  resisti- 
vity outside  the  corer  barrel  were  fabricated  at  the  Naval  Underwater 
Systems  Center,  New  London,  Conn, 
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Figure  40.   Photograph  of  the  Microlateralog  Pad  Electrode 
Array  and  the  Miniature  Guard. 
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Figure  41.  Photograph  of  the  Inner  Corer  Ring 
Electrode  Array. 
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Electronic  Measuring  Component 

The  Electronic  Measuring  Component  consists  of  an  electrical 
conductivity  signal  conditioner,  DC/AC  inverter,  a  constant  voltage 
transformer  and  a  thermistor  assembly.  All  of  the  above  components 
are  located  in  a  pressure  case  (see  Figure  42),  with  the  exception 
of  the  thermistor  assembly  which  is  located  in  a  separate  housing 
that  is  threaded  into  the  upper  end-cap  of  the  pressure  case. 

The  pressure  case  with  the  electronics  will  be  mounted  on  the 
DOSP.   The  pressure  case  is  cylindrical  with  an  inside  length  of 
24  inches  and  inside  diameter  of  6  inches.   The  top  end-cap  is  tapped 
to  allow  the  threading  of  the  electrical  connectors  (manufactured  by 
Vector  Cables  Corporation,  5615  Lawndale,  Houston,  Texas)  and  the 
thermistor  housing.   Through  these  electrical  connectors  the  elec- 
tronic instrumentation  receives  12  volts,  1  amp.  D.C,  power  from  the 
DOSP  batteries.   In  addition,  signals  are  transmitted  and  received 
from  electrode  arrays  through  these  connectors,  and  a  D.C.  data  sig- 
nal is  sent  to  the  FM  telemetry  system  for  transmission  to  the  sur- 
face.  Both  pressure  case  end-caps  are  removable  and  the  electronic 
components  are  mounted  on  trays  which  are  connected  to  the  top  end- 
cap.  The  pressure  case  wall  thickness  is  1/2  inch  for  the  cylinder 
and  1-1/2  inches  for  the  end-caps.   Since  fabrication  is  out  of 
6061-T6  aluminum,  a  thickness  diameter  ratio  of  0.083  will  permit  a 
hydrostatic  pressure  of  5000  psig  with  a  safety  factor  of  1.30.   The 
end-caps  have  the  same  pressure  capabilities.   The  pressure  case  is 
tested  repeatedly  to  3000  psig  in  the  University  of  Rhode  Island 
hydrostatic  test  tank  before  actual  sea  tests.   The  procedure  for 
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Figure  42.   Pressure  Case  and  Electronic  Measuring  Components: 

Signal  Conditioner,  DC/AC  Inverter,  Constant  Voltage 
Transformer  and  a  Thermistor  Assembly. 
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safety  testing  pressure  cases  on  the  DOSP  is  more  thoroughly 
discussed  in  Lewis. 

The  conductivity  signal  conditioner,  manufactured  by 
Honeywell  Inc.,  Fort  Washington,  Pa.,  (Model  description  number 
552022-1002-10Q-000),  is  a  solid  state  null  balance  A.C.  amplifier. 
It  supplies  the  voltage  to  drive  the  sensor  array,  and  compares  the 
resultant  output  of  the  sensor  voltage  electrodes  with  a  stable  ref- 
erence  voltage,   A  differential  voltage  is  used  to  regulate  the 
drive  current.   The  signal  conditioner  produces  an  output  linearly 
proportional  to  the  conductivity  of  the  measured  medium.   The  out- 
put signal  range  is  0  to  5  volts  D.C.  and  is  sent  to  the  FM  tele- 
metry system. 

The  signal  conditioner  operates  as  a  potentiometer.  An  A.C. 
current  I  is  applied  through  two  sensor  current  electrodes,  CI  and 
C2  (see  Figure  A3).  The  resulting  voltage  drop,  V  ,  across  two 
measuring  electrodes,  Ml  and  M2,  is  compared  with  a  fixed  A.C.  ref- 
erence voltage,  VD.  When  they  differ,  the  resulting  error  voltage, 
K 

V  ,  alters  the  sensor  drive  current  I  •   The  amplifier  system  then 
c  c 

changes  its  output  to  maintain  V  equal  to  V^.   Since  V  is  held 
.  OR         o 

constant,  the  drive  current  (I  )  is  (by  Ohm's  Law)  inversely  pro- 
portional  to  the  medium's  resistance  and  directly  proportional  to 
its  reciprocal  conductivity. 

Because  of  the  large  measuring  error  which  would  be  intro- 
duced by  polarization  of  the  sensor  electrodes,  A.C.  rather  than  D.C. 
current  is  employed.   The  voltage  required  for  the  signal  conditioner 
is  from  107  to  127  volts  A.C.  at  50  or  60  Hz.   This  is  the  same 
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Figure  43.  Skematic  Diagram  of  Electronic  Measuring  System 
and  the  Inner  Corer  Ring  Electrode  Array. 
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frequency  supplied  to  the  sensor  electrodes  and  is  high  enough  to 
prevent  polarization.   The  A,C.  power  requirement  of  the  signal 
conditioner  makes  an  inverter  necessary.   The  inverter,  a  12  volt 
D.C.  input, square  wave  60  Hz  output,  is  manufactured  by  Fork 
Standards,  Inc.,  West  Chicago,  111.  (Model  number  32A057d).   Be- 
fore the  signal  is  sent  on  to  the  signal  conditioner,  a  constant 
voltage  transformer  (produced  by  Solar  Electric  Co.,  Chicago,  111, 
catalog  number  30498),  is  used  to  smoothe  the  signal  out.   The 
inverter  and  constant  voltage  transformer  allow  the  use  of  the 
DOSP  battery  power  supply.   This  permits  the  Electrical  Resistivity 
Measuring  System  to  be  used  on  the  ocean  floor  in  water  depths  of 
5000  feet. 

The  thermistor  assembly  and  housing  (type  number  H66),  is 
manufactured  by  Fenwal  Electronics  Inc.,  Framingham,  Mass.   The 
nominal  pressure  rating  of  the  housing  is  10,000  psig.   The  ther- 
mistor automatically  corrects  the  conductivity  readings  to  a  reading 
for  25  Centigrade.   The  thermistor  is  connected  to  an  adjustable 
compensation  network  in  the  signal  conditioner  which  alters  the  ref- 
erence voltage  to  match  sensor  output  versus  input  voltage  charge 
with  the  medium's  temperature.   The  range  of  adjustment  in  the  sig- 
nal conditioner  is  from  1.4  to  2,5  percent  per  degree  Centigrade, 
For  most  marine  conditions,  2,5  percent  per  degree  Centigrade  is 
chosen  since  this  is  the  temperature  coefficient  usually  used  for 
sea  water.   For  shallow  water  tests  the  thermistor  was  embedded  into 
the  acrylic  corer  barrel  next  to  the  brass  electrodes.   The  ther- 
mistor was  wired  to  the  cable  wire  and  sealed  with  epoxy  cement. 
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FM  Multiplex  Telemetry  Component 

The  development  of  a  telemetry  component  for  the  DOSP  was 

41 
conducted  by  Barry  '  under  contract  with  the  Naval  Underwater  Systems 

Center.   The  system  was  designed  to  command  the  DOSP  from  a  remote 

location  on  a  surface  ship,  and  receive  and  transmit  data  signals 

from  the  DOSP  in  ocean  depths  to  6000  feet.  All  signals  are  sent 

and  received  via  coaxial  cable  by  means  of  a  frequency  multiplexed 

(FM)  technique. 

The  FM  telemetry  system  for  the  DOSP  can  be  divided  into  three 
distinct  spectra:   touch  tone  signals  (25  to  115  Hz),  the  slowly 
varying  channels  (370  to  15.5  kHz),  and  the  acoustic  channels 
(150  to  450  kHz).  Resistivity  is  one  of  the  slowly  varying  data 
channels  with  a  channel  center  frequency  of  1300  Hz, 

The  Electrical  Resistivity  Measuring  Systemfs  signal  condi- 
tioner provides  the  necessary  voltage  level  (0  to  5  volts)  to  drive 
the  voltage  controlled  oscillator  of  the  telemetry  system.   After 
the  signals  have  been  frequency  modulated  they  are  multiplexed  with 
acoustic  data  and  touch  tone  commands.   This  technique  allows  all 
data  to  be  transmitted  simultaneously. 


PROCEDURE 

Laboratory  Tests 

Laboratory  tests  were  conducted  after  the  Electrical  Resis- 
tivity Measuring  System  was  fabricated  in  order  to  test  its  perform- 
ance prior  to  sea  tests.   For  analysis  and  calibration  all  three 
components  (electrode  array,  electronic  circuits  and  FM  telemetry 
data  link)  were  connected  together  and  data  recorded.   The  full 
range  voltage  (0  to  5  volts  D.C.)  signal  that  is  accepted  by  the 
FM  telemetry  data  link  was  supplied  by  the  signal  conditioner.   This 
full  range  was  obtained  by  increasing  the  salinity  of  the  solution 
that  the  electrode  array  was  sampling.   The  full  scale  conductivity 
range  that  may  be  sampled  is  from  0  to  60,000  ^ mho/cm. 

The  inner  corer  ring  electrode  array  was  calibrated  by 
filling  the  corer  tube  with  standard  solutions  of  known  conductivity 
while  the  Microlateralog  pad  and  miniature  guard  mode  electrode  ar- 
rays were  immersed  in  the  standard  solutions.   The  voltage  output 
for  each  electrode  array  in  the  standard  solutions  was  read  on  the 
signal  conditioner. 

When  the  Electrical  Resistivity  Measuring  System  was  used 
independently  of  the  DOSP  for  testing  in  the  laboratory  or  in 
shallow  water,  the  FM  multiplex  telemetry  component  was  no  longer 
required.   Under  these  conditions  power  was  supplied  either  by  D.C. 
battery  or  by  standard  115  volts,  A.C.  60  Hz  that  have  passed  the 


160  ' 


161 


Inverter.   Thus  wiring  to  the  electrode  array  and  a  volt  meter  to 
read  the  signal  conditioner  output  voltage  were  the  only  items 
required  to  make  the  conductivity  measurements. 

Ottawa  sand  and  Providence  silt  were  used  in  the  laboratory 
tests  in  order  to  compare  the  results  of  measurements  obtained  with 
the  different  electrode  arrays  with  the  data  from  earlier  laboratory 
experiments.   The  procedure  for  testing  the  inner  corer  ring  elec- 
trode array  was  as  follows:   a  sodium  chloride  solution  of  known 
conductivity  was  first  poured  into  the  corer  tube  which  had  been 
capped  to  prevent  leakage.  A  measurement  of  the  conductivity  of 

this  solution  (k  )  was  then  taken  in  the  corer  tube,  A  known  weight 
c 

of  Ottawa  sand  was  then  added  slowly.   The  sand  settled  out  rapidly 

in  a  loosely  packed  configuration.   The  length  and  the  conductivity 

of  the  sand-water  matrix  (k  )  were  measured  and  recorded  for  later 

m 

porosity  calculations.   On  completion  of  this  measurement,  the  corer 
tube  was  tapped  five  times  to  densify  the  sand.   Conductivity  meas- 
urements were  again  taken  and  the  length  of  the  sand-water  matrix 
was  again  recorded.   This  procedure  was  repeated  until  the  conducti- 
vity and  length  measurements  became  constant.   The  sand  was  then  in 
a  closely  packed  configuration. 

The  Microlateralog  (pad)  and  the  miniature  guard  mode  elec- 
trode array  (probe)  was  tested  by  using  a  large  acrylic  tube  chamber 
five  feet  long  and  one  foot  in  diameter  filled  with  40  inches  of 
Providence  silt  and  10  inches  of  saline  water  (30  parts  per  thousand). 
The  pad  and  the  probe  were  placed,  one  at  a  time,  in  the  water  direct- 
ly above  the  sediment  interface  and  the  water  conductivity  (k  )  was 

c 
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measured.   In  turn  both  the  pad  and  probe  were  lowered  into  the  silt 

arid  at  six  inch  increments  the  conductivity  of  the  sediment  (k  )  was 

*  m 

measured.   The  larger  size  of  the  pad  required  greater  force  to 
penetrate  the  sediment;  consequently  it  was  only  inserted  to  a  depth 
of  one  foot.   The  probe  was  inserted  the  entire  length  of  the  silt 
column.   Laboratory  analysis  of  the  marine  sediments  was  conducted 
in  order  to  compare  the  results  with  those  of  the  in  situ  sea  tests. 
Water  content,  void  ratio,  shear  strength,  specific  gravity,  organic 
content,  pH,  grain  size  analysis,  liquid  limit  and  plastic  limit 
were  determined.   In  addition  X-ray  dif fractograms  were  obtained  for 
the  clay-silt  size  fraction  of  each  marine  sample.   Details  of  the 
above  analysis  are  presented  in  Appendix  C. 

The  interstitial  water  salinity  was  obtained  after  the  water 
content  of  the  marine  sample  was  analyzed.   The  dried  sample  was 
mixed  with  twice  as  much  distilled  water  as  was  originally  present. 
The  slurry  was  repeatedly  mixed  and  allowed  to  settle  out,  then  a 
conductivity  measurement  of  the  solution  was  taken  with  a  Beckman 
conductivity  meter.   Since  the  number  of  ions  in  solution  remained 
constant  and  conductivity  changes  are  linear  with  dilution  change, 
the  original  conductivity  of  the  interstitial  water  can  be  determined 
by  doubling  the  conductivity  reading  of  the  dilute  solution.   These 
measurements  were  compared  with  in  situ  and  laboratory  measurements 
of  bottom  water  salinity. 

Part  of  the  marine  sediment  was  diluted  into  a  slurry  using 
bottom  water  obtained  from  the  core  sites  and  then  was  re-consolidated 
to  pressures  of  0.25,  1.60,  A. 00  and  16.00  psig  in  the  redesigned 
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laboratory  equipment  described  earlier.   This  allowed  a  Formation 
Factor-porosity  relationship  to  be  obtained.   Comparisons  were  also 
made  with  both  the  in  situ  Formation  Factor  and  porosity  values 
obtained  from  laboratory  analysis  of  the  marine  sediment  from  the 
Electrical  Resistivity  Measuring  System  corer. 

Sea  Tests 

Shallow  water  sea  tests  were  conducted  in  Narragansett  Bay 
at  the  following  locations,  which  were  chosen  for  convenience  and 
for  the  bottom  type  variety  found  there: 

Station  A  -  41°  34,i'n   -   71°  25.4'w 
Station  B  -  41°  29.5*N   -   71°  24.7'w 
Station  C  -  41°  22.7'n   -   71°  30.6'w 
The  inner  corer  ring  electrode  array  was  tested  at  Station  A. 
Station  A  had  a  sandy  silt  bottom  with  a  7  to  18  percent  clay  frac- 
tion. Since  the  depth  of  the  water  was  approximately  20  feet,  the 
corer  barrel  was  attached  to  two  ten  foot  sections  of  1^  inch  dia- 
meter galvanized  pipe.  The  procedure  at  this  station  was  as  follows: 
The  research  vessel  was  moored  with  bow  and  stern  anchors  and  the 
station  location  was  confirmed  by  a  visual  fix.  The  salinity  and 
temperature  of  the  water  column  were  measured  by  a  Beckman  model 
RS5-3  in  situ  salinometer- temperature  meter.  While  the  corer  barrel 
was  lowered  slowly  to  the  bottom  the  conductivity  was  continuously 
measured.   Since  the  FM  telemetry  component  was  not  needed  in  shal- 
low water,  a  voltmeter  was  used  to  read  voltage  output  directly  from 
the  signal  conditioner.  As  the  corer  barrel  penetrated  the  bottom, 
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conductivity  was  measured  every  six  inches.   The  corer  barrel  was 
then  pulled  out  and  continuous  conductivity  readings  were  made  un- 
til the  core  was  stored  on  deck.   Conductivity  readings  were  again 
taken  in  the  laboratory. 

Station  B,  the  site  of  the  future  Department  of  Ocean  Engineer- 
ing underwater  laboratory,  had  a  silty  sand  bottom  with  a 
A  to  9  percent  clay  fraction.   The  depth  of  the  water  was  approxi- 
mately 40  feet.   Both  the  inner  corer  ring  electrode  array  and. the 
Microlateralog  pad  array  were  tested  at  this  station.   A  six  foot 
section  of  one  inch  diameter  galvanized  pipe  was  securely  mounted  to 
the  Microlateralog  pad  which  was  then  carried  down,  positioned,  and 
inserted  into  the  bottom  by  two  trained  scuba  divers.   The  ring 
electrode  array  in  the  corer  barrel  was  also  taken  down  by  the  scuba 
divers.   An  attempt  was  made  to  systematically  survey  the  area. 
Using  the  pattern  shown  in  Figure  44,  Microlateralog  penetrations 
were  made  at  radial  distances  of  25  and  50  feet  from  where  the  corer 
was  inserted.   The  numbers  circled  in  Figure  44  were  the  positions 
that  were  measured  before  the  scuba  divers'  air  supply  was  depleted. 
The  procedure  was  the  same  as  that  used  for  Station  A,  except  that 
the  scuba  divers  handled  the  corer  barrel  and  the  Microlateralog  pad. 
White  tape  was  placed  at  six  inch  intervals  on  the  corer  barrel  (see 
Figure  45)  and  every  12  inches  on  the  Microlateralog  pad  to  allow 
visual  observation  of  penetration  depth.   The  positions  had  been  pre- 
marked  with  weighted  cans.   The  divers  held  the  pad  in  the  water 
prior  to  penetration  to  enable  a  bottom  water  conductivity  measure- 
ment.  The  Microlateralog  pad  was  inserted  into  the  bottom  at 
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Figure  44.   survey  Pattern  Planned  for  Narragansett  Bay 
Station  B. 
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Figure  45.  Corer  Barrel  with  white  Tape  Marks 
Every  six  Inches. 
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one  foot  increments.   A  communication  line  notified  the  divers  (by 
a  number  of  pulls)  that  the  measurement  was  complete,  and  then  the 
probe  was  inserted  an  additional  foot.   Since  the  Microlateralog  pad 
did  not  have  a  thermistor,  it  was  assumed  that  the  temperature  of 
the  sediment  was  the  same  as  the  temperature  of  the  bottom  water. 
From  the  later  temperature  measurement,  the  conductivity  readings 
were  corrected  to  25  Centigrade. 

Station  C,  a  shallow  water  area  near  Succotash  Point, 
Jerusalem,  R.  I.,  had  a  sandy  bottom.   The  water  was  0  to  5  feet 
deep.   Both  the  inner  corer  ring  electrode  array  and  the  Micro- 
lateralog pad  array  were  tested  at  Station  C.   The  pad  was  attached 
to  the  six  foot  galvanized  pipe  used  at  Station  B.   The  procedure 
used  was  the  same  as  for  the  other  two  stations,  except  that  the 
instrumentation  was  left  on  shore,  the  corer  barrel  and  the  pad  were 
carried  out  by  a  wader,  and  measurements  were  made  at  random.   The 
sand  bottom  was  very  hard,  consequently  the  core  obtained  was  ap- 
proximately one  foot  and  the  Microlateralog  pad  could  only  be  in- 
serted to  6  and  12  inches. 


RESULTS  AND  DISCUSSION 

Laboratory  Tests 

The  results  of  the  calibration  of  the  total  system  (all  three 
components)  are  presented  in  Figure  46.   Figure  46  shows  that  the 
0  to  5  volt  D.C.  output  signal  which  was  sent  through  the  FM  tele- 
metry component  and  6000  feet  of  coaxial  cable  was  transformed  line- 
arly to  a  range  of  -1.5  to  1.5  volts  D.C.   Calibration  curves  for 
the  inner  corer  ring  electrode  array,  the  Microlateralog  pad  array 
and  miniature  guard  mode  array  are  shown  in  Figure  47.   These  curves 
indicate  that,  for  each  electrode  array,  an  individual  linear  rela- 
tionship exists  between  the  conductivity  of  the  medium  sampled  and 
the  signal  conditioner  output  voltage.   Consequently  Equation  7  can 
be  written  as: 


r    k    V 
FF  =  -*  =  -£  =  J=L  ,  (7) 

r    k    V  ■  ' 
cram 


where  V  is  the  signal  conditioner  output  voltage  when  the  electrode 

array  is  sampling  the  continuous  phase  (water), 
and   V  is  the  signal  conditioner  output  voltage  when  the  electrode 

array  is  sampling  the  mixture  (sediment). 
Therefore  if  these  two  voltages  represented  in  Equation  7  are  meas- 
ured with  temperature  compensation,  the  Formation  Factor  of  a  sedi- 
ment can  be  readily  obtained. 
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Figure  46.     Calibration  of  FM  Telemetry  Channel  Showing 

Voltage  Input  from  the  Signal  Conditioner  and 
the   Response   Recorded  Through  FM  Telemetry 
Channel. 
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Figure  47.  Calibration  Curves  for  the  Inner  Corer  Ring 

Electrode  Array,  Microlateralog  Pad  Electrode 
Array  and  the  Guard  Mode  Electrode  Array. 
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The  results  of  laboratory  testing  of  the  inner  corer  ring 
electrode  array  with  Ottawa  sand  compared  well  with  the  results 
of  the  laboratory  equipment  described  earlier  (see  Figure  48).   The 
data  from  the  inner  corer  ring  electrode  array  are  shown  to  be  with- 
in the  +  2  percent  porosity  error  lines  of  the  line  of  best  fit 
determined  earlier  for  Ottawa  sand  by  47  Formation  Factor-porosity 
tests. 

In  Figure  49  Formation  Factors  obtained  with  the  miniature 
guard  mode  electrode  array  and  the  predicted  porosity  (obtained 
from  Figure  22)  are  plotted  versus  depth  for  Providence  silt.   A 
sediment  sample  was  obtained  from  the  surface,  the  center,  and  the 
bottom  of  the  Providence  silt  column  for  porosity  determination.  As 
shown  in  Figure  49  the  value  agrees  well  with  the  predicted  value. 
Three  measurements  were  taken  with  the  Microlateralog  pad  array 
probe  to  a  depth  of  one  foot.   These  three  Formation  Factors  along 
with  the  predicted  porosity  are  also  shown  in  Figure  49  as  a  function 
of  depth,  and  agree  well  with  other  data  on  Providence  silt. 

Sea  Tests 

In  Figure  50  ten  in  situ  Formation  Factors  for  Station  A  and 
five  in  situ  Formation  Factors  for  Station  B  are  plotted  as  a  func- 
tion of  the  porosity  determined  in  the  laboratory.   The  lines  of 
best  fit  which  were  determined  from  28  measurements  for  Station  A 
and  11  measurements  for  Station  B  by  the  redesigned  laboratory  equip- 
ment are  also  shown  in  this  figure.  Two  thirds  of  the  in  situ  Forma- 
tion Factors  fell  within  +  2  percent  of  the  lines  of  best  fit  while 
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Figure  48.   Formation  Factor-Porosity  Data  Points  of  Ottawa 

Sand  Obtained  With  the  Inner  Corer  Ring  Electrode 
Array  and  the  Line  of  Best  Fit  from  Earlier  Laboratory 
Data  of  Ottawa  Sand. 
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Figure  49.      Formation  Factors  Versus   Depth  Obtained  With  the 
Miniature  Guard  Mode  Electrode  Array,   the  Corres- 
ponding Predicted  Porosity  and  the  Measured  Values 
of  Porosity. 
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Figure  50.   In-situ  Formation  Factors  Obtained  with  the  Inner 
Corer  Ring  Electrode  Array  Versus  the  Measured 
Laboratory  Porosity  Values  for  Station  A  and  B 
and  the  Lines  of  Best  Fit  for  Stations  A  and  B 
Obtained  From  the'  Redesigned  Laboratory  Equipment, 
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one  third  were  within  +  A  percent.   This  same  information  is 
presented  versus  depth  of  the  core  in  Figure  51. 

There  are  several  factors  that  could  have  caused  this  vari- 
ation.  The  variation  was  not  caused  by  assuming  the  interstitial 
water  conductivity  was  the  same  as  the  bottom  water  conductivity, 
since  as  shown  in  Figure  52  the  interstitial  water  salinity  deter- 
mined in  the  laboratory  on  the  cores  taken  compare  well  with  the 
in  situ  bottom  water  salinities.   For  example  the  interstitial 
water  salinity  of  the  second  core  taken  from  Station  A  ranged  from 
31.0  to  31.8  parts  per  thousand  while  the  values  obtained  by  the 
inner  corer  ring  electrode  array,  the  in  situ  salinometer  and  the 

laboratory  conductivity  bridge  were  31.2,  31.1,  and  31.0  respective- 

42 
ly.   These  results  agree  with  the  work  of  Kullenberg  *  and 

41  12 

Smirnov,   and  justify  the  assumption  used  by  Kermabon  et  al.  '  that 

Formation  Factors  could  be  determined  by  using  bottom  water  salinity 
in  place  of  interstitial  water  salinity. 

The  fact  that  one  third  of  the  Formation  Factors  fell  outside 
the  +  2  percent  error  lines  would  not  appear  to  be  caused  by  the 
equipment  since  the  laboratory  tests  with  this  equipment  did  not 
have  such  a  wide  variation.   The  presence  of  shells  is  another  pos- 
sible cause  of  variation.   However  sample  disturbance  before  the 
porosity  of  the  sample  was  determined  in  the  laboratory  appears  to 
be  the  most  reasonable  explanation.   Although  disturbance  was  intended 
to  be  monitored  by  continuous  conductivity  readings  during  core  re- 
trieval, this  did  not  occur  at  Stations  A  and  B  due  to  an  error  . 
introduced  by  the  temperature  thermistor.   The  thermistor  was  de- 
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Figure  51.  Measured  Porosity  and  Predicted  Porosity  Versus 
Depth  for  the  Cores  Taken  at  Stations  A  and  B 
With  the  Inner  Corer  Ring  Electrode  Array. 
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Figure   52.      Measurements   from  Three  Methods  of   Determining 
Bottom  Water  Salinities  and  the   Interstitial 
Water  Salinity  of  the  Cores  Taken  at  the  Three 
Stations. 
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signed  to  automatically  correct  the  conductivity  readings  to 
25  Centigrade,  and  the  error  was  introduced  whenever  the  temper- 
ature sensing  element  embedded  in  the  acrylic  corer  barrel  was  dif- 
ferent from  the  temperature  of  the  sediment.   This  was  observed  as 
the  cores  were  retrieved  from  Stations  A  and  B  after  the  cores  were 
removed  from  the  refrigerator  in  the  laboratory.   This  error  could 
be  eliminated  and  disturbance  adequately  monitored  in  the  following 
ways: 

1.  Measure  the  conductivity  of  the  core  sample  after  it  has 
reached  room  temperature, 

2.  Measure  the  conductivity  of  the  core  sample  in  the  refrig- 
erator (or  ship's  reefer)  after  the  sediment  has  reached 
the  refrigerator  temperature. 

3.  Bypass  the  thermistor  in  the  electronic  circuit;  measure 
the  temperature  of  the  sample  and  correct  the  conductivi- 
ty measurement  to  25  degrees  Centigrade. 

If  the  thermistor  was  bypassed  and  the  temperature  corrected 
to  25  degrees  Centigrade  by  taking  the  temperature  manually,  the 
conductivity  changes  upon  retrieval  could  be  considered  due  to 
disturbance.   If  this  were  repeated  in  the  laboratory,  any  changes 
in  the  conductivity  readings  would  be  most  likely  due  to  disturbance, 
Disturbance  indicated  by  a  settling  of  about  one  inch  was  observed 
on  the  core  from  Station  B, 

Proof  that  disturbance  of  the  samples  before  the  porosity  was 
determined  could  have  caused  the  wider  variation  in  porosity  predic- 
tion of  the  sea  tests  is  shown  in  Figure  53.  The  sediment  column 
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Figure  53.   Formation  Factor-Porosity  Relationship  at  Different 

Stages  of  Sediment  Disturbance  in  the  Inner  Corer  Ring 
Electrode  Array  and  the  Line  of  Best  Fit  of  the  Same 
Sediment  Obtained  by  Using  the  Redesigned  Laboratory 
Equipment. 
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height  decreased  2.2  cm.  and  the  porosity  decreased  from  41  percent 
to  38  percent  after  the  corer  from  Station  C  was  hit  and  dropped. 
As  seen  in  Figure  53,  Formation  Factors  that  are  not  corrected  for 
this  sample  disturbance  lie  on  the  +  4  percent  error  line.   The  line 
of  best  fit  was  determined  for  13  measurements  by  the  redesigned 
laboratory  equipment.   Five  Formation  Factor-porosity  values  were 
determined  at  different  stages  of  disturbance  after  the  core  stood 
for  24  hours  at  room  temperature.   Since  the  sediment  column  height 
was  taken  at  each  disturbance  stage,  the  porosity  was  determined  by 
using  the  diameter  of  the  corer  barrel  and  the  weight  of  the  soil. 
A  sample  calculation  is  presented  in  Appendix  C.   Figure  53  shows 
that  if  this  disturbance  is  taken  into  account  the  porosity  predic- 
tions lie  within  +  2  percent  of  the  line  of  best  fit.   This  shows 
that  the  inner  corer  ring  electrode  detected  the  disturbance  when 
the  error  due  to  temperature  differences  was  eliminated.  When  the 
porosity  determinations  were  correlated  with  resistivity  measure- 
ments taken  after  disturbance,  the  Formation  Factors  agreed  well  with 
the  line  of  best  fit. 

An  important  advantage  of  the  Electrical  Resistivity  Measur- 
ing System  corer  was  discovered  when  the  core  cutter-catcher  assembly 
did  not  work  and  the  core  was  lost  at  Station  B.   The  loss  of  the 
core  was  immediately  detected  by  the  fact  that  the  conductivity  read- 
ing changed  to  that  of  the  surrounding  sea  water.   Thus  in  deep  sea 
operations  another  core  can  be  taken  without  bringing  the  corer  to 
the  surface. 

The  inner  corer  ring  electrode  array  can  be  used  to       !ne 
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the  salinity  of  the  water  column.   At  Station  A  the  conductivity 
readings  of  the  water  column  indicated  a  salinity  of  30.5  parts  per 
thousand  while  the  in  situ  salinometer  read  30.2  parts  per  thousand. 

The  predicted  porosity  determined  by  the  Formation  Factor 
measurements  made  by  the  Microlateralog  pad  was  obtained  by  using 
the  line  of  best  fit  for  Station  B  shown  in  Figure  50.  These  pre- 
dicted porosity  values  are  shown  as  a  function  of  depth  (to  5  feet) 
for  six  positions  at  Station  B  in  Figure  54.   For  comparison,  poros- 
ity as  determined  in  the  laboratory,  of  the  two  foot  core  obtained 
from  the  Electrical  Resistivity  Measuring  System  corer  is  also  shown 
in  Figure  54.  The  porosity  values  do  not  agree  as  well  as  expected. 
For  example  at  a  depth  of  one  foot,  the  predicted  porosity  value  of 
position  one  (25  feet  from  the  core)  was  72  percent  while  the  value 
determined  on  the  core  was  66  percent.  However  the  predicted  poros- 
ity determined  by  the  Formation  Factor  measurements  made  by  the 
Microlateralog  pad  for  six  positions  at  Station  C  and  the  line  of 
best  fit  for  Station  C  (see  Figure  53)  was  in  good  agreement  with  the 
porosity  determined  in  the  laboratory.  The  average  predicted  value 
was  42  percent  for  six  inch  pad  penetration  and  40  percent  for  one  foot 
penetration  while  the  porosity  values  determined  in  the  laboratory 
ranged  from  38  to  41  percent.   Since  the  penetration  was  limited  to 

one  foot  in  the  hard  sandy  bottom  of  Station  C,  a  vibrating  or  pneu- 

36 
matic  hammer  similar  to  Borowczyk's  and  Krolikowski' s   Universal 

Probe  would  have  to  be  used  for  deeper  penetrations. 

The  Formation  Factors  determined  in  situ  and  the  laboratory 

determinations  of  porosity,  water  content,  void  ratio,  cone  penetra- 
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Figure  54.  Predicted  Porosity  Versus  Depth  for  Six  Positions 
at  Station  B  Obtained  by  Using  the  Microlateralog 
Pad  and  th«  Measured  Laboratory  Porosity. 
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meter,  organic  content,  Atterberg  limits  and  pH  data  on  the  two 
cores  from  Station  A  and  the  one  core  from  Station  B  are  shown  in 
Figures  55  and  56  respectively.   The  in  situ  Formation  Factor  varied 
inversely  with  the  porosity  and  water  contents.   These  results  agree 
with  the  findings  on  the  marine  sediments  discussed  in  the  first 
section  of  this  paper.   Although  the  laboratory  analysis  presented 
in  these  figures  help  define  the  sediment  type  at  these  stations,  no 
other  soil  property  correlations  can  be  made  from  this  data.   If  a 
larger  number  of  cores  were  taken,  it  might  be  possible  to  correlate 
resistivity  with  another  soil  property  such  as  shear  strength,  bulk 
density  and  acoustic  properties. 
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Figure   55. 


In  Situ  Formation  Factors  Porosity,  Water  Content,  Void 
Ratio,  Cone  Penetrameter,  Organic  Content,  Atterberg 
Limits,  and  pH  of  .Two  Cores  from  Station  A.   (Organic 
content,  Atterberg  Limits  and  pH  were  determined  from 
top,  bottom  and  middle  core  sections) 
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(Organic  content,  Atterberg  Limits  and  pH 
were  determined  from  top,  bottom  and  middle 
core  sections) 
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Figure   56, 


In  situ  Formation  Factor,    Porosity,   Water  Content,   Void 
Ratio,   Cone   Ponetrameter,   Organic  Content,   Atterberg 
Limits  and  pit  of  the  Core  from  Station  B.       (Organic 

content,   Atterberg  Limits  and  pH  were  determined   from 
ton.    bottom  e  k3  middle  core    sechionsl 


CONCLUSIONS 

1.     The  Electrical  Resistivity  Measuring  System  that  was  designed 
and  fabricated  was  successfully  tested  in  the  laboratory  and  in 
Narragansett  Bay.   The  system  was  tested  by  connecting  the  electrode 
array,  the  electronic  circuits,  and  the  FM  telemetry  data  link  with 
6000  feet  of  coaxial  cable  in  order  to  simulate  the  requirements  for 
deep  ocean  testing  on  the  Deep  Ocean  Sediment  Probe  (DOSP).  Three 
electrode  arrays  were  fabricated  and  tested:   the  inner  corer  ring 
electrode  array,  the  Microlateralog  pad  electrode  array  and  the  mini- 
ature guard  mode  electrode  array.  The  predicted  porosity  values  ob- 
tained with  the  system  in  the  laboratory  were  in  good  agreement 
(+  2  percent  of  the  lines  of  best  fit)  with  values  determined  earlier 
with  the  redesigned  laboratory  equipment.  The  system  was  tested  in 
the  shallow  water  of  three  locations  in  Narragansett  Bay.   Using  the 
inner  corer  ring  electrode  array  most  of  the  data  obtained  was  in 
good  agreement  with  the  data  obtained  later  in  the  laboratory  with 
the  redesigned  laboratory  equipment.  The  data  that  was  not  in  good 
agreement  (+  4  percent  of  the  line  of  best  fit)  could  have  been  due 
to  disturbance  of  the  sample  before  the  porosity  was  determined  in 
'  the  laboratory.  The  Microlateralog  pad  electrode  array  needs 
further  testing.  Though  the  porosity  prediction  using  the  line  of 
best  fit  was  in  good  agreement  with  laboratory  porosity  values  at 
Station  C,  a  wider  variation  was  found  at  Station  B. 
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2.  The  conductivity  readings  of  the  inner  corer  ring  electrode 
array  on  the  water  column  and  the  bottom  sea  water  converted  to  a 
salinity  very  similar  to  the  readings  of  an  in  situ  salinometer, 
30.5  and  30.2  parts  per  thousand  respectively.   The  salinity  of  the 
bottom  water  as  determined  by  the  inner  corer  ring  electrode  array 
was  in  good  agreement  with  the  interstitial  water  salinity  deter- 
mined in  the  laboratory,  31.2  and  31.1  parts  per  thousand  respec- 
tively. 

3.  The  Electrical  Resistivity  Measuring  System  was  found  to  be 

an  excellent  monitoring  system.   The  System  indicated  that  the  bottom 
was  reached  by  showing  a  marked  change  in  conductivity  readings.   In 
the  same  way  the  loss  of  a  core  was  immediately  detected.   This  abil- 
ity would  be  especially  valuable  for  a  deep  ocean  corer.   Tests  in 
the  laboratory  with  a  sediment  in  the  Electrical  Resistivity  Meas- 
uring System  corer  showed  that  the  resistivity  measurements  detected 
sample  disturbance. 

4.  The  Electrical  Resistivity  Measuring  System  costs  approximately 
$1,500,  not  including  the  FM  telemetry  data  link  component  which  is 
shared  by  other  systems  aboard  the  DOSP.  This  relatively  low  cost 
does  not  include  the  50  man  hours  required  to  connect  the  system,  but 
does  include  all  the  electronic  circuits,  electrode  arrays  (including 
the  corer  barrel),  pressure  housing,  thermistor  housing,  electrical 
connectors  and  wiring.   Since  the  electronics  of  the  Electrical  Resis- 
tivity Measuring  System  and  the  redesigned  laboratory  equipment  are 
the  same,  module  components  of  these  systems  can  be  interchanged* 
This  duplicity  affords  quick  back-up  spare  parts  whenever  necessary. 
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The  Electrical  Resistivity  Measuring  System  and  its  components 
could  be  adapted  to  a  probe,  coring  device  or  any  deep  sea  plat- 
form.  It  could  be  used  in  the  deep  ocean  or  in  shallow  coastal 
waters.   If  hard  sandy  bottom  is  encountered,  a  vibrating  probe 
would  be  required  for  deep  penetration. 


FURTHER  STUDIES 

1.  The  results  of  this  study  indicate  that  electrical  resisti- 
vity measurements  could  be  used  to  determine  the  cation  exchange 
capacity  of  individual  clays.  For  example  clays  with  a  high  cation 
exchange  capacity  could  be  expected  to  have  Formation  Factors  of  less 
than  one  at  low  interstitial  water  conductivities.  Clays  with  a  low 
cation  exchange  capacity  would  always  have  Formation  Factors  greater 
than  one  and  the  Formation  Factors  would  remain  constant  regardless 
of  the  interstitial  water  conductivity  (assuming  porosity  remained 
the  same).   The  Formation  Factors  of  the  individual  clays  could  be 
correlated  with  the  cation  exchange  capacity  obtained  by  standard 
methods.  Once  established,  the  electrical  resistivity  method  would 
be  easier  and  less  time  consuming  than  the  standard  method. 

2.  The  development  of  an  induction  type  Electrical  Resistivity 
Measuring  System  should  be  pursued  for  in  situ  sediment  measurements. 
This  type  of  measuring  system  would  eliminate  direct  metal  electrode 
contact  of  the  sediment  and  could  result  in  greater  accuracy. 

3.  The  combination  of  the  Resistivity,  Acoustic,  Nuclear,  Thermal 
and  Geochemical  measuring  systems  on  the  DOSP  offer  a  unique  opportu- 
nity to  further  study  and  determine  accurately  in  situ  properties 
such  as  porosity,  water  content,  sound  velocity  and  bulk  wet  density. 
In  addition  the  ability  of  the  Electrical  Resistivity  Measuring  Sys- 
tem to  monitor  such  situations  as  touch  down,  core  loss,  and  sample 
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disturbance  make  the  system  even  more  valuable  for  deep  water. 
These  side  assets  should  be  further  studied  in  future  deep  sea 
tests. 

A.     The  use  of  the  Electrical  Resistivity  Measuring  System  in 
pollution  detection  should  be  studied  further.   For  example  a 
continuously  monitored  water  pollution  detection  device  could  be 
simply  attached  to  a  buoy.   The  cell  would  consist  of  two  small 
tubes  each  having  an  inner  ring  electrode  array.  One  tube  contain- 
ing water  of  desired  qualities  would  be  sealed.  The  other  tube, 
adjacent  to  the  first,  would  be  a  flow-through  tube  sensing  the 
environmental  water.  The  instrumentation  would  work  on  the  Formation 
Factor  concept,  i.e.  if  the  environmental  water  exactly  met  standards 
the  Formation  Factor  would  be  one.  Any  large  variation  from  unity 
would  give  warning  of  a  possible  polluting  situation. 

Perhaps  a  similar  Formation  Factor  device  could  be  used  to 
measure  the  cleansing  or  infusion  of  pollutants  in  bottom  sediments. 
These  residual  pollution  effects  in  bottom  sediments  could  be  moni- 
tored continuously  and  over  long  periods  of  time. 

5.     New  electrode  arrays  similar  to  the  probe  described  by 

27 
Bouma  et  al.   could  be  tested  for  a  wide  range  of  uses.   For  example, 

these  probes  could  be  inserted  to  considerable  depths  (20  feet)  at 

various  positions  in  an  estuary  such  as  Narragansett  Bay,  and 

seasonal  variations  in  porosity  could  be  detected. 

This  same  electrode  array  placed  on  the  legs  of  an  oil  rig 

platform  or  Texas  tower  could  serve  as  a  warning  system  for  such 

phenomena  as  bottom  erosion  and  scour* 
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6.  The  ability  of  electrical  resistivity  measurements  to  predict 
other  submarine  soil  properties  in  situ  should  be  pursued.   For  ex- 
ample, undrained  shear  strength  measured  in  situ  by  vane  shear  or  a 
penetrameter  could  be  correlated  with  electrical  resistivity  measure- 
ments. Also  relationships  could  be  developed  between  electrical  re- 
sistivity measurements  and  acoustic  properties  of  soils  or  geochemical 
properties  such  as  Eh  and  pH. 

7.  The  possibility  of  attaching  an  electrical  resistivity  elec- 
trode array  to  a  telescoping  or  vibrating  deep  ocean  floor  penetra- 
ting device  (50  feet  or  more)  merits  further  study  since  properties 
such  as  porosity,  water  content  and  density  could  be  predicted  to 
these  depths.  This  knowledge  would  be  very  useful  to  ocean  engineer- 
ing projects  such  as  ocean  floor  foundations. 

8.  A  sample  could  be  monitored  by  the  Electrical  Resistivity 
Measuring  System  while  it  is  returned  to  the  surface  from  the  deep 
ocean  if  a  second  temperature  sensing  element  was  provided.   Since  a 
thermistor  is  located  in  the  tip  of  one  of  the  DOSP's  probes,  a  sec- 
ond temperature  of  the  water  column  could  be  used  to  adjust  the  con- 
ductivity measurement  to  a  correct  value,  and  the  disturbance  of  a 
sample  during  retrieval  could  be  detected.  This  method  of  monitoring 

disturbance  could  be  checked  by  the  nuclear  gauge  suggested  by 
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Rose  and  Roney   or  an  acoustic  corer. 

9.  Electrical  resistivity  measurements  could  be  used  to  determine 
sediment  grain  size.   This  could  be  accomplished  by  fabricating  a 

30  inch  long,  one  inch  diameter  plastic  tube  with  three  or  four  sets 
of  inner  ring  electrode  arrays  placed  vertically  along  the  inner  wall 
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of  the  tube. 

A  slurry  containing  a  known  amount  of  soil,  water  and  de- 
flocculating  agent  would  be  added  to  the  cell  tube  and  electrical 
resistivity  measurements  would  be  made  at  definite  time  increments. 
As  the  soil  particles  settled,  the  resistivity  of  the  lower  array 
would  increase  while  the  resistivity  of  the  upper  arrays  would  de- 
crease.  If  the  same  assumptions  used  in  a  hydrometer  analysis  are 
applied  to  this  method,  the  electrical  resistivity  measurement  could 
then  be  related  to  particle  size.   This  method  would  be  easier  than 
the  hydrometer  method  since  the  measurements  could  be  made  electri- 
cally and  a  graphical  record  could  be  made  automatically. 
10.    Electrical  resistivity  measurements  could  be  used  to  monitor 
consolidation  in  situ.  This  could  be  accomplished  by  placing  an 
electrode  array  under  an  object  or  structure  with  the  electrodes' 
current  field  in  the  direction  of  the  sediment.  Changes  in  the 
porosity  due  to  the  weight  of  the  structure  would  be  detected  by  the 
electrical  resistivity  measurements.  Since  the  cost  of  electrode 
arrays  is  low,  several  could  be  placed  to  monitor  changes  both 
continuously  and  periodically. 
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APPENDIX  A 
DETAILED  DEVELOPMENT  OF  SOME  EQUATIONS 

1.  Development  of  "Parallel"  Equation 

Assume  that  the  medium  consists  of  a  series  of  plates,  all  par- 
allel to  the  direction  of  current  flow.  Then  the  plates  may  be  thought 
of  as  individual  resistors  in  parallel.   (Some  of  the  plates  are  sheets 
of  interstitial  water.)  To  find  the  total  conductivity,  the  individual 
conductivities  are  summed,  but  weighted  by  their  volumes: 

Vm  "  koVo  +  kdVd  •  ;  ' 

where : 

Vm  is  the  volume  of  the  mixture  (sediment), 

V  is  the  volume  of  the  continuous  phase  (water), 

V^  is  the  volume  of  the  dispersed  phase  (soil). 

dividing  by  Vm  , 

^  =  kcn  +  kd  (1-n), 
where  n  is  the  porosity,  and  dividing  by  k  , 

!2»  =  ^  =  n  +  lCs  (l-n), 
c 

or  the  Formation  Factor  (FF)  is, 

FF~i  ^  =  n  +  Rg1  (l-n) . 

2.  Development  of  "Series"  Equation 

Assume  medium  consists  of  series  of  plates,  all  perpendicular  to 
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direction  of  current  flow.  To  find  the  conductivity  of  the  mixture,  the 
individual  resistivities  are  weighted  according  to  their  volumes  and 
summed: 

rmVm=  rcV  rdvd  ' 


dividing  by  V  , 
m 


r  =  r  n  +  r,  (1-n), 

m    c     a 


dividing  by  r„  , 


/ 

I 


-  -  n  +  -  <l-„), 
o       c 


or, 


FF  =  ^  =  n  +  Rd  (1-n). 

3.  Detailed  Derivation  of  the  Resistivity  Equations  for  the  Wenner  and 
the  Schlumberger  Configurations  of  Electrodes 

The  most  direct  approach  to  measurement  of  the  value  of  sediment 
aggregate  resistivity  is  by  the  four-electrode  method.  The  four-electrode 

system  of  measuring  earth  resistivity  is  quite  old,  having  been  developed 
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by  Wenner   of  the  National  Bureau  of  Standards  in  1915.  Referring  to 

Figure  A-l  a  known  current  is  passed  between  the  outer  two  electrodes, 
and  the  difference  in  potential  between  the  inner  two  electrodes  is 
measured.  The  resistivity  can  be  calculated  from  knowledge  of  the  current, 
the  potential  difference,  and  the  geometric  arrangement  of  the  electrodes. 
Assuming  a  point  current  source  in  a  homogenous  medium,  the  Laplace 
equation  may  be  written  as: 
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V  v  =  o  . 


If  there  is  no  angular  dependence  (i.e.  the  resistivity  of  the  sediment 
is  not,  in  general,  perfectly  isotropic)  the  Laplace  equation  can  be 
written  as 

(  d2V/dR2  )  +  (2/R)  (dV/dR)  =  0  , 

where  V  is  the  electric  potential  at  a  distance  R  from  the  point  source. 
The  solution  to  this  equation  is  well  known  and  is 

V  -  (  Cx  /  R  )  +  C2  . 

To  determine  the  constants  C,  and  C«  ,   assume  that  as  R  extends  to  infin- 
ity, V  tends  to  zero.  Thus,  C«  ■  0  ,  and  so 

V  -  C  /  R  .  (1) 
By  definition,  the  resistivity  of  a  medium  is 

r  =  -(A/I)  (dV/dR)  ,  (2) 

where  r  is  the  resistivity,  I  is  the  current,  and  A  is  the  area  through 
which  it  flows.  Now  from  (1), 

dV/dR  »  -C/R2  .  (3) 

If  it  is  assumed  that  a  point  is  the  source,  then  the  current  flows 

through  a  sphere  surrounding  the  source,  and  the  area  of  this  sphere  is 

2 

simply  4TTR  .  Substituting  this  value  and  equation  (3)  into  equation 

(2),  it  becomes: 


A4 


r  -  -(4TTR  /I),  (-C/R  )  =  4ttC/I  , 


or 


C  =  rI/4-rr  .' 

Substitution  of  this  value  into  (1)  gives 

V  -  rI/4jTR  .  (4) 

A  current  source  implies  a  current  sink,  and  the  potential  of  such  a 
sink  is  just  the  negative  of  (4) .  By  superposition,  the  potential  at  a 
point  p,  V(p),  is 

V(p)  =  V+(p)  +  V~(p)  , 

where  V  (p)  and  V*(p)  are  the  potentials  at  point  p  due  to  the  current 
source  and  the  current  sink,  respectively.  Thus, 

V(p)  =  (rlMxr)  (R^1  -  R~X)  ., 

where  R,  and  R«  are  the  distances  frcra  the  point  p  to  the  current  source 
and  sink,  respectively. 

In  practical  situations,  one  measures  the  difference  in  potential 
between  two  points  a  finite  distance  from  the  source  and  the  sink,  as  in 
the  figure  below. 
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source         point  1  point  2         sink 

OX  X  0 

<— RX > 

< R3 — > 

< R4 > 

IfAVis  defined  as: 

AV    =     (potential  at  point  1)  -  (potential  at  point  2)   • 

r 

Then 

AV  -  -&-  J"    -JL    +    _J^-  -     —L.    -    _JL  .      (5) 

^         4-rr  ]_    Ri  R4         R2  R3  _ 

The  choice  of  the  electrode  spacing  is  somewhat  arbitrary.  The  two  most 
common  are  the  "Wenner"  and  the  "Schlumberger"  configurations.  These  are 
defined  by: 

Wenner  Schlumberger 
R1                    a  a 

R,  2a  3a 

R„  2a  3a 

R4  a  a 

where  a  is  a  convenient  unit  of  spacing. 
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Equation  (5)  then  reduces  to 

Wonner  r    =    4-j7a(Av/I) 

Schlumberger  r    =    3]Ta(AV/I) 
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APPENDIX  B 
DETAILS  OF  EQUIPMENT  AND  CALIBRATION 

The  basic  features  of  the  laboratory  and  in  situ  equipment,  with 
descriptive  photographs,  have  been  discussed  in  the  text.  Also  calibra- 
tion procedures  were  discussed  and  in  some  cases  calibration  curves  were 
presented.  Thi3  section  will  present  some  detailed  drawings  of  thi3 
equipment,  block  diagrams  of  electronics,  some  flow  patterns  developed 
by  the  electrode  arrays  and  calibration  curves  not  previously  given  in 
the  text. 

1.  Drawings  of  Cells  and  Electrode  Arrays. 

A.  Simple  Cell  Apparatus. 

Figures  7,  8  and  9  in  the  text  are  photographs  of  the  simple  cell 
apparatus.  Figure  B-l  is  a  drawing  of  the  simple  cell  with  two  types  of 
top  electrodes  used  in  testing. 

B.  Modified  Apparatus. 

Figures  11  and  12  in  the  text  are  photographs  of  one  sediment  cell 
and  chamber  respectively.  A  detailed  drawing  of  one  cell  tube  is  given 
in  Figure  B-2.  Other  components  of  the  modified  apparatus  such  as  sup- 
port frame,  piston  rod,  piston  head  and  bottom  cap  can  be  obtained  from 
plan  no.  G-120  on  file  in  the  University  of  Rhode  Island  Engineering 
Instrument  Shop. 

C.  Redesigned  Apparatus. 

A  photograph  of  the  redesigned  apparatus  is  presented  in  Figure 
14.  A  drawing  of  the  redesigned  cell  and  pressures  vessel  is  shown  in 
Figure  B-3. 
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Figure  B-l    Simple  Apparatus  Test  Cell  with  Variable  Sample  Length 
Electrode  and  Normal  Top  Electrode. 


B2 


8i/2u 


+  .000 

-.010" 

01 


±. 


r 

1.00" 

i 


1.750  0D—^ 


16 


o    o— 


o  - 


o    o- 


o— 


NOTE*-   MILL  FLAT    FOR 
PASSAGE    OF    WIRES 


+  .00  J" 
-.002" 

.437"     BORE 


#76  DRILL 
flO   HOLES) 


T 
.7874"  (20mm) 


■\ 


.7874"  (20  mm) 
.7874"  (20mm) 

-I 

.7874"  (20  mm) 

A 
.250" 


Figure  B-2    Modified  Apparatus  Test  Cell, 
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Figure  B-3    Redesigned  Apparatus  Pressure  Vessel  and  Test  Cell. 
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D.  Inner  Corer  Ring  Electrode  Array. 

In  the  text  Figure  36  shows  a  descriptive  drawing  of  the  inner 
corer  ring  electrode  array  and  Figure  39  is  a  photograph  of  the  same 
array.  A  drawing  giving  the  dimensions  of  this  array  is  shown  in  Figure 
B-4. 

E.  Microlateralog  Electrode  Array  and  the  Miniature  Guard  Mode  Elec- 
trode Array. 

Descriptive  and  dimensional  detail  is  shown  in  the  photographs  of 
these  two  arrays  in  Figure  40.  , 

2.  Block  Diagrams  of  the  Electronic  Instrumentation. 

A.  The  electronic  equipment  used  in  the  simple  and  modified  apparatus 
can  be  seen  in  Figure  13.  Figure  B-5  gives  the  block  diagram  of  this 
electronic  equipment. 

B.  The  block  diagram  of  the  Honeywell  Signal  Conditioner  (model  no. 
552022-1002-100-000)  is  shown  in  Figure  B-6.  The  signal  conditioner  was 
the  basic  electronic  component  used  in  the  redesigned  laboratory  apparatus 
and  the  in  situ  Electrical  Resistivity  Measuring  System. 

C.  Figure  42  in  the  text  shows  the  electronic  components  of  the  Elec- 
trical Resistivity  Measuring  System.  Figure  B-7  is  a  block  diagram  of  the 
system's  electronic  components. 

3,  Flow  Patterns  of  Electrode  Arrays. 

A.  Simple  cell  electrode  array  for  the  current  flow  between  top  and 
bottom  electrodes  resulted  in  the  flow  pattern  shown  in  Figure  B-8.  If 
the  current  was  allowed  to  flow  between  the  top  and  bottom  platinum  point 
electrodes  the  flow  pattern  shown  in  Figure  B-9  developed. 
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Figure  B-4    Inner  Corer  Ring  Electrode  Array. 
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Figure  B-5    Electronic  Block  Diagram  of  the  Simple  and  Modified 
Apparatus. 
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Figure  B-6    Electronic  Logic  and  Block  Diagram  of  the  Honeywell,  Ino, 
Signal  Conditioner  Model  Number  552022-1002-100-000. 
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Figure  B-7    Electronic  Block  of  In  Situ  Measuring  System. 
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B.  The  flow  pattern  developed  in  the  modified  apparatus  was  similar 
to  Figure  B-9. 

C.  The  cell  of  the  redesigned  apparatus  developed  a  flow  pattern 
similar  to  Figure  97  in  the  text. 

D.  The  flow  pattern  for  the  inner  corer  ring  electrode  array  and  the 
microlateralog  pad  electrode  array  are  presented  in  the  text  in  Figures 
37  and  38  respectively.  The  flow  pattern  for  the  miniature  guard  mode 
electrode  array  is  shown  in  Figure  B-10. 

4.  Calibration  of  Equipment. 

A.  Using  equation  46  of  the  text  and  measuring  the  resistivity  of 
standard  solutions  with  the  Beckman  Conductivity  Meter  the  cell  constants 
for  the  various  electrode  configurations  of  the  simple  and  modified  test- 
ing cells  were  obtained.  The  various  solutions  used,  with  measured  resis- 
tivity, resistance  and  cell  constants  for  two  electrode  configurations 
are  given  in  Table  B-l.  Also  given  are  the  geometric  cell  constants 
obtained  from  equation  46  by  dividing  the  cross  sectional  area  of  the  cell 
by  the  potential  electrode  spacing. 

B.  Figure  B-ll  is  the  calibration  curve  for  the  redesigned  testing 
apparatus.  This  calibration  curve  was  obtained  by  filling  the  cell  with 
solium  chloride  solutions  of  0.1,  0.15,  0.2,  0.5  and  0.6  normality  and 
measuring  the  signal  conditioner  output  voltage. 

C.  Figure  47  in  the  text  shows  the  calibration  curve  for  the  inner 
corer  ring  electrode  array,  the  microlateralog  pad  electrode  array  and  the 
miniature  guard  mode  electrode  array. 
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APPENDIX  C 
TEST  METHODS  AND:  DATA 

1.  Test  Methods  Used  on  Sediments. 

Figure  C-l  lists  the  general  properties  of  the  sediments  deter- 
mined in  this  investigation  and  also  identifies  the  references  where  the 
procedural  method  used  is  outlined.  Data  obtained  by  these  methods  are 
given  in  Table  3  and  Figures  27,  55  and  56  of  the  text. 

Porosity  (n)  and  void  ratio  (e)  can  be  determined  from  the  water 
content  (w)  and  specific  gravity  (Gg)  of  a  saturated  sediment  by  the 
following  method: 


and, 


or, 


e 


n  = 


n 


Gsw     , 

e 

e  +  1     ' 

Gsw 

Gcw  +  1 

B 

(CI) 
(C2) 


TABLE  C-l      ' 
SEDIMENT  PROPERTY  AND  PROCEDURAL  METHOD  REFERENCE  USED  IN  TESTING. 


Property 

1.  V/ater  Content 

2.  Liquid  Limit 

3.  Plastic  Limit 

4.  Specific  Gravity 

5.  Grain  Size  and  Hydrometer 

6.  Sediment  Classification 

7.  Descriptive  Grain  Size  Measurements 

(05o  and  60) 

8.  pH        * 

9.  Organic  Content 

10.  Undrained  Shear 

11.  Mineralogy  by  X-ray  Diffraction 

Technique 


Procedural  Method  Reference 

Lambe  (46) 
Olmsted  (47) 
Lambe  (46) 
Lambe  (46) 
Lambe  (46) 
Lambe  (46) 
Shepard  (48) 

Greweling  and  Peech  (49) 
Ackroyd  (50) 
Hansbo  (51) 
Cullity  (52) 


CI 


2.  Experimental  Methods  and  Data  Obtained. 
A.  Simple  Cell  Apparatus. 

The  porosity  was  determined  in  the  simple  cell  by  measuring  the 
length  of  the  sediment  sample  in  the  test  cell  (L ),  the  diameter  of  the 
cell  (Dg),  the  weight  of  soil  in  the  test  cell  (Wg)  and  the  specific  grav- 
ity of  the  soil  (G  ).  Then  the  total  volume  of  the  sample  (VL)  is, 

VT  =?■  D2L   =  V  +v   ,  (C4) 

1      4     S  S        8     V 

where :  , ' 

V  is  the  volume  of  the  soil, 

8 


V  is  the  volvme  of  the  voids, 


and, 

V      «      s 

Then, 

V              t 

and  the 

porosity  (n)   is, 

(C5) 


-  V  (C6) 

6  , 


(C7) 

The  procedure  described  in  the  text  was  used  to  measure  the  resis- 
tivity of  the  sediment  and  interstitial  water.  Formation  Factors  were 
determined  by  using  equation  3  in  the  text.  Table  C-2  lists  Formation 
Factor  -  porosity  values  obtained  for  various  sediments  with  the  simple 
cell  apparatus. 

B.  Modified  Apparatus. 

The  porosity  was  determined  in  the  modified  cell  by  measuring  the 
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TABLE  C-2 

FORMATION  FACTOR  -  POROSITY  DATA  OBTAINED  WITH  THE  SIMPLE  CELL 
FOR:  OTTAWA  SAND,  GLACIAL  SAND  (1-a),  GLACIAL  SAND  (1-b) 
GLACIAL  SAND  (2),  NARRAGANSETT  BAY  SAND, 
PROVIDENCE  SILT,  KAOLINITE  CLAY,  AND  ILLITE  CLAY. 


OTTAWA  SAND 

Formation 

Porosity  (%) 

Factor 

40.5 

3.54 

40.0 

3.69 

40.4 

6.00 

40.9 

3.74 

34.9 

4.26 

34.0 

4.44 

35.4 

4.44 

35.4 

4.36 

41.1 

3.49 

41.6 

3.50 

42.4 

3.48 

42.1 

3.42 

34.6 

4.24 

34.2 

4.24 

35.3 

4.31 

36.8 

4.20 

40.5 

3.73 

35.7 

4.56 

32.3 

4.82 

29.9 

5.03 

27.3 

5.09 

35.7 

4.43 

31.2 

4.69 

28.7 

4.88 

26.0 

5.06 

41.2 

3.43 

35.9 

4.18 

33.8 

4.50 

33.3 

4.54 

32.8 

4.61 

40.7 

3.50 

34.0 

4.26 

41.1 

3.48 

35.0 

4.27 

41.3 

3.56 

35.8 

4.25 

41.9 

3.52 

GLACIAL  SAND  (1-a) 

Formation 


Porosity  (%) 


Factor 


46.0 

3.19 

46.7 

3.14 

46.5 

3.23 

47.6 

3.10 

39.5 

3.87 

39.9 

3.83 

40.4 

3.97 

43.9 

3.47 

39.6 

4.04 

35.7 

4.46 

31.2 

4.71 

29.9 

4.80 

39.6 

4.00 

36.7 

4.29 

33.5 

4.46 

31.2 

4.63 

40.5 

3.88 

GLACIAL  SAND   (1-b) 

Formation 

Porosity  (%) 

Factor 

46.8 
47.0 
46.9 
46.9 
44.0 
44.6 
44.6 
43.9 
43.9 
41.4 
39.6 
35.7 
32.3 
31.2 


3.14 
3.08 
3.17 
3.11 
3.21 
3.21 
3.35 
3.18 
3.43 
3.57 
3.80 
4.25 
4.63 
4.80 
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TABLE  C-2   (CONT.) 


OTTAWA  i 

5A1ID 

Formation 
Factor 

GLACL* 
Porosity 

J,  SAND  (1-b) 

Porosity  (%) 

(%) 

Formation 
Factor 

35.9 

4.13 

38.7 

3.83 

42.7 

3.55 

36.7 

4.16 

35.9 

4.32 

34.6 

4.53 

39.4 

3.76 

31.2 

4.77 

37.6 

4.05 

35.6 

4.32 

34.4 

4.38 

NARRAGANSETT  BAY 

33.8 

4.44 

STATION 

C 

41.1 
33.3 

3.49 
4.54 

Porosity 

(%) 

Formation 
Factor 

40.5 

4.21 

GLACIAL  SAND  I 

35.7 

4.97 

Porosity  (%) 

Formation 
Factor 

33.5 
32.3 
31.2 

5.50 
5.66 
5.77 

45.6 

3.27 

35.7 

4.94 

41.4 

3.50 

34.6 

5.26 

38.8 

3.72 

33.5 

5.47 

37.6 

3.98 

31.2 

5.73 

•   36.5 

4.20 

39.1 

4.53 

34.4 

4.58 

38.5 

4.55 

34.0 

4.67 

37.0 

4.73 

45.1 

3.47 

34.8 

5.02 

40.7 

3.65 

38.2 

4.08 

36.7 

4.29 

PROVID 

SILT 

35.5 
34.3 
43.9 

4.55 
4.91 
3.62 

Porosity 

(%) 

Formation 
Factor 

42.4 

3.79 

60.5 

2.58 

36.4 

4.25 

59.9 

2.86 

35.2 

4.50 

57.7 

3.08 

33.0 

4.86 

57.5 

2.89 

45.1 

3.62 

41.3 

3.84 

38.2 

4.18 

KAOLI 

CLAY 

36.6 
35.0 
34.4 

4.35 
4.48 
4.75 

Porosity 

(%) 

Formation 
Factor 

45.5 

3.56 

91.0 

1.20 

40.2 

3.93 

89.7 

1.24 

38.1 

4.23 

88.0 

1.25 

36.6 

4.51 

87.7 

1.28 

35.0  " 

4.74 

84.0 

1.40 

33.4 

5.05 

78.7 

1.53 
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TABLE  C-2  (COOT.) 


GLACIAL  SAND  (2) 


ILLITE  CLAY 


Porosity  (%) 

Formation 
Factor 

42.3 

3.74 

37.7 

4.32 

33.5 

5. CO 

28.7 

5.54 

27.3 

5.69 

34.6 

4.85 

32.3 

5.19 

29.9 

5.62 

27.3 

5.77 

37.9 

4.02 

Porosity  (%) 


Formation 
Factor 


81.0 

1.60 

80.0 

1.63 

78.7 

1.76 

..  7 
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height  of  the  sediment  column  at  the  end  of  the  test  (h  )  and  the  weight 

of  the  soil  (W  ) .  Then  knowing  the  cross-sectional  area  of  the  sample 
s 

(A)  and  the  specific  gravity  of  the  soil  (G  )  the  total  volume  (V  )  can 

s  T 

be  written  as: 


VT  «  h^  ,  (C8) 


and, 


or. 


Vv  -  VT  "  Vs  '  ,  (C6) 


vv  -  VT  -  £  *  (09) 

8 


and  dividing  by  V„, 


v 


V         w 

,  i  .  Is .  (cio) 


T         V  hf 


or, 


n,  -  1  -  ^  .  <CU) 


where, 


5L 

GsA  ' 


(C12) 


and  n,  is  the  porosity  of  the  sediment  column  of  height  h..  The  porosity 
at  other  consolidation  pressure  height  increments  can  be  obtained  by  using 
equation  Cll  and  the  corresponding  height  recorded  at  each  pressure 


C6 


increment.  Then  the  water  content  (w)  and  void  ratio  (e)  can  be  obtained 
by  the  following  equations: 

n    ,  (C13) 

1  -  n 


and. 


W  -  §-  .  (CI) 

s 


Electrical  resistivity  measurements  of  the  sediment  were  made  at 
each  pressure  increment  and  Formation  Factors  were  determined  by  assuming 
the  resistivity  of  the  interstitial  water  was  that  of  the  chamber  water. 
Table  C-3  contains  the  incremental  height  of  the  sediment  column  and  the 
corresponding  consolidation  pressure,  Formation  Factor,  porosity,  void 
ratio,  water  content,  and  C-constants  (equation  C12)  for  illite  and  kool- 
inite  clays  obtained  with  the  modified  testing  apparatus. 
C.  Redesigned  Apparatus. 

Porosity  was  determined  at  the  beginning  and  at  the  end  of  a  test 
by  measuring  the  water  content  and  specific  gravity  of  the  soil  and  using 
equation  C3.  These  porosity  values  were  then  compared  to  Formation  Factor 
values  determined  by  the  method  described  in  the  text.  Figure  29  in  the 
test  is  an  example  of  the  recorded  data  obtained  with  the  redesigned 
equipment.  If  the  total  test  time  (t^)  is  known,  along  with  the  time 
prior  to  rapid  decrease  in  Formation  Factor  values  (t, )  and  the  time  For- 
mation Factor  values  again  became  relatively  constant  (t~),  and  if  the 
corresponding  Formation  Factors  at  these  two  times  (t,  and  t„)  is  also 
known,  then  the  recorded  curves,  such  as  Figure  29,  can  be  produced. 
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TAELE  C-3 

INCREMENTAL  SEDIMENT  SAMPLE  HEIGHT  WITH  CORRESPONDING  CONSOLIDATION 
PRESSURE,  FORMATION  FACTOR,  POROSITY,  VOID  RATIO,  WATER  CONTENT,  AND 
C-CONSTANT  FOR  KAOLINITE  AND  ILLITE  CLAY  OBTAINED  WITH  THE 
MODIFIED  TESTING  APPARATUS. 


(Kaolinite  Clay) 


•HO  ^» 

iH  <U  -ho  «-»  *4  ^ 

«tfW  -M  Cl«  few  v-^+j 

G  4-»  T>    O  O  *—  *»,  O  '  (dO 

OP!  -HM  -H  -M  w  4J  -M  w 

See  o  w  «a  o  o  o  ^  o  c 

M^    U  WW  E  4-»  O  f)  -H  :Q)-M  O— > 

o  ti  *-*  co  Ro  ^  --»  -h+j  -m  e:  od 

CO)  O    M  O    tf  OC  0(0  (dO  It) 

H  CO  .3  UP<  PhUa  -               Pk«~'  >C«  ^O  O  *-" 


20.5 

^ 

1.59 

76.8 

3.31 

127.18 

4.76 

18.3 

.25 

1.83 

74.0 

2.84 

109.41 

4.76 

15.9 

.40 

1.94 

70.1 

2.34 

90.00 

4.76 

14.4 

.80 

2.13 

66.9 

2.02 

77.88 

4.76 

13.3 

1.65 

2.42 

64.2 

1.79 

69.00 

4.76 

13.0 

2.00 

2.54 

63.4 

1.73 

66.57 

4.76 

12.4 

3.00 

2.72 

61.6 

1.60 

61.72 

4.76 

20.3 

m. 

1.56 

76.4 

3.23 

124.17 

4.80 

20.1 

.20 

1.68 

76.1 

3.19 

122.60 

4.80 

17.9 

.40 

1.81 

73.2 

2.73 

104.95 

4.80 

14.8 

.80 

2.05 

67.6 

2.08 

80.14 

4.80 

14.3 

2.00 

2.14 

66.4 

1.98 

76.11 

4.80 

21.2 

- 

1.64 

77.7 

3.48 

133.93 

4.73 

20.3 

.23 

1.77 

76.7 

3.29 

126.61 

4.73 

18.0 

.42 

1.83 

73.7 

2.81 

107.89 

4.73 

14.4 

.85 

2.17 

67.2 

2.04 

78.62 

4.73 

13.2 

1.65 

2.46 

64.2 

1.79 

68.88 

4.73 

12.2 

3.30 

2.79 

61.2 

1.58 

60.74 

4.73 

11.4 

6.00 

3.02 

58.5 

1.41 

54.24 

4.73 

16.3 

— 

1.69 

74.0 

2.85 

109.75 

4.23 

14.4 

.25 

1.97 

70.6 

2.40 

92.47 

4.23 

12.8 

.40 

2.21 

67.0 

2.03 

77.92 

4.23 

11.5 

.80 

2.55 

63.2 

1.72 

66.10 

4.23 

10.8 

1.65 

2.70 

60.8 

1.55 

59.61 

4.23 

10.5 

3.00 

2.81 

59.7 

1.48 

56.92 

4.23 

16.2 

mm 

1.69 

76.5 

3.26 

125.50 

3.80 

14.1 

.25 

1.83 

73.0 

2.71 

104.25 

3.80 

12.6 

.40 

2.04 

69.8 

2.32 

89.07 

3.80 

11.1 

.80 

2.36 

65.8 

1.92 

73.89 

3.80 

10.6 

1.00 

2.39 

64.2 

1.79 

68.83 

3.80 
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TABLE  C-3  (CONT.) 

(Kaolinite  Clay) 
Pc        FF        n         e  w 


16.1 

1.63 

74.6 

2.94 

112.94 

4.09 

14.4 

.25 

1.84 

71.6 

2.52' 

96.95 

4.09 

12.5 

.40 

2.14 

67.3 

2.06 

79.09 

4.09 

.11.3 

.80 

2.46 

63.8 

1.76 

67.80 

4.09 

.  10.5 

1.70 

2.68 

61.1 

1.57 

60.28 

4.09 

16.8 

— 

1.67 

74.8 

2.96 

113.93 

4.24 

14.7 

.25 

1.92 

71.2 

2.47 

94.88 

4.24 

13.2 

.40 

2.22 

67.9 

2.11 

81.28 

4.24 

12.0 

.80 

2.55 

64.7 

1.83 

70.39 

4.24 

11.0 

1.65 

2.84 

61.5 

1.59 

61.32 

4.24 

10.5 

3.30 

2.89 

59.6 

1.48 

56.79 

4.24 

15.9 

- 

1.70 

73.4 

2.76 

106.11 

4.23 

14.5 

.25 

1.83 

70.8 

2.43 

93.38 

4.23 

13.2 

.40 

2.10 

68.0 

2.12 

81.68 

4.23 

.  11.9 

.80 

2.32 

64.5 

1.81 

69.74     . 

4.23 

11.0 

1.65 

2.59 

61.5 

1.60 

61.56 

4.23 

10.7 

2.00 

2.65 

60.5 

1.53 

58.83 

4.23 

16.2 

— 

1.65 

74.9 

2.98 

114.62 

4.07 

11.8 

.80 

2.41 

65.5 

1.90 

73.05 

4.07 

10.6 

1.65 

2.79 

61.6 

1.60 

61.71 

4.07 

20.6 

— 

1.80 

75.9 

3.15 

121.28 

4.96 

20.3 

.25 

1.80 

75.6 

3.09 

118.95 

4.96 

15.8 

.40 

2.03 

68.6 

2.19 

84.06 

4.96 

14.4 

.80 

2.29 

65.6 

1.90 

73.20 

4.96 

13.3 

1.65 

2.56 

62.7 

1.68 

64.67 

4.96 

12.5 

3.30 

2.78 

60.3 

1.52 

58.47 

4.96 

11.6 

6.60 

3.10 

57.2 

1.34 

51.49 

4.96 

(Illite  Clay) 

h 

Pc 

FF 

n 

e 

w 

c 

19.7 

. 

2.05 

70.2 

2.36 

84.59 

5.91 

18.8 

.40 

2.16 

68.6 

2.18 

78.30 

5.91 

17.7 

.80 

2.28 

66.6 

1.99 

71.47 

5.91 

16.9 

1.65 

2.38 

65.0 

1.86 

66.56 

5.91 

16.6 

2.00 

2.50 

64.0 

1.78 

63.72 

5.91 

20.1 

- 

2.12 

70.2 

2.42 

86.60 

.  5.98 

19.5      ". 

.40 

2.14 

69.4 

2.27 

81.29 

5.98 

18.3 

.85 

2.20 

67.3 

2.06 

73.77 

5.98 

17.3 

1.70 

2.28 

65.4 

1.89 

67.75 

5.98 

16.2 

3.30 

2.78 

63.1 

1.71 

61.29 

5.98 

15.9 

4.00 

2.90 

62.4 

1.66 

59.48 

5.98 

19.3 

- 

2.08 

70.8 

2.42 

86.91 

5.64 

18.6 

.40 

2.20 

69.7 

2.30 

82.45 

5.64 

17.1 

.70 

2.39 

67.0 

2.03 

72.77 

5.64 
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TABLE  C-3  (COiYf.) 

(Illite  Clay) 
Pc        FF        n         e  w 


19.6 

1.94 

70.5 

2.39 

85.66 

5.78 

18.8 

.20 

1.97 

69.3 

2.26 

80.91 

5.78 

18.2 

.45 

2.20 

68.2 

2.14 

76.87 

5.78 

17.4 

.  .85 

2.28 

66.8 

2.01 

72.11 

5.78 

16.6 

1.65 

2.40 

65.2 

1.87 

67.15 

5.78 

15.6 

3.30 

2.63 

62.9 

1.70 

60.77 

5.78 

14.7 

6.00 

2.93 

60.7 

1.54 

55.36 

5.78 

18.5 

.20 

2.09 

69.2 

2.25 

80.69 

5.69 

16.3 

.85 

2.16 

65.1 

1.86 

66.83 

5.69 

15.4 

1.65 

2.28 

63.1 

1.71 

61.16 

5.69 

14.6 

3.30 

2.44 

61.0 

1.57 

56.13 

5.69 

14.2 

5.00 

2.60 

59.9 

1.50 

53.61 

5.69 

18.2 

.20 

2.14 

68.4 

2.16 

77.41 

5.76 

17.7 

.40 

2.16 

67.4 

2.07 

74.30 

5.76 

16.9 

.85 

2.24 

65.9 

1.93 

69.32     . 

5.76 

16.1 

1.65 

2.33 

64.2 

1.80 

64.34 

5.76 

15.2 

3.30 

2.56 

62.1 

1.64 

58.74 

5.76 

14.4 

6.00 

2.80 

60.0 

1.50 

53.76 

5.76 

19.5 

.20 

1.95 

69.6 

2.29 

82.22 

5.92 

19.2 

.40 

1.98 

69.2 

2.24 

80.40 

5.92 

17.0 

.85 

2.05 

65.2 

1.87 

67.08 

5.92 

16.4 

1.65 

2.20 

63.9 

1.77 

63.45 

5.92 

16.2 

2.00 

2.27 

63.5 

1.74 

62.24 

5.92 

17.7 

.20 

2.08 

67.9 

2.11 

75.65 

5.69 

17.3 

.40 

2.16 

67.1 

2.04 

73.13 

5.69 

16.5 

.85 

2.21 

65.5 

1.90 

68.10 

5.69 

15.6 

1.65 

2.36 

63.5 

1.74 

62.42 

5.69 

14.7 

4.00 

2.64 

61.3 

1.58 

56.76 

5.69 

20.8 

.25 

2.00 

68.7 

2.15 

77.06 

6.61 

20.7 

.40 

2.13 

68.1 

2.13 

76.34 

6.61 

20.6 

.85 

2.14 

67.9 

2.12 

75.85 

6.61 

18.3 

1.65 

2.21 

63.9 

1.77 

63.39 

6.61 

17.2 

3.30 

2.39 

61.6 

1.60 

57.42 

6.61 

16.3 

6.00 

2.63 

59.5 

1.47 

52.55 

6.61 

20.8 

.20 

2.09 

68.4 

2.17 

77.62 

6.57 

20.2 

.40 

2.19 

67.5 

2.08 

74.37 

6.57 

19.9 

.80 

2.19 

67.0 

2.03 

72.71 

6.57 

19.6 

1.00 

2.40 

66.5 

1.98 

71.09 

6.57 

20.5 

.20 

1.98 

68.2 

2.15 

77.01 

6.51 

18.6 

.40 

2.09 

65.0 

1.86 

66.56 

6.51 

18.4 

1.00 

2.22 

64.6 

1.83 

65.47 

6.51 

18.0 

2.00 

2.26 

63.8 

1.76 

63.25 

6.51 
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Table  C-4  contains  the  data  obtained  with  the  redesigned  equipment,  the  total 
test  time,  (t_),  times  t.  and  t?,  along  with  the  corresponding  Formation 
Factors,  porosity,  void  ratios  and  water  contents  for  these  times.  Also 
in  this  figure  is  the  length  of  the  sediment  column  at  the  beginning  and  at 
the  end  of  a  test,  the  water  collected  in  the  reservoir  and  the  consolidation 
pressure  used  in  each  test. 

D.  Statistical  Technique  and  Methods. 

Statistical  analysis  of  data  was  accomplished  with  a  statistical 
package  of  computer  programs  called  **STATPACK  (53).  The  computer  programs 
are  available  at  remotely  located  terminals  and  prior  knowledge  of  computer 
or  programming  is  not  necessary.  The  analyses  of  STATPACK  used  in  this 
study  were: 

Transformation    to  obtain  logarithmic  functions  of  Formation 

Factor  and  porosity  data. 


Regression 


to  obtain  a  functional  relationship  among  For- 
mation Factor  (FF)  and  porosity  (n) .  The  rela- 
tionship considered  here  is  of  the  form: 


Log  FF  *=  log  a  -  m  log  n, 

where  FF  is  called  the  dependent  variable  and 
n  the  independent  variable.  The  slope  of  the 
line  is  (m)  and  the  intercept  value  is  (a). 
Line  of  best  fit  and  the  coefficients  m  and  a 
were  determined  by  the  method  of  least  squares. 
Other  statistics  such  as,  correlation  coeffic- 
ient, sura  of  squares,  attributed  to  the 
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regression,  sum  of  the  squares  of  deviations 
from  the  regression,  F-values  for  analysis  of 
variance,  standard  error  of  estimate,  standard 
error  of  regression  coefficient,  corapxited  t- 
value,  and  residuals  were  also  obtained. 

Scatter  Diagram    to  obtain  Formation  Factor  -  porosity  printed 

scatter  diagrams  of  data.  The  line  of  best  fit 

+ 
and  the  -  2  percent  porosity  error  line  were 

also  placed  on  the  scatter  diagrams. 

Tables  C-5  and  C-6  along  with  Figure  C-l  outline  the  procedure  used 
on  all  sediments  analyzed.  Table  C-5  shows  the  values  of  Formation  Factor  - 
porosity  data  obtained  for  Providence  silt.  Also  given  are  the  values  of 
the  logarithmic  transformation  of  the  Formation  Factor  and  porosity.  Table 
C-6  shows  the  intercept,  regression  coefficient,  standard  error  of  regres- 
sion coefficient,  computed  t-value,  correlation  coefficient,  standard  error 
of  estimate  and  the  analyssi  of  variance  for  the  regression  technique. 
Figure  C-l  is  a  scatter  diagram  of  the  Formation  Factor  -  porosity  data. 

The  solid  line  is  the  line  of  best  fit  and  the  two  dashed  lines  are  the 

•f- 

-  2  percent  porosity  error  lines. 
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Figure  C-2 


Scatter  Plot  of  the  Data  in  Table  C-5  with  Log  Formation 
Factor  as  the  Horizontal  Variable  and  the  Log  Porosity  as 
the  Vertical  Variable. 
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nd  equipment  for  this  purpose  was  designed,  fabricated  and  tested.   For 
ation  Factor- porosity  relationships  determined  in  the  laboratory  for 
hree  clays  (kaclinite,  illite,  and  montmorillonite)  ,  Providence  silt, 
our  sands,  :nd  four  marine  sediments  showed  that  porosity  was  predic 
lithin  +_  I  percent.   The  Formation  Factors  ranged  from  1.1  to  5.9  while 
orosity  ranged  from  26  to  95  percent.   The  particle  size  and  distribu- 
tion influenced  the  electrical  resistivity  of  these  sediments  independen- 
f  porosity  while  particle  shape  did  not,   The  laboratory  equipment  was 
conomical,  suzie,    easy  to  operate  and  could  be  used  to  determine  permea- 
bility, tortuosity  and  void  ratio-log  consolidation  pressure  of  sedi- 
ents  . 

The  Electrical  Resistivity  Measuring  System  designed  for  the  Deep 
n  Sediment  Probe  consists  of  three  interchangeable  electrode  arrays 
he  electronic  circuit  and  the  FM  telemetry  data  link  with  6000  fe< 
.oaxial  cable.   The  predicted  porosity  values  obtained  with  the  System 
jn  the  iaboratorv  were  within  +  2  percent  of  the  line  of  best  fit  obt 

teste* 
ett  b?.y    and  most  of  the  data  agreed  well        e  value  obtained  later 
the  laboratory  equipment  on    the  same  se<  i  nt . 
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Unclassified 

The  inner  corer  ring  electrode  array  read  continuous 
conductivity  through  the  water/sediment  column  thus  per- 
mitting a  value  of  salinity  to  be  obtained.   The  bottom 
water  salinity  obtained  in  this  way  was  in  good  agreement 
with  the  interstitial  water  salinity  of  the  cores.   The 
system  detected  a  core  loss,  touch  down,  and  sample  dis- 
turbance.  This,  relatively  low  cost  system  allows  accurate, 
quick  in  situ  porosity  determination  and  shows  promise  in 
other  areas  such  as  the  monitoring  of  pollution,  the  pre- 
diction of  other  bottom  sediment  properties  and  as  a  warning 
system  on  ocean  structures. 
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